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Preface 



The area of science to which this book is devoted 
is very young. Its basic ideas were formulated 
as recently as 1957 in a paper by two British 
scientists, S. R. Broadbent and J. M. Hammers- 
ley. During the mid-1950’s Broadbent was work- 
ing at the British Goal Utilization Research 
Association on the design of gas masks for use 
in coal mines. He came across an interesting 
problem and presented it to the mathematician 
Hammersley. 

The principal element of a mask is filled with 
carbon granules through which the gas must 
flow. Carbon contains pores that are connected 
together in an intricate manner to form a sort of 
complicated maze. A gas can enter the pores by 
being adsorbed on their inner surface. It was 
found that if the pores are wide and well con- 
nected, the gas penetrates deep into the carbon 
filter. Otherwise the gas cannot get beyond the 
outer surface of the carbon. The motion of a gas 
through the maze is a new type of process, which 
differs from diffusion. 

Broadbent and Hammersley called the phe- 
nomenon percolation , and the theory underlying 
processes of this type is referred to as percolation 

theory . 
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Since Broadbent and Hammersley published 
their pioneer paper 28 years ago, it was discovered 
that percolation theory can be used to interpret 
an exceptionally wide variety of physical and 
chemical phenomena. The electrical properties 
of disordered systems, such as amorphous semi- 
conductors, crystalline semiconductors with im- 
purities, and materials formed as mixtures of 
a dielectric and a metal, are probably the best 
understood application of percolation theory. 

The phenomena best described by percolation 
theory are critical phenomena. They are charac- 
terized by a critical point at which some of the prop- 
erties of the system undergo abrupt changes. 
Critical phenomena include sfecond-order phase 
transitions (e.g. the transition of a metal from 
its normal to its superconducting phase when its 
temperature^fis lowered). The physics of all 
critical phenomena is very unusual, but there 
are some common features, the most important 
of which is that in the neighborhood of the critical 
point, the system appears to break into blocks 
which differ in their properties, with the size 
of the individual blocks growing until the system 
approaches the critical point. The blocks are 
quite randomly shaped. In some phenomena the 
whole configuration changes chaotically because 
of thermal motion, while in other phenomena the 
configuration may be frozen in time but changes 
from specimen to specimen. The blocks are in 
complete disorder, so that no regularity is dis- 
cernable from an instantaneous photograph. How- 
ever, this geometry, which can be called the 
geometry of disorder , has quite definite properties 
“on the average”. 
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Actually, geometry is inseparable from physical 
properties. For instance, the physical properties 
of a crystal are determined by the geometry of 
its lattice. Likewise, the “geometry of disorder” 
determines a number of properties of a system 
in the vicinity of a critical point. The most 
interesting feature is that owing to the large size 
of the blocks the geometry is virtually independ- 
ent of the atomic structure of the material and 
thus possesses properties common to a number 
of quite dissimilar systems; hence, the univer- 
sality of the physical properties that we find in 
the neighborhood of critical points. 

This type of relation between physics and geo- 
metry can be traced in percolation theory, and 
in fact, making this relationship clear is the 
main objective of this book. Percolation theory 
is formulated in terms of simple geometric images, 
such as wire nets, spheres or crystal lattices. 
The theory does not operate with the concept of 
temperature, and this makes it possible to clarify 
the idea of critical phenomena for readers not 
familiar with statistical physics. 

Percolation theory, as a theory of critical 
phenomena, is not yet a mathematically rigorous 
science. A large number of important proposi- 
tions have not yet been proved, and certain 
questions have not been answered. However, if 
a rigorous proof exists but appears to be com- 
plicated, I decided to replace it here by argu- 
ments that will explain the result rather than 
prove it. However, an effort was made to separate 
those propositions that have not yet been proved 
from those that have. 

The book offers a detailed presentation of the 
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Preface 

theory of percolation and its various applications. 
However, a definition of what is understood by 
percolation theory and what processes it de- 
scribes is deferred until the last page of the book. 
The definition must encompass so many complex 
concepts that to give it at the beginning would 
be meaningless. Nearly every chapter concen- 
trates on one specific problem whose analysis 
leads to a problem in percolation theory. It is 
assumed that, having read several chapters, the 
readers will become aware of the common ground 
the various percolation problems have, and also 
become aware of the relevance of the title of the 
book. 

As a rule, the problems I chose represent im- 
portant applications of percolation theory. How- 
ever, some (such as laying out an orchard in 
Chapter 5 and the propagation of rumors in 
Chapter 11) are illustrative or even slightly 
humorous. 

The fundamentals of elementary probability 
theory required for understanding the material 
are also covered. Chapter 1 gives a general notion 
of probability and random variables. Chapter 2 
introduces the rules for the addition and multi- 
plication of probabilities and defines the distri- 
bution function. An “easy reading” version of the 
book is obtained by skipping Chapter 2 and the 
other chapters and sections marked with two 
asterisks, although the reader will not then be 
able to follow the derivation of certain of the 
quantitative results or a few of the exercises. 
Nevertheless, this should not hamper the under- 
standing (even if slightly less comprehensive) 
of the other chapters. 



I consider the exercises in the text important. 
As a rule, the exercises are quite simple, and 
I recommend that the reader do them without 
looking at the answers section beforehand unless 
this is specifically advised. 

B. I. Shklovskii played an important role in 
the creation of this book because together we dis- 
cussed its structure and title, and then he read 
the manuscript. I am extremely grateful to him 
for his contribution. 

I am also grateful to my colleagues L. G. Asla- 
mazov, N. B. Vasilyev, Yu. F. Berkovskaya, and 
M. E. Raikh because they read the manuscript 
and suggested a number of useful changes. 

I am especially grateful to my wife N. I. Efros 
for carrying the heavy burden of preparing the 
manuscript for publication. 

A. Efros 
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Part I 

Site Percolation Problem 



Chapter 1 

Percolation Threshold 



Two Pundits Shred a Wire Mesh 

It is not frequent nowadays for a scientific journal 
to publish a report on experiments done with, 
for example, a piece of ordinary wire mesh bought 
in the nearest hardware store. A paper written 
by two American physicists Watson and Leath, 
published in the Physical Review in 1974, was 
definitely not the first one in the realm of per- 
colation theory; nevertheless, it is with this study 
that we start our story. 

The piece of wire mesh used by Watson and 
Leath was a square containing 137 X 137 = 
18 769 nodes, or sites, with the neighbor sites 
separated by a spacing of 1/4 inch = 6.35 mm. 
The scientists soldered copper electrodes to two 
opposite sides of the square and connected the 
mesh to an electric resistance-measuring circuit 
(Fig. la). Then they started to block some sites 
in the mesh and monitored the electric resistance 
as a function of the fraction of blocked sites. As 
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Fig. 1. Schematic of the experiment of Watson and Leath: 
(a) an initial wire mesh (the number of sites in the figure 
is greatly reduced); (b) a piece of wire mesh with blocked 
?y es (blocked sites are shown by black circles, and non- 
nlocked Bites by open (white) circles); (c) a black site sig- 
nines that the four wires coming to the site are not in 
contact, and a white site signifies intact contacts (no elec- 
tnc current flows through black sites in any direction, 
while it flows through white sites in any direction). 
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shown in Fig. \b, c, a site was blocked by simply 
Cll “ m ", a11 four wires joining at this site. 

1 be choice of a new site to be blocked was made 
randomly among those still in the mesh. In 
principle, this could be done by writing the 
coordinates of each site on separate slips of 
paper, putting all these slips into a hat, stirring 
well the contents, and then extracting them one 
by one However such a procedure, as well as 
any o her mechanical method of sequence allot- 
ment, is highly impractical if the number of sites 
is large; consequently, the scientists used a ran- 
dom sequence of site coprdinates generated in 
a computer. Later we shall describe how a com- 
puter can be “made to generate” random numbers, 
but for the time being we can think that a com- 
puter was replaced, with no loss to clarity by 
a hat-and-paper-slips combination. 

Obviously, the electric conductance of the 
mesh decreased with increasing number of blocked 
sites. (Electric conductance is defined as 
a quantity reciprocal of resistance. Resistance 
is measured in ohms, and electric conductance 
in inverse ohms (ohm- 1 ).) Furthermore, if x 
denotes the ratio of nonblocked sites to the total 
number of sites (137 2 ), the electric conductance 
vanished at a certain value of x that hereafter 
we call the threshold (critical) value, or percolation 
threshold, and denote by z c . This vanishing 
occurred when the last path connecting the left- 
and right-hand electrodes had been cut. One of 
the objectives of the experiment was to deter- 
mine x c . It was found that x c = 0.59. 

Probably, the first question that must be 
answered is whether the variable x 0 is random, 



irre producible from one experiment to another, 
or x c is quite definite. Let us assume that the 
experiment is repeated with another square piece 
of wire mesh and with a different random sequence 
of blocked sites. It is common sense to expect 
that the configuration of blocked and intact 
sites being quite different at each stage of the 
second experiment in comparison with the first, 
the cutting of the last path connecting the elec- 
trodes must also occur at a different value of x, 
so that the second experiment must give a new x c . 
This is certainly correct. 

The threshold value x c in the experiment we 
are now discussing is a random quantity, or 
variable. Quantities of this sort will be encoun- 
tered throughout the book, and so it is useful 
to answer from the very beginning the following 
question: 

What Is a Random Variable? 

In mathematics a variable is said to be random 
if the values it takes and the frequency at which 
it takes them are known, but what particular 
value the variable will assume in each particular 
case is unknown (and cannot be known in the 
framework of the given mathematical prob- 
lem). 

Here is a classical example of a random vari- 
able: a cube (a die) with numbered faces is being 
thrown onto a table. The numeral read off the 
topmost face is a random variable. Such a vari- 
able is said to be discrete because it assumes only 
certain discrete values (in the chosen example 
these are six numerals: 1, 2, 3, 4, 5, 6). It would 
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be impossible to predict the specific numeral 
obtained in each given trial (i.e. in each throw), 
but the probability of getting a specific number 
(e.g. 4) is predictable. Assume that the number 
of trials was Q, and face 4 appeared in cases. 
The ratio QJQ is called the relative frequency of 
recording a given value of the random variable 
(here, the numeral 4). If the total number of 
trials is not very large, this ratio varies: in 
another series of Q trials the ratio QJQ may be 
substantially different. However, as the number 
of trials Q increases, these fluctuations progres- 
sively diminish. The relative frequency of record- 
ing a given value of the random variable ap- 
proaches the limit which is called the probability 
of this value. 

Let us denote the probability of obtaining 
face 4 by P (4). If the cube is “honest” (not 
“loaded”), that is, if all its faces are equivalent, 
the quantity P (4) is easily predicted. On the 
average, any one of the six faces of the cube will 
turn up topmost the same number of times, so 
that QJQ = QJQ = . . . = 1/6 provided Q is 
large. Therefore, the probabilities of recording 
the six numerals are identical and equal to 1/6. 

We thus find that after a large number of 
throws, chance melts into the background and 
clears the stage for a regularity, namely, the 
symmetry of cube’s faces. 

Mean Value and Variance 

Let us return to the experiment with the wire 
mesh. We concluded that since the experiment 
operated with a random sequence of blocked 
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sites, the critical concentration x c at which the 
current between the left- and right-hand elec- 
trodes is interrupted is also a random variable, 
and it would be impossible to predict the value 
of x c in each specific experiment. 

The theoretical approach to the situation 
could be a study of the “mean” properties of the 
quantity x c , that is, of the properties apparent 
after a sufficiently large number of experiments 
run under identical conditions. These conditions 
are, first, the total number of sites in the wire 
mesh, JfT {JfT = 137 2 in the experiment described 
above), and second, the properties of the generator 
of random numbers that prescribes a random se- 
quence of blocked sites. The requirement that the 
properties of the generator should not change 
from one experiment to another does not mean 
at all that the sequences of blocked sites must be 
identical. (In this case all values of x c would be 
identical!) All we need is to use the same method 
of generating the random sequence of sites to be 
blocked in all experiments (e.g. a hat with slips 
of paper). 

Having conducted Q experiments with a wire 
mesh containing jfT sites, we obtain Q values x t , 
where the subscript i is the number of the exper- 
imental run. For example, x 15 denotes the x c 
obtained in the fifteenth series of measurements. 
The most important of the mean values is the 
arithmetic mean, or the average, x Ql obtained 
by adding up all x t and dividing the sum by 
the number Q of experiments: 

x 0 =- 3+^2+ .. . -1 -xq 

Q 



2 * 



(1) 
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The bar over a and the subscript Q denote averag- 
ing over the results of Q experiments. The quanti- 
ty is still a random variable. If we conduct 
another series of Q experiments under the same 
conditions and calculate the resultant x Q , this 
x Q will be somewhat different. However, the 
larger the number of experiments, Q, within the 
series, the smaller the differences between the 
averages representing different series. The point 
is that random fluctuations of x t cancel each other 
in the long run, so that as Q increases, the arithmetic 
mean x Q tends to a quite definite value independent 
of Q but dependent on the conditions under which 
the experiments were run. This limiting value is 
called the mean value of a random variable. (In 
probability theory this limiting value is also 
called the expectation of a random variable, but 
this term will not be used in this book.) 

The mean value of the percolation threshold 
of a wire mesh made up of JT sites will be denoted 
by X 0 (JT). The quantity x 0 (JT) is not a random 
variable hut a certain quantity. Its dependence 
on JT is a regular property worth thinking 
about. 

An important characteristic of a random vari- 
able x c is also found in the deviations 8 t of the 
values x t from the mean value: 

6 i = *i — *c (JT) (2) 

The deviations 8 t vary from one experiment to 
another, and we want to select a quantity char- 
acterizing the properties of 6* “on the average”. We 
cannot select the arithmetic mean for this quantity 
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because it tends to zero when Q °o- Indeed, 



. ,.+8q *i+*2+ •••+*<? 

Q 



■x c (jfT) 



But the greater Q is, the less the , d f 1 ^ ere . Ilc ®. 
between the first term on the right-hand side 
of this equality and the second term; this proves 
the proposition made above. This result is ob- 
tained because the individual values of x t neces- 
sarily fall on both sides of the mean value, so 
that on the average the deviations cancel out. 

We could also take the arithmetic mean ol the 
nonnegative quantity |8,|; however, the stand- 
ard procedure is to calculate the variance o (jy ) 
which is the arithmetic mean of the squares 
of deviations (for Q -+■ oo) which, of course, are 
also nonnegative quantities: 



^ + 6 2+--'+ 6 Q 



The quantity 5 (JT) = [S 2 (X)l 1/2 is called the 
root-mean-square deviation of a random vanab e, 
or rms deviation . It is the rms deviation o («# ) 
that characterizes the typical deviation °*the 
values Xi from their mean value x c (jy)» Obvi- 
ously, the quantity 8 (JT) is a l so a function <>f 
the total number of sites, JT , in the wire 

Strictly speaking, x c is a discrete quanti y 
because it is obtained by dividing the number 
of nonblocked sites by the total number of sites, 
JT, and therefore, it assumes only those values 
that convert to integers when multiplied by JV . 
Let us denote all possible distinct values of the 
random variable x Q by x 
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By definition, <?i + <?2 + • • • = (?> so that 
the definition of P ( x h ) yields 

P (x,) + P(x 2 )+ ... = l (7) 



The sum of the probabilities of all the values 
that a random variable can assume equals unity. 



Why a Large Wire Mesh? 

It was very simple to calculate, in the problem 
of “honest” cube, the probability for the random 
variable to assume a specific value. The properties 
of the random variable x c are incomparably more 
complex. 

It will be shown at the end of the next chapter 
how to solve the problem for a square network of 
four nodes (sites) (2x2 sites, JT — 4). The 
result of the solution: the random variable 

x c can assume only two values, namely, 1/2 
and 1/4. The first of these is assumed with prob- 
ability P (1/2) = 2/3, and the second with 
probability P (1/4) = 1/3. According to for- 
mulas (5) and (6) (where the sums consist of only 
two terms), x c (4) = 5/12 and 8 (4) = V 2/12. 

With sufficient patience, we can also solve 
the problem with 3x3 sites (JT = 9). Experi- 
ence shows that the required efforts increase 
enormously if the side of the square is augmented 
by only one site. At the same time, the networks 
whose properties are of especially high interest 
contain a very large number of sites (e.g. 10 15 ). 
Such networks can serve as models of films con- 
sisting of atoms. Indeed, the spacings between 
the atoms in condensed materials (liquids and 
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blocked sites differ from one another the less, 
the greater the number of sites in the wire mesh. 

Why should this be so? The point is that a 
sufficiently large wire mesh manifests all or 
nearly all the possible configurations of intact 
and blocked sites. These configurations appear 
to change places in different experimental runs. 
Consequently, the role of chance is the smaller, 
the greater JT is. An infinite net contains an in- 
finitely large number of large subnets, so that 
randomness^ ceases to play any role at all, and 
the value of x c is not a random variable but a certain 
quantity equal to 

x c = lim x c (jf) 
jir -+oo 

It is this limiting value that is in fact defined 
as percolation threshold, and it was for finding 
this threshold that Watson and Leath have con- 
ducted their experiment. Otherwise why should 
they have taken a wire mesh with nearly 19 000 
sites? They could have taken a 2 X 2 wire mesh! 

Now let us formulate the most important 
result of this chapter: 

The concept of a sharply defined percolation 
threshold independent of the choice of a random 
sequence of blocked sites used in the experiment 
holds in an infinitely large system. No sharply 
defined threshold exists in a finite system, but 
there is the so-called critical region with width of 
the order of 6 (JIT) into which the values of x c , 
which are recorded in the majority of experimental 
runs with different random sequences, fall. As the 
size of the system increases, this region contracts 
to a point. 
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Obviously, there can be no end to improving 
the accuracy of x c in the subsequent decimal 
places (see Exercise 4). 



Exercises 

1. Define a discrete random variable a as the 
numeral on the topmost face of a cube after 
a throw. Find the mean value of a. 

2. Define -percolation threshold as a value of x 
at which percolation sets in not from left to 
right, but from top to bottom. Will this change 
the results of individual experimental runs, 
x c (JT), xj Assume the network to be square. 

3. Answer the same question but take into 
account that percolation threshold is defined as 
the minimum value of x at which percolation 
exists both from left to right and from top to 
bottom. 

4. Answer the same question but in the case 
when percolation threshold is defined as the 
maximum value of x at which there is no per- 
colation both from left to right and from top 
to bottom. 

5. Make use of formula (8) and calculate the 
root-mean-sqnare deviation in the conditions 
of the Watson-Leath experiment (JT = 137 2 ). 
What accuracy can be expected if only one ex- 
perimental run was carried out? 

Note. In principle, the result of a single experi- 
ment may differ very much from the mean value 
x c (JT\ However, using the distribution function 
given below for percolation thresholds (see for- 
mula (6) in Chapter 2), it can be proved that 
the probability for the result of a randomly se- 



28 



Part I. Site Percolation Problem 



leaed experiment to fall in the interval from 

c + _ . 1 0 x o {J^) + 8 is roughly 0.7. The 

E tbf 1S ’ the 1 malIer the “ ty P ica l deviation” 
irom tne mean value is. 



Chapter 2** 

Basic Rules 

for Calculating Probabilities. 

Continuous Random Variables 

This book is devoted to the laws governing disorder 

rand Wldely the conce P ts of probability and 
random quantity, or variable. These concepts 

rhJnt been ?f rtIy , : intr °duced in the preceding 
hapter, so that the reader who is not inclined 
to delve into the mathematical side of the prob- 
-V s ^ elcome t0 ski P Chapter 2, as well as 
w/twi subsequent obapters and sections marked 

To follow m aSt , er f i . sks - ^e reader, who wants 
follow the solution of a number of beautiful 
mathematical problems given in the book and to 
3 mor . e Profound picture of percolation 
theory, must know the rules for the addition 

“ «hrcSpS! ,on of prob » bm “ e5 

Events and Their Probabilities 

probabilit y is used not only when 
we deal with numerical values taken up by a 
random variable. Any experiment with a random 
outcome can be discussed in these terms. The 
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distinct results of experiments are referred to 
as events. The relative frequency of an event 13 
delined as the ratio of the number of trials that 
ted to this event , to the total number 
nrobability of an event is defined as the limit 

which the relative frequency of the event approaches 

when the number of trials tends to infinity. 

FxamDle Red, green, and blue balls in iden- 
ticaSheVs ere pit in’ a box. The halls are then 
stirred and one of them is extracted in a grab bag 
fashion. What is the probability of the event con- 
sisting in extracting a red ball? C0 ?.P^ 
the experiment with dies, here events differ not 
in quantity but in quality (the color of balls). 
However, the argument runs along the same line. 
The number of balls being identical ineachcolor, 
a red ball will be chosen in 1/2 l of : all trials Con- 
seuuentlv the sought probability equals 1/d. 
ffi ^probability of drawing a blue ball or a green 

b8 By definition, probability is a quantity that 
varies from zero to unity. An even wi l ^ ze 
probability is, for instance, the extraction of 
a blue ball from a box with only red balls, lhe 
probability of drawing a red ball from the same 
box is unity. An event with probability equal to 
unity is said to be a certain event, not a random 

^The concept of probability plays an enormous 
role in clarifying the regularities m the wor d 
of random processes. A law is qnite often v rtu- 
ally buried under chance. Imagine that you y 
to find a regularity in birthrates of boys and gi 
Lorn the record of one ^ a ^it y R You see a ran- 
dom sequence of the type BGBBBBGBGG.... 
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Addition of Probabilities 

Events are said to be incompatible if they cannot 
be observed in the same trial. For instance Z 

SeTiS ,f ‘ n dra ' ,i " s * , red l «|l « i-iorn! 

a blue ball her 6 eve V t consisting in drawing 
nrnhlom 4 becau , se > by the conditions of the 

either re’dZ/ h/ 16 bal1 Can be drawn in any trial: 
ing in getting the 6 ’ ° r ^ re , en ' ^he events consist- 
of our rehe numerals 5 and 2 in one throw 
w e are incompatible. 

prlabim” Pr ° VI! tW ° imp ° r,ant properties of 

1. Addition rule. The probability for one (no 

Tvfnts e l Zuau\h eVent 4 ° f several incompatible 
events equals the sum of the probabilities of these 
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events . Suppose that we want to find the proba- 
bility for a cube to show either 3 or 4. The number 
of trials in which these two outcomes were 
obtained is equal to the number of trials that 
gave 3 and the number of trials that gave 4. 
By definition, the sought probability will be 
found if this sum is divided by the total number 
Q of trials, and the limit is found for Q oo. 
The limit of each term of the sum divided by Q 
equals the probability of obtaining one of the 
numerals of interest, the sought probability 
indeed equals the sum of the probabilities of 
each of them. Therefore, the probability of 
obtaining either 3 or 4 is 1/6 + 1/6 = 1/3. The 
probability of obtaining either 1 or 2 or 3 or 4 
is 1/6 + 1/6 + 1/6 4- 1/6 = 2/3. And the prob- 
ability of getting either 1 or 2 or 3 or 4 or 5 or 6 
is 1/6 + 1/6 + 1/6 + 1/6 + 1/6 + 1/6 = 1. This 
result is a particular case of the second 
property. 

2. Let us refer to the set of incompatible events 
that covers all possible outcomes of a given trial 
as to the complete system of events. For instance, 
in the experiment with a cube the complete 
system of events is made up of the events con- 
sisting in drawing the numerals 1, 2, 3, 4, 5, 6. 
The second property states: 

The sum of the probabilities of the events forming 
the complete system equals unity. By virtue of the 
first property, this sum equals the probability 
for one of the events forming the complete sys- 
tem to occur. But by the definition of the com- 
plete system, one of the events necessarily does 
take place. (In the example with a cube this 
means that one numeral out of the six possible 
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numerals has to turn up.) An event that cannot 
but occur is a certain event, and its probability ’ 
equals unity. This proves the second property. 

(In the case of a cube it states that the sum of the 
probabilities for the six possible numerals equals 
unity.) ^ 

As far as the probabilities of different values 
that can be assumed by a random variable are 
concerned, this property was formulated as for- 
mula (7) of Chapter 1. 

In some cases the reservation citing the incom- 
patibility of events may prove quite substantial 
when the rule for the addition of probabilities 
is applied. Consider the following example. 

Example. Five marksmen fire simultaneously 
at one target. The marksmen have identical 

f the tar £ et with the probability 

ot l/o. What is the probability of at least one 
marksman hitting the target? 

We must find the probability of at least one 
ot the following five events (no matter which) 
to take place: the target is hit by the first marks- 
man, it is hit by the second one, and so forth. 

A luring thought is to resort to the rule for the 
addition of probabilities. This rule states that 
the probability of one of the marksmen to hit 
the target equals the sum of the probabilities: 

P = 1/3 + 1/3 + 1/3 + 1/3 + 1/3 - 5/3 

The result is blatantly absurd. The sum is greater 
than unity, and this is nonsense. Where did we > 
make a mistake? Recall that the rule for addition 
is formulated only for incompatible events. Is 
it excluded that several marksmen hit the target 
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simultaneously? Of course, they can. This is 
a typical example of compatible events. Hence, 
the rule for addition is inadmissible here. 

To solve the marksmen problem, we have 
to make use of the rule for the multiplication 
of probabilities formulated in the next section. 

Multiplication of Probabilities 

We again throw a cube. The question that has 
to be answered is the following. We throw the 
cube twice and obtain two numbers. What is 
the probability for these numbers to be 6 and 4, 
and precisely in this order: first 6 and then 4? 

The approach to the solution is standard. We 
make Q trials (each consisting of two throws) 
and determine the number of trials giving the 
desired result. First we select all those trials 
in which the first throw gave 0 regardless of the 
result of the second trial. This problem is already 
familiar. All the faces of the cube being identical, 
the number 6 (or any other number from 1 to 6) 
appears on the first throw in 1/6 of the trials, that 
is, at the first stage we select Q 1 = Ql 6 trials. 
(The number Q is assumed to be quite large, so 
that random deviations from Q x are small.) 
Now we want to select those trials in which the 
second throw gave 4. Each number appears on 
the second throw again with equal probability. 
Consequently, the number 4 was found in 1/6 
of the trials. Hence, the number of trials in 
which 4 followed 6 is Q 2 = l/6*l/6*(?, and the 
probability of this event is 

Q 2 /Q = 1/6- 1/6 = 1/36 
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Now we introduce complications. Let a trial 
consist of three throws, and let us find the prob- 
ability of the trial to yield three numbers in 
a predetermined order, for instance, 4, 5, 1 
or 6, 6, 6. Arguing in a similar manner, we shall 
find that the number of trials giving the sought 
result is Q s = 1/6 -Q 2 = 1/6 • 1/6 • 1/6 • <?, and the 
probability of this result is 

QJQ = 1/6- 1/6* 1/6 = 1/216 

Let us analyze another example. Assume that 
one bicycle per every ten thousand bicycles 
manufactured by a bicycle plant has an unsound 
front-wheel axle, and two have unsound rear- 
wheel axles. Therefore, the probability for a 
bicycle taken at random to have an unsound 
front- wheel axle is 1/10 000, and that for the 
rear- wheel axle is 2/10 000. Assume now that 
the front- and rear-wheel axles are manufactured 
in different workshops, so that a defect in one 
of them does not increase or diminish the prob- 
ability^ of a defect in the other. We want to 
find the probability for a bicycle picked at ran- 
dom to have two unsound axles. We have to 
argue as in the preceding cases. Among Q bicycles, 
select those with unsound front-wheel axle. Their 
number is Q ! 10 000. Among these select the 
bicycle with unsound rear-wheel axle. This gives 
(<?/10 000)* (2/10 000). The sought probability is 
(1/10 000). (2/10 000) - 2.10- 8 . 

In both examples we fixed the probabilities 
of several events and wanted to find the proba- 
bility for these events to occur jointly , that is, 
in the same trial. The results obtained can be 
formulated in a general statement: 
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Then the probability for a bicycle selected 
randomly from the total annual output of the 
plant to have two defective axles equals the 
probability of the front-wheel axle to be defec- 
tive, that is, equals 1/10 000 and not 2-10~ 8 . The 
multiplication rule for probabilities is thus 
valid only for independent events. 

The multiplication rule yields a straightforward 
solution of the problem of five marksmen for- 
mulated in the preceding section. Let us recall 
it: five marksmen shoot simultaneously at a 
target, and the probability for each of them to 
hit the target is 1/3. What is probability of at 
least one marksman hitting the target? 

The solution is most easily obtained if we 
calculate the probability for all five marksmen 
missing (let us denote this probability by P 0 ). 
Since the missing or hitting the target by each 
marksman must be regarded as independent 
events, the probability P 0 equals the product of 
the probabilities for each marksman to miss. 
The event consisting in a marksman hitting 
the target and the event consisting in this marks- 
man missing it form a complete system of events. 
The sum of the probabilities of these two events 
is unity. If the probability for a marksman to 
hit the target equals 1/3, the probability for 
him to miss it is 1 — 1/3 = 2/3. The probability 
for all five marksmen to miss is 

P 0 = 2/3- 2/3 -2/3 -2/3 -2/3 = (2/3) 5 

The event consisting in all five marksmen mis- 
sing and the event consisting in at least one of 
them hitting the target form a complete system 
of events. The sum of the probabilities of these 



events equals unity. Therefore, the sought 
piobability P satisfies the equation 



P. 4- P, = 1 



whence 

p = i _ p 9 = 1 - (2/3) 5 « 0.87 



Exercises 

1. Jp' sites of a network consisting of JP sites 
are blocked. What is the probability of a ran- 
domly selected site to be blocked? nonbiocked. 

2. Find the probability for three consecutive 
throws of a cube to give (a) 1, 2, 3 in any order; 
(b) 1, 2, 2 in an arbitrary order. 

3. Workshop A manufactures high-quality 
parts at a probability 0.8, and workshop B 
manufactures them at a probability 0.9. Three 
parts made by workshop A and four parts made 
by workshop B were taken at random. Find the 
probability for all seven parts to be good. 

Percolation Threshold in a 2x2 Network 

The excerpts from probability theory given above 
are quite adequate for an analysis of a percola- 
tion problem in a square network consisting of 
four sites (JT =4). 

Figure 2 gives the schematic of the experiment 
with a 2 X 2 network. The figures numbering 
the four sites are written on separate slips of 
paper, the slips are put in a hat, and the contents 
of the hat are well shaken. Assume now that the 
first extracted piece of paper has a numeral 1 , 
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and the site 1 is blocked (Fig. 2b). (The line of 
argument and the final results do not change 
in the least if the first to be blocked is a site 
with a different number. The point is that all 
sites in a four-sited network occupy equivalent 




Fig. 2. Calculations for a 2 X 2 network: 

(a) an initial network; ( b ) one network site is blocked; 
(c), ( d ) and (e) two sites are blocked. In case (c), current is 
interrupted only after a third site has been blocked, so 
that x 0 = 1/4. In cases (d) and (e), current is interrupted 
after a second site has been blocked, so that x c = 1/2. 
The three cases (c), ( d ) and (e) are equally probable. 



positions.) If the next site to be blocked is site 2 , 
the current will not be cut (Fig. 2c): it will 
flow through the lower wire. After the third 
site is blocked ( 3 or 4 ), the current is certainly 
interrupted, and one should state that the crit- 
ical fraction of nonblocked sites is 1/4. If, how- 
ever, the second to be blocked is site 3 or 4 , the 
current is interrupted, and the critical fraction 
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. . linfl t0 be 1/2 (Fig. 2d, e). The percolation 
threshold x 0 is therefore a fiscal- 

SKhe probability P for ttt jariablo to aa- 
’"Saltut'Sort' ‘Sich site wM be Uj 

affsa, afsistd * p 

probability 8 of site 3 or 4 being the second. After 
site 1 has been blocked, all three remammg sites 
have equal probabilities of being Moeked at the 
riovt step The sum of the three probabilities 
Saif unity because there three evente form a 
complete system. Hence, each of these pro 

bl The S probability for site 2 to be tiocked^next 
is thus 1/3. But if site 2 is the next then x c - U4. 

Hence, the probability of x c = 1/4 ’ 

is P (1/41 = 1/3- Now we want to find the prob 
ability for either site 3 or 4 to flocked at 

the second step. By virtue of tlie add /^ abi V 
for probabilities, it equals the sum of probabU 
ities l/3 + 1/3 = 2/3. This is precisely the 
probability for a; c to take on the value 1/ , 
therefore, P (1/2) - 2/3. We must have 

P (1/2) + P (1/4) = 1 

because only two values of z c are possible. Indeed, 

P (1/2) + P (1/4) = 2/3 -|- 1/3 = 1 

The mean value of the percolation threshold 
ine mean j;„ tQ ui ft According to for- 
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mula (5) of Chapter 1, 

x c (4) = 1/2 -P (1/2) + 1/4 -P (1/4) 

= 1/2 -2/3 + 1/4 -1/3 = 5/12 

This figure differs quite significantly from the 
threshold value x c = lim x c (,;/" ) that was found 

JT-+ OO 

to be, as we have mentioned already, roughly 
0.59. 

There are no difficulties in calculating the 
variance of percolation threshold. According 
to formula (6) of Chapter 1, 

6 2 (4) = (1/2 - 5/12) 2 • 2/3 + (1/4 - 5/12) 2 -1/3 

= 1/72 

The root-mean-squarc deviation is 
8(4) = 1^ 2/12 



Exercise 

4. Retrace the argument, assuming that the 
Prst site to be blocked was site 3. 



Continuous Random Variables 

So far we have been discussing only discrete ran- 
dom variables. However, there are also continuous 
random variables that can assume any value with- 
in a certain interval of the numerical axis. 

Assume that a random variable a is allowed 
to take on any value y in the range from A to B 
(A sg: y < B) f but some of these values occur 
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frequently, while others only rarely. In order 
to describe this in terms of mathematics, the 
distribution function f ( y ) of a random variable a 
is introduced. 

The main property of distribution iunction is 
as follows: if the points A t and B t fall within the 
interval (A, B) and A, < B x , the probability 
for the value of the random variable to be within 
the interval A, < y < Z?i is equal to the area 




Fig. 3. 



encompassed by the graph of the function f (y), the 
abscissa axis, and the perpendiculars erected at 
points A x and B i (this area is shaded in Fig. 3 
The reader familiar with integral calculus will 
recognize that this probability (let us denote 
• i i 71 / A D \\ ic rvlvDll T» v 'tli ft formula 



B i 

P (A„ B ,) = J / ( y ) dy 

Ai 

Since all the values of the random variable 
are within the interval {A, B) and one of these 
values is necessarily assumed by the variable, 
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the sought area equals unity. In other words, 

B 

P(A, B)=^f(y)dy=l (1) 

A 

Sometimes this equality is called the condition 
of normalization of distribution function. 

A figure whose area is given by the integral (1) 
is called a curvilinear trapezoid (see Fig. 3). 
If the interval (A x , B x ) is so small that the dis- 
tribution function within this interval remains 
practically unchanged, the curvilinear trapezoid 
can be successfully replaced by a rectangle with 
height / (y ± ), where y x is an arbitrary point within 
the interval (A u B x ). Then 

P (A Xl B x ) = / (y x ) A (2) 

where A = B x — A x is the width of the interval 
(A X1 B x ). 

In the mathematical literature the function 
/ (y) is referred to as the probability density . As 
we see from formula (2), if the interval is suffi- 
ciently narrow (otherwise this formula is not 
valid at all!), the probability for the random 
variable to fall within the interval is directly 
proportional to the interval width. The function 
/ (y) is the probability divided by the interval 
width, or it is the probability per unit length 
of the interval, or in other words, it is the prob- 
ability density. Nevertheless, physicists often 
prefer the term “distribution function”. 

Formulas (5) and (6) of Chapter 1 for the mean 
value and variance are rewritten for a continuous 



random variable in the form 



B 

a=* J yf{y) d y 

A 

B _ 

52 — j (i i — a) 2 f (y) dy 

A 



(3) 

(4) 



where a is the mean value of the continuous ran- 
d Tet ,a us a constder an example of distribution 
^“uniform distribution. A contimums random 




y 



Fig. 4. 

with equal probabilities and cannot assume any 
o roftt.ancde with area /„■ 1, where /„ is the vam 
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interval width is 1. The normalization condition 
dictates that /„• 1 =1, that is, f 0 = 1. Therefore, 

1 at 

0 at 1 and y<i 0 ^ 
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by fjr(y)- It will be more convenient to define y 
not as the threshold value as such, but as its 
deviation from the mean value x 0 ( J ’). Then 
fjjriy) A is the probability for the threshold 
value found in a specific experiment to deviate 



Exercise 

5. A continuous random variable a assumes 
with equal probability all values from — 1 to +1. 
Find the probability for it to fall within the 
interval from — 3/4 to — 1/4. 

Percolation Threshold 

as a Continuous Random Variable 

Strictly speaking, percolation threshold is a dis- 
crete random variable because all the values that 
it can assume convert to integers after being 
multiplied by the total number JT of sites. 
However, if JT is very large, the difference be- 
tween the nearest allowed values of this 
random variable is very small (it equals 
t ^ s-1 ). Therefore, in this most important 
case of a very large number of sites, the percola- 
tion threshold x c can be considered with high 
accuracy as a continuous random variable that 
takes on all possible values within a certain 
interval of the numerical axis. Then the quantity 
x c must be characterized by the distribution 
function / ( y ). In this section we describe the 
shape of / (y) at large values of ,jp'\ 

The distribution function for threshold values 
.r c must depend on the number of network sites , 
/, with which we experiment. Therefore, it is 
more correct to denote the distribution function 




Fig. 5. Functions The number of sites, Jf increases 
with increasing curve number. The quantity A is half- 
width of curve 3 at half-height marked by the dashed 
line. 

from the mean value x c (jT) by a quantity that 
falls within a small interval A in the neigh- 
borhood of y. By definition, the mean value 
calculated by using the function fjf(y) and tor “ 
mula (3) equals zero. 

Figure 5 represents the function / ^(y) tor three 
different values of JT • As we see from the graph, 
the distribution function becomes sharper as the 
number of sites, JT , increases. This means that as 
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JT grows, the deviations from the mean value 
(we remind the reader that it was assumed equal 
to zero!) become gradually less probable. Accord- 
ing to formula (1) of the preceding section, the 
areas under all bell-shaped curves must be iden- 
tical. With increasing \J\T the maximum height 
of the curves increases, but the width diminishes. 
The width of a bell-shaped curve can be defined 
as the distance between the points at which this 
curve intersects the horizontal line drawn at a 
distance from the abscissa axis equal to half the 
maximum height of the curve (see Fig. 5). We 
denote this width, usually called half-width, 
by Ajp. 

The probability for the values of percolation 
threshold to fall beyond the half-width of the 
curve is lower than that of the most probable 
value of the threshold by a factor not less than 2. 
Hence, the half-width characterizes a typical 
spread of percolation thresholds, including the 
deviations whose probability is half that at the 
maximum of the curve / 

Remember that the root-mean-square deviation 
carries essentially the same information (Chap- 
ter 1). It does not define a deflection whose prob- 
ability is exactly half as large as the maximum 
probability, but at the same time characterizes 
the typical spread of values of percolation 
thresholds. 

The quantities A jp and 5 jp are proportional to 
each other for any bell-shaped curve, but the 
proportionality factor depends on the shape of 
the curve. Computer calculations demonstrate 
that the percolation threshold distribution func- 



Ch. 2. Basic Rules for Calculating Probabilities 

tion is Gaussian (named after the great mathe- 
matician Karl Friedrich Gauss). This function 
has the following form: 

fjrW-SJ 7S“ P <6> 

where exp a = 6 a , and 6 ^ 2.72 is the base of 
the natural logarithm. It is plotted in Fig. 5 for 
different values of 8^. As the function is sym- 
metric with respect to the point y = 0 at which 
it reaches the maximum, the half-width A jp can 
be found from the relation (see Fig. 5) 

Making use of formula (6), we obtain 
A jr = 2(21n2) 1/2 5 ^ 

By virtue of formula (8) of Chapter 1, the quanti- 
ty 8 jp vanishes as a power function when JT 
oo . This means that as the number of sites 
grows infinitely, the half-width of percolation 
threshold distribution function tends to zero, 
that is, the function itself degenerates to a sharp 
peak. All values of the percolation threshold, 
except one, have zero probability. In this con- 
nection, we repeat the most important state- 
ment of the preceding chapter: as JT oo, per- 
colation threshold converts from a random vari- 
able into a certain quantity. 
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Exercise 

6. (For readers familiar with integral calculus!) > 
Substitute the function f jf(y) defined by for- 
mula (6) into formulas (3) and (4) and prove _ 
that the mean value a calculated by means of 
this function is zero, and the variance 6 2 equals 



Chapter 3 

Infinite Cluster 

This chapter also deals with the site problem of 
percolation theory, but this time we formulate 
it in a different language: the language of clus- 
ters. Furthermore, we discuss a different object: 
instead of a network with blocked sites we dis- 
cuss a doped ferromagnetic, i.e. that with im- 
purity atoms. This is a much more complicated 
object, and so it must be described at least briefly. 

Permanent Magnet 

The reason why iron, nickel, cobalt, and some 
other materials can form permanent magnets 
is probably known to almost everyone. The 
explanation of this phenomenon is that the 
atoms of which such materials are composed are 
themselves elementary magnets. They possess 
magnetic moments . 

The magnetic needle of the compass is a well- 
known system possessing a magnetic moment. 
A magnetic moment is a vector. The needle of 
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the compass has southern and northern poles, 
and its magnetic moment is directed from the 
southern to the northern pole. The external mag- 
netic field makes the compass needle to turn 
so that it is oriented along the magnetic lines 
of force. Of course, any magnetic moment rotates 
in an external field in the same way. The com- 
pass needle produces an external magnetic field. 




Fig. 6. Current loop and its magnetic moment. 

Any magnetic moment produces a completely 
similar magnetic field. 

It was found as early as the beginning of the 
19th century that the source of magnetism is the 
motion of electric charges, that is, the electric 
current. The magnetic moment is produced by 
the current. The magnetic moment p of a planar 
current loop shown in Fig. 6 is given by the for- 
mula 

p = (He) IS 

where I is the electric current, S is the loop area, 
and c is the speed of light (in the CGS system of 
units). The orientation of the vector is perpen- 
dicular to the plane of the loop and such that 



50 



51 



Part I. Site Percolation Problem 

the current flows counterclockwise if we look 
in the direction pointed by the vector’s arrow. 

If the system consists of more than one current 
loop, we can use the heat-to-tail method for the 
addition of the magnetic moments of the loops 
and find the total magnetic moment of the system. 

What is the origin of the atomic magnetic 
moment? Any atom is known to consist of a heavy 
nucleus and an electron shell. The magnetism 
of solids originates precisely from the magnetic 
moment of the shell (the atomic nucleus also 
can have a magnetic moment, but it is approx- 
imately a thousandth of that of the shell). 

First, the magnetic moment of the shell grows 
from the motion of electrons around a heavy 
nucleus. This rotation can be put in correspond- 
ence with a certain current / and an effective 
area S. Besides, quantum mechanics ascribes 
to each electron an additional magnetic moment 
called the spin moment . This last moment is in 
no way related to the characteristics of motion 
of the electron, but represents its inherent prop- 
erty. However, the spin moment produces 
a magnetic field just as ordinary moment does. 
Most often the net magnetic moment of the 
electron shells of the atoms of which a solid 
consists is zero. However, in some materials 
such as iron, nickel, cobalt, and some others 
the electron shells possess a nonzero magnetic 
moment. 

The magnetic moments of neighboring atoms 
in a solid interact with one another. In principle, 
this interaction is similar to the interaction 
between two compass needles placed close togeth- 
er. Each needle produces a magnetic field acting 
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on the other needle. However, the situation 
becomes significantly more complicated because 
the interaction takes place not in the vacuum. 
The outer electron shells of atoms decisively 
affect the character of the interaction, up to 
reversing the directions of the applied forces. 

Experiments show that in some materials the 
interaction between magnetic moments is such 




Fig. 7. Fragment of the crystal lattice of a ferromagnetic. 
The arrows show the orientations of magnetic moments. 

that the forces between them make the moments 
align in the same direction. These materials 
are called ferromagnetics (Fig. 7). 

If the magnetic moments of all atoms are 
oriented in the same direction, the total magnetic 
moment oM equals the arithmetic sum of individual 
moments: M = | \JT , where JT is the number of 
atoms in the solid, and \i is the magnetic moment 
per atom. 

As the size of a body increases, its magnetic 
moment grows proportionally to the volume of 
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the body (the number of atoms, JT , is propor- 
tional to volume). A specific characteristic of 
magnetic properties, that is, the quantity inde- 
pendent of the size and dependent only on the 
properties of the atoms composing the body, is 
the spontaneous magnetization M. It is defined 
as the magnetic moment of a unit volume, that 
is, it is found as the total moment M divided 
by the body volume V : 



where v 0 — V/JT is the volume per atom. 

The word “spontaneous” emphasizes that the 
magnetization M is not induced by an external 
magnetic field but appears because of internal 
forces. A permanent magnet is just such a body 
in which the spontaneous magnetization is non- 
zero. This magnetization produces a magnetic 
field in the medium surrounding the magnet 
(or in the vacuum). 

The unit of spontaneous magnetization in the 
GGS system of units is one gauss (1 G). For in- 
stance, in iron kept at very low temperatures 
M = 1740 G. From this figure we can find the 
magnetic moment p per atom. It is approximately 
2.2 of the spin magnetic moment of the electron. 
The fact that the moment p is found to be of the 
order of the spin moment confirms the correctness 
of our concept of the nature of spontaneous mag- 
netization. 

Thermal motion destroys the magnetic ordering, 
so that there is a critical temperature called 
Curie temperature (or point), above which the 
spontaneous magnetization is zero. For instance, 
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the Curie temperature of iron is 770 °G. Iron 
cannot form permanent magnets at higher tem- 
peratures. 

Doped Ferromagnetics 

Now let us consider a material which is a solid 
solution (mixture) of magnetic and nonmagnetic 
(i.e. having a zero magnetic moment) atoms. 
This is a crystal in whose lattice sites magnetic 
and nonmagnetic atoms sit, their arrangement 
being quite random. 

Let us assume that the interaction between 
the magnetic moments of the atoms decreases 
with distance so fast that we have to take into 
account only the interaction between nearest 
neighbors. This means that if two magnetic 
atoms are at neighboring lattice sites, their 
moments are necessarily parallel, but if they are 
separated by at least one nonmagnetic atom, 
their moments’ orientations are arbitrary: they 
“know nothing” about each other. 

The question that we want to pose now is 
whether the spontaneous magnetization arises 
in the presence of nonmagnetic atoms, and how 
many nonmagnetic atoms are needed to destroy 
the spontaneous magnetization. It will be shown 
below that the answer to this question reduces 
to solving the site percolation problem that 
has been formulated in Chapter 1. 

Let us make some definitions. Two magnetic 
moments will be said to be connected with each 
other if they are neighbors or if they are connected 
via a string of neighboring magnetic atoms (Fig. 8). 
The phrase “are neighbors” signifies that the 
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atoms are nearest neighbors. In a square lattice 
shown in Fig. 8 the nearest neighbors are the 
neighbors along the horizontal and vertical 
directions, but not the neighbors along the dia- 
gonals. An ensemble of connected atoms is said 
to form a cluster. This definition carries the 
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Fig. 8. Fragment of a plane lattice with magnetic (open 
circles) and nonmagnetic (black circles) atoms. The mag- 
netic atoms form one four-atom cluster, one two-atom 
cluster, and five one-atom clusters. The dashed lines de- 
lineate the cluster boundaries. The moments of diiferent 
clusters can have different orientations. 



following meaning. The magnetic interaction 
makes the connected atoms orient their magnetic 
moments in the same direction. As a result, 
each cluster has a resultant magnetic moment 
proportional to the number of atoms of which 
it consists. Moreover, we agreed that magnetic 
atoms which are not nearest neighbors do not 
interact at all. Consequently, the atoms belong- 
ing to different clusters do not interact with 
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one another, and therefore, the relative orienta- 
tion of the magnetic moments of different cluster 

iS £fuTdetTe F by x^the fraction of magnetic 
atoms that is, the ratio of the number of mag- 
netic atoms to the total number of lattice sites. 
By definition, x varies within the interval from 

First we assume that the fraction of magnetic 
atoms is very small (x «1). Obvmusiy in thm 
case they are primarily separated (as raisins in 
a bun). A cluster made of two magnetic atoms 
is a rare event, a cluster of three atoms is sti 
rarer and so on. This statement is extremely 
important for what follows and will be proved 
mathematically a little farther below. So far 
we suggest that those who refuse to recognize 
this statement as obvious should take it fo 

gr At e z«l the number of clusters is therefore 
approximately equal to the number JT of mag- 
netic atoms, and hence, grows proportionally 
to ks the total number of lattice sites increases. 
Ho«veMhe magnetic moments of the* Rustem 
“know nothing” about one another, and hence, 
are randomly oriented with respect to eac 
other (see Fig. 8). In order to find the net mag 
net c moment J of the system, we need to add 
up the magnetic moments of individual atoms 
Zy the head-to-tail method. By virtue of random 
directions, these moments cencel nut so that 
the spontaneous magnetization is ultimately 
zero. \e have thus found that the spontaneous 
magnetization vanishes at low concentrations of 
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Formation of an Infinite Cluster 

Now consider the case in which almost all atoms 
are magnetic. Obviously, a small admixture of 
nonmagnetic atoms does not cancel the spontane- 
ous magnetization but only diminishes it. Let 
us discuss this effect in terms of clusters. At 
x ™ 1 JV* atoms belong to a single cluster. 

666666 

666666 

646*66 

A i A B A A 

o 6 • o • 6 

6 6 6 if 6 6 

A A A A A A 

Fig. 9. Fragment of a plane lattice with magnetic (open 
circles) and nonmagnetic (black circles) atoms at a high 
concentration of magnetic atoms. With the exception of 
atom B, all magnetic atoms belong to a single cluster and 
have identically oriented magnetic moments. 

If x only slightly differs from unity, some atoms 
drop out of this cluster. This occurs, first, because 
some atoms are replaced by nonmagnetic atoms 
(atoms A in Fig. 9), and second, because some 
magnetic atoms form isolated clusters (atom B 
in Fig. 9) with their own orientation of magnetic 
moment. Nevertheless, when x is close to unity, 
a single cluster survives, permeating the whole 
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lattice however large it may be. This cluster is 
said to be infinite. 

Of course, this concept acquires a rigorous 
meaning only for an infinite system. Let us take 
a large series of samples with a given number of 
magnetic atoms, each containing the same total 
number of atoms, and in each sample let us 
find a cluster with the maximum number of mag- 
netic atoms. Now we average the number of 
magnetic atoms belonging to the maximum cluster 
over all samples of the series and denote the 
result of averaging by J\T max * Therefore, ,AT m ax 
is the average number of atoms in the largest 
cluster. The quantity jT max is a function of 
both JT and x. The existence of an infinite cluster 
has the following corollary: at a given value of x 
the ratio J max // tends, when J\T increases 
infinitely, to a nonzero limit 

lim Jte±=P(x) 

JT-* 00 ^ 

The fraction of atoms P ( x ) belonging to the larg- 
est cluster does not depend on the number j\P 
of atoms if JT is sufficiently large but depends 
on x. And as J\T goes to infinity, the quantity 
JT max also tends to infinity. It is for this reason 
that we speak of an infinite cluster. 

Only one infinite cluster can exist in a system. 
Assume that not only the average number of 
atoms in the largest cluster is prescribed at 
the given values of JT and x but also the average 
number of atoms in the next-largest ^cluster. 
Let us denote this last number by JT max* By 
definition /; a x < JT m ax- The statement that 
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only one infinite cluster exists in a system signi- 
fies that 

lim 

JT-+CO 

at any value of x. This means that two clusters 
permeating the system necessarily merge some- 
where and transform into a single cluster.* 
We have thus obtained that at a sufficiently 
high concentration x of magnetic atoms, a certain 




Fig. 10. Graphs of functions P (x) and a (x)/o (1). Although 
the two functions vanish at a common point, at the 
critical value x 0 they behave quite differently, for reasons 
to be discussed in Part III of the book. 

fraction of these atoms belong to one cluster, 
and they all have a common direction of atomic 
magnetic moments. This means that the spontane- 



* Strictly speaking, the statement on the survival of a 
single infinite cluster has not been proved. There are 
arguments supporting the validity of this conjecture to a 
greater or lesser extent, but it would be more correct to 
say that specialists simply take it for granted. 
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ous magnetization is nonzero: 

M = ^-P(x) 

V 0 

Recall now that only small clusters exist at 
a low concentration x of magnetic atoms. In this 
case an increase in the number of sites, J\T , results 
only in the growth of the number of small clusters, 
but not in a larger number of particles in each 
cluster. Then 

lim =0, that is, />(*) = 0 

JT* °° J 

We thus have to conclude that there is a critical 
concentration x c at which an infinite cluster is 
formed, x c satisfying the inequality 0 < x c < 1. 
At this very concentration x c the spontaneous 
magnetization appears, and the function P (a;) 
becomes distinct from zero (Fig. 10). Consequent- 
ly, a material cannot form a permanent magnet if 
the fraction of nonmagnetic atoms is greater than 
1 _ Xc (the fraction of magnetic atoms is less 
than x G ). 

Exercise 

1. Find the function P (x) for x not very differ- 
ent from unity. 

Site Percolation Problem Revisited 

What remains to be done now is to say that from 
the point of view of the critical concentration x c 
the problem of network conduction and the prob- 
lem of doped ferromagnetic are identical. 
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The electric conduction problem can be refor- 
mulated just as easily in cluster terms. We then 
only want to replace the concept “nonmagnetic 
atom” in all definitions by the term “blocked site”. 

Figure 8 shows a configuration of magnetic 
(open circles) and nonmagnetic (black circles) 




Fig. 11. The same configuration as in Fig. 8, but with 
nonblocked sites substituted for magnetic atoms. 

atoms. Let us perform the above-described repla- 
cement in this configuration and convert it 
from a doped ferromagnetic to a wire mesh with 
cut-out knots. This is done by removing in 
Fig. 8 the arrows indicating the directions of 
magnetic moments and then by tracing the wires 
so that they connect the lattice sites (Fig. 11). 

Figure 11 clearly shows the main property of 
clusters applied to the network problem. If 
a potential difference is applied to any pair of 
sites within one cluster, a closed circuit with electric 
current is formed. (Of course, this property is 
meaningful only for clusters containing at least 
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two sites.) When a potential difference is applied 
to any pair of sites belonging to different clusters , 
a closed circuit is not formed and electric current 
is zero. If x < x c , the system contains only clus- 
ters composed of a finite number of sites, so 
that as the size of the system increases, the cur- 
rent between the lateral electrodes will neces- 
sarily be interrupted sooner or later. But if 
x > x G , a very large system cannot but include 
into its lateral faces some sites belonging to an 
infinite cluster. This infinite cluster will create 
the electric conductivity o (x) distinct from zero 
and independent of the size of the system. 

Let us return to Fig. 10 which shows the func- 
tions P (x) (the fraction of the sites belonging 
to an infinite cluster) and a (x)/o (1) (cr (1) is 
the electric conductivity at x = 1, i.e. with 
no blocked sites). Both functions vanish at the 
same point that we first identified as the per- 
colation threshold and later the point at which 
an infinite cluster is born. 

We were thus always dealing with a problem 
of percolation theory that is referred to as the 
site percolation problem. If we were interested 
in the value of x c for a “plane ferromagnetic”, 
then we could look up the result of the experiment 
with the wire mesh and say that x c = 0.59. 
However, actual ferromagnetic materials crystal- 
lize into three-dimensional, not plane, lattices. 
An example of a three-dimensional lattice is the 
primitive cubic lattice whose unit cell has already 
been shown in Fig. 7. 

The problem of wire mesh electric conductance 
is easily generalized to the three-dimensional 
case. Imagine a cube consisting of numerous 
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cells and made by welding together pieces of wire, 
as shown in Fig. 12. We can solder two metal 
plates on the opposite faces of this cube and, as 
in Fig. 1, make an electric circuit and study its 
electric conductance as a function of the number 
of blocked sites. The blocking of each site inter- 
rupts the contact between six wires entering 
this site. As in the two-dimensional case, there 




Fig. 12. Simple cubic lattice. 



is a critical concentration x c of nonblocked sites 
below which the electric conductance vanishes. 

The problem of doped ferromagnetic and the 
related concept of an infinite cluster were equally 
valid for plane and three-dimensional lattices. 
The critical concentration x c of magnetic atoms 
at which an infinite cluster appears is at the same 
time the threshold of percolation from one cube 
face to the opposite face in a sufficiently large 
cube. It must be borne in mind that the quantity 
x c is strongly dependent on the type of the lat- 
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tice. It was found to be 0.59 for a square lattice, 
but x c = 0.31 for a primitive cubic lattice. 
(For details see Chapter 6.) 

Clusters at a Low Concentration 
of Magnetic Atoms ** 

The conclusions drawn in the preceding sections 
were based to a large extent on the statement that 
at a low concentration x of magnetic atoms these 
atoms are usually single, two-atom clusters are 
rare, three-atom clusters are even rarer , and so on. 
Let us prove this statement. 

We introduce a function P<ji{x) giving the 
probability for a randomly chosen atom to belong 
to a cluster consisting of not less than M atoms. 
This means that the atom chosen at random is 
(a) magnetic, and (b) connected to not less than 
o/H — 1 other magnetic atoms. Let us calculate 
the function P^ifo) f° r = i and M = 2. 

The function P x (x) is the probability for a 
randomly chosen atom to be magnetic. This 
probability equals x (see Exercise 1 to Chapter 2 
where the word “nonblocked” must be replaced 
by “magnetic”, and the word “blocked” by ‘non- 
magnetic”): 

P x (x) = x (!) 

The function P 2 {x) equals the probability for 
a randomly chosen atom to be magnetic and to 
have another magnetic atom among its nearest 
neighbors. Obviously, these two events are inde- 
pendent, so that the sought probability can be 
represented by the product of the probabilities 
of these two events. The first of them (the prob- 
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ability for the atom to be magnetic) equals x 
and thus 

P 2 (x) = xW (x) (2) 

where W (x) is the probability for at least one 
magnetic atom to be found among the nearest 
neighbors of the atom. The function W (x) 
depends on what kind of lattice we consider. 
Let us limit the analysis to the square lattice 

O 

o . ; 

3 0 1 

O 

4 

Fig. 13. 

in which each atom has four nearest neighbors 
(Fig. 13). We need to find the probability for 
at least one of the atoms 1, 2, 3, 4 to be mag- 
netic. 

This problem is best solved as follows. The 
event consisting in all four atoms being nonmag- 
netic and the event consisting in at least one 
of the four atoms being magnetic form a complete 
system of events. The sum of the probabilities 
of these two events equals unity. We denote the 
probability of the first event by W 0 ; the probabil- 
ity of the second event is just the sought quantity, 
W. As we have said, W + W 0 = 1. The proba- 
bility for atom 1 to be nonmagnetic equals 1 — x. 
The probability for atom 2 or 3 or 4 to be non- 
magnetic also equals 1 — x. The events consisting 
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in different atoms being nonmagnetic are inde- 
pendent. Consequently, the probability for all 
four atoms to be nonmagnetic equals the product 
of all four probabilities: W 0 = (1 #) 4 . This 

yields W = 1 — W 0 = 1 — (1 — xf . By virtue 
of formula (2), 

P 2 (x) = x [1 - (1 - xf\ (3) 



If x <Cl, the expression for P 2 (x) can be sim- 
plified by dropping the terms with high powers 
of x . Using the binomial formula, we obtain 

1 — (1 — xY = 4# — 6# 2 + 4.z 3 — x 4 



Note that if x <C 1» the ratio of each next term 
to its predecessor is small: 



6j 2 

Ax 






and 




Therefore, if x <§; 1, we can write to a good accu- 
racy that 

1 — (1 — xY ~ 4 x 
whence 

P 2 (x) ~ 4x 2 (4) 

A comparison of formulas (1) and (4) shows that 
at x 1 the ratio 



P 2 (x) 
p iW 



4^<1 



(5) 



that is, the probability that a randomly chosen 
atom belongs to a cluster of two or more atoms 
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is much less than the probability for it to form 
a one-atom cluster. 

A simpler derivation of formula (4) for P 2 (x) 
can be given, immediately taking into account 
the condition a:<l. It is relatively simple with 
this derivation to calculate the functions PqM{x) 
for oM > 2 (see Exercise 3). The derivation is as 
follows: a cluster consisting of more than two 
atoms, one of them being atom 0 , necessarily 
incorporates either atom 1 or 2 or 3 or 4 . The 
probability for atoms 0 and 1 to belong to one 
cluster equals the probability of both of them 
being magnetic, and is given by the product 
of the probabilities for each of these atoms to 
be magnetic, that is, x*x = x 2 . The same can 
be said about the probability for a cluster to be 
formed by atoms 02, 03 or 04. All these pro- 
babilities equal x 2 . The probability for at least 
one of these events to realize equals the sum of 
the probabilities, that is, Ax 2 , hence we arrive 
at formula (4). 

This conclusion holds only for a; <1. Only 
under this condition can we use the rule for the 
addition of probabilities. Indeed, the addition 
rule holds only for incompatible events. But the 
event consisting in atoms 0 and 1 being magnetic 
is compatible with the event of atoms 0 and 2 
being magnetic. The simultaneous realization 
of these events means that all three atoms 0,1, 
and 2 are magnetic, and hence, form a three-atom 
cluster. The probability of the compound event 
equals the product of the probabilities for all 
three atoms to be magnetic, that is, equals 
X*X'X = x 3 . 

If x <Cl, this probability is much less than 
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the calculated probability of the formation of 
a two-atom cluster. Therefore, the probability 
for the events to occur jointly can be neglected, 
and the events can be treated as incompatible. 
This justifies the derivation given above, under 
the condition x <<1. 

In fact this means that if x 1 , the probability of 

the formation of a three-atom cluster can be neglected 
in calculating P 2 (x). 

Therefore, in the case of a: <1 the function 
P 2 (x) actually coincides with the probability 
for a site chosen at random to belong to a cluster 
consisting of two (and not more than two!) 
atoms. Correspondingly, the function P 3 (x) 
describes a three-atom cluster. It is proportional 
to x 3 and is small in comparison with P 2 (#). 
The general result is that the function P ^(x) 

contains the powers of x not smaller than x^ , 
and that for x -<1 we should have P ^(x) <C 

^ ( x )• 

We thus found that if a site chosen at random 
in the case of x <C 1 is magnetic, it almost cer- 
tainly will form a one-site cluster. The probabil- 
ity for it to belong to an M - site cluster sharply 
drops with increasing <M. 

Exercises 

2. Find P 2 {x) for the primitive cubic lattice 
shown in Fig. 12. Find it for an arbitrary lattice 
in which each atom has z nearest neighbors. 

3. Find P s (#) for a square lattice, making use 
of the condition x <^1. 

4. Find P z (x) for a square lattice, without 
employing the condition 
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Chapter 4 

Solution 

of the Site Percolation Problem 
by Monte Carlo 
Computer Techniques 

The Monte Carlo method is one of the most wide- 
spread methods of solving percolation theory 
problems. The aim of this chapter is to give 
a general idea of this method, to explain in detail 
how the main element of the method — the gene- 
rator of random numbers — operates, and to give, 
in conclusion, a concrete computer program that 
makes it possible to determine the percolation 
threshold of the site problem. The very first 
question that arises now is: 

Why Monte Carlo? 

“But what is zero? You see that croupier, the 
curly-headed one, the chief one, showed zero now? 
And why did he scoop up everything that was on 
the table? Such a heap, he took it all for himself. 
What is the meaning of it?” 

“Zero, Granny, means that the bank wins all. 
If the little ball falls on zero, everything on the 
table goes to the bank....” 

“You don’t say so! And I shall get nothing?” 
“No, Granny, if before this you had staked on 
zero you would have got thirty-five times what 
you staked.” 

“What! Thirty-five times, and does it often 
turn up? Why don’t they stake on it, the fools?” 
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“There are thirty-six chances against it, Gran- 
v » 

“What nonsense. Potapitchl, Potapitch! Stay, 
Pve money with me here. 

She took out of her pocket a tightly packed 
purse, and picked out of it a friedricli d or. 

“Stake it on the zero at once. „ 

“Granny, zero has only just turned up, l 
■said, “so now it won’t turn up for a long time. 
You will lose a great deal, wait a little, anyway. 
“Oh, nonsense; put it down!” 

“As you please, but it may not turn up again 
till the evening. You may go on staking thous- 
ands; it has happened.” 

“Oh, nonsense, nonsense. If you are afraid ol 
the wolf you shouldn’t go into the forest. What? 
Have I lost? Stake again!”* 

This excerpt from Dostoevsky s Gamblei 
describes the most exciting game °f chance ot 
the last century: the roulette. It should be noted 
that from the standpoint of probability theory, 
the inexperienced ecstatic grandmother reveals 
more common sense than the Gambler w 10 is er 
consultant. The probability of getting zero will 
not diminish in the least if it turned up on le 
preceding round. There is no sense in waiting, as 
the Gambler advises. This frequently encounter- 
ed misconception is presumably based on the 
misunderstanding of the fact that the probabil- 
ity of getting zero two times running^ is small. 
But this does not mean at all that if zero turned 
up once, the probability to get it the second time 

* “Gambler” in The Short Novels of Dostoevsky , Dial 
Press, 1945, New York; translated from the Russian by 
Constance Garnett, 
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on the next throw is less than on the first. Of 
course, the probability remains absolutely the 
same. The city of Monte Carlo in the Monaco 
principality earned its fame as the world capital 
of roulette. It was this city that gave its name 
to one of the most powerful among the modern 
numerical methods in mathematics. 

And what is it that this method and the roulette 
share? It is the fact that the main element of the 
Monte Carlo method is that very revolving 
wheel that decides the fates of people, destroying 
some and rewarding others in numerous casinos 
oi Monte Carlo. Actually, mathematicians have 
greatly improved it. There is no revolving wheel, 
but^a standard computer program that is called 
the random-number generator”. But this does not 
alter the principal point. From the mathematical 
viewpoint the wheel of the roulette game is 
nothing but a random-number generator. 

What Is the Monte Carlo Method? 

As a rule, the term “Monte Carlo method” is 
applied to any mathematical technique essent- 
ially based on a random-number generator. 

Usually a modern computer has a standard 
program that generates random numbers dist- 
ributed randomly within the interval from zero to 
unity, that is, it “plays out” the values of a con- 
tinuous random variable that assumes with equal 
probability all values within the interval (0, 1). 

Each time the program is addressed, it out- 
puts one such number with a predetermined 
number of decimal places that depends on the 
computer model. 
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The simplest use of the Monte Carlo, method 
is for calculating integrals. Imagine, for instance, 
that we need to calculate the volume within a 
closed surface of a complicated shape. Let us 
choose a cube that we know for certain to in- 
clude the whole surface (Fig. 14). Now we obtain 




Fig. 14. To the Monte Carlo determination of a pear’s 
volume. 

from the random-number generator a set of points 
distributed uniformly within the cube. This is 
done as follows. Assume that the cube edge is 
L long, and all the three coordinates of the points 
within the cube vary from zero to L (see Fig. 14). 
Addressing three times the random-number gen- 
erator, we obtain three numbers i/j, j/ 2 , Us within 
the interval (0, 1). These numbers will give us 
the coordinates of the first point within the cube 
by the formulas X , = Ly lf Y t = Ly 2 , Z x — /d/a- 
Repeating this procedure Q times, we obtain y 
points that on the average uniformly fill the 
cube. Let Q x be the number of points that hap- 
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Jlktr b t f j U W . 1 , thin 1 the surface. The points being 
debuted nnifomly, the number Q, character- 
lzes the volume enclosed by the surface. Name- 

equals L 3 (? 1 /$ ffiCientIy large ’ the S0llgIlt volume 

t A t a a °f y *i is available that makes it possible 

obtain* the H nu “ ber of Points, Q, required to 
obtain a result with a prescribed accuracy. An 

pineal approach can be used: we repeat the 
Sw2? ei ? Se 7 raI times ’ each time using a 

tb7rp«!lt f numbers ’ and then compare 

the results. If the results fluctuate within the 

prescribed accuracy, everything is fine and the 
answer is correct. The Monte Carlo technique 
has very important advantages over conventional 
integration methods in the case of multidimen- 
sional spaces (multiple integrals). 

are aases , in which the Monte Carlo 
method is the only one practicable. Imagine 

wan ,* t0 study behavior of a system 
consisting °f an enormous number of particles, 

f “k? 16 ’ the behavior of a ? as - In Principle 
tins problem must be solved by the methods of 

P . hy 'l 1CS; however, these methods are 
found to be ineffective if the interaction between 

? SI 6 ' j ® trong ( this happens when densitv 
; and temperature low). Then the proper- 
t es of the gas are studied by computer simula- 
ion. Tlie nurab er of gas particles participating 
in the simulation is determined by the size of 

f n °f“ PUt t- memo f. y ‘ The memory must store the 
ntormation on the coordinates of all particles. 

I he simulation consists in randomly picking 
up one of the particles which then travels a ran- 
dom distance. (This means that the coordinates 



of this particle in the computer memory change.) 
Then another particle is chosen randomly, and 
so on. The potential energy of the interaction 
between the particles in the gas depends on their 
mutual arrangement. The energy is calculated 
from the very beginning and then is recalculated 
after each displacement. The probabilities for 
particles to be displaced by a specific distance 
are chosen to agree with the potential energy 
in such a way that the model system “live” on 
the average as the real system does. 

As a result, the computer memory stores what 
resembles “instantaneous photographs” of the 
gas taken at successive time moments. The 
photographs include the coordinates of all gas 
particles and thus make possible the calcula- 
tion of mean thermodynamic characteristics, such 
as pressure, heat capacity, and some others. 

The very procedure of simulation much re- 
sembles a game which is ruled by a rigid code 
that includes calling the roulette, that is, the 
random-number generator. Slight deviations from 
the rules, or a “crooked” roulette, make some 
configurations of atoms in the gas appear more 
often than others; this affects the result of aver- 
aging and invalidates the answer. 

Random-number generators are used not only 
in the Monte Carlo method hut also in the so- 
called analogue experiments exemplified by the 
experiment with a wire mesh described in Chapter 
1. We have mentioned earlier that the random 
sequence of blocked sites required for this ex- 
periment was generated on a computer. In order 
to choose the next site, it is necessary to address 
the program and obtain a random number y. 
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We must multiply it by the total number JT of 
sites and then add unity to the product. Then 
we have to take the integral part of JT y + 1 . 
This gives an integer that falls in the prescribed 
interval from 1 tOe^. In fact, such numbers may 
appear more than once. But there is no harm in 
it. If it is found that the site with the number 
just generated was blocked earlier, it is neces- 
sary to address the computer for a new random 
number and convert it into the number of the 
site. 

Somewhat later we shall describe a program 
with which we can calculate percolation threshold 
via the Monte Carlo method, but now we shall 
discuss the principal element of this method, 
namely, the random-number generator. 

How to Think Up a Random Number 

We thus need random numbers distributed uni- 
formly in the range from zero to unity. So far 
the problem is not yet rigorously formulated. 
We need to know how many decimal places are 
required in each number. Let us assume that 
only two decimal places are necessary. Then take 
a telephone directory, open it at an arbitrary 
page, and copy a sequence of the last two figures 
of each telephone number, preceding these numer- 
als with “0.”. This will give an acceptable two- 
digit table of random numbers. And what are 
we to do if ten-digit random numbers are required? 
Looks like a computer is then unavoidable. 

Come to think of it, the very idea that a com- 
puter can generate random numbers may appear 
strange. Indeed, any computer follows the algo- 
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rithm that was fed into it, that is, it performs pre- 
cisely the actions programmed by a man. How 
then can a chance element be introduced into 
its performance? 

In fact, there is no chance element in the pro- 
gram of the random-number generator . The princi- 
ple of its functioning isjas follows. When the pro- 
gram is addressed for the first time, some num- 
ber y 0 must be fixed. Then a rigorously defined 
sequence of operations transforms this number 
into a new number 

Ui = O/o) (!) 

where O is a specially chosen function or a se- 
quence of operations that transform y 0 into y v 
It is this function that determines the algorithm 
of random-number generation. In its turn, the 
number y x is the starting point for generating 
the next number y 2 by the same recipe: 

Vt = ® O/i) (2) 

Obviously, the function O is such that all 
numbers y ly y 2 , . . ., y n satisfy the inequalities 
O^z/n^l- This sequence of numbers repre- 
sents the sought sequence of random numbers. 

It can be readily ascertained that the thus ge- 
nerated sequence cannot be infinite. Indeed, 
any computer operates only with numbers con- 
taining a limited number of decimal places. 

Their number is limited. (There are only 10 2 
two-digit numbers and 10 n rc-digit numbers.) 

Therefore, sooner or later a number y n will 
coincide with an earlier generated number, say, 
y n _ L . Then the sequence will be repeated: y n +i 
will coincide with y n ~h+x% so on* 
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Consequently , the sequence generated by formulas 

(1) and (2) is inevitably periodic . In view of 
this fact, such numbers are not truly random and 
are referred to as pseudorandom (i.e. “as if ran- 
dom”, or resembling random numbers). 

Nevertheless, they can be used as random 
numbers provided the sequence of numbers re- 
quired to solve a given problem is shorter than 
the period L of the sequence. 

In its turn, the period L is determined by the 
number of decimal places with which the com- 
puter operates (i.e. the number of memory cells 
allotted to each number) and by the quality 
of the algorithm (i.e. by the properties of the 
function O that enters into formulas (1) and 

(2) ). 

The development of a good random-number 
generator is a very difficult problem. As a rule, 
“off-the-cuff” generators prove to be of poor qual- 
ity. Some specific generators are discussed below. 

The Mid-Square Method 

Historically, this was the first ever method of 
random-number generation by a computer. It 
was suggested in 1946 by the brilliant mathe- 
matician John von Neumann. The method makes 
it possible to generate random numbers with 
an arbitrary number of decimal places, corre- 
sponding to the capabilities of the computer. 
The method is extremely simple. Assume that 
we need four-digit numbers. Let us choose the 
first number X 0 arbitrarily. For instance, 
X 0 — 8219. Raise it to the second power. This 
gives an eight-digit number 67551961. Extract 
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the four middle digits: 5519. The second number 
of the sequence is X x = 5519. Now raise 5519 to 
the second power, obtaining 30 45936 1. I he 
! third random number is then X 2 = 4593. If 
the first digits among the middle group are zeros, 
the resultant number has a smaller number ol 
•nonzero digits. For instance, XI — 21 09564 J, 
X 3 = 956. Squaring X 2 3 , we can obtain an eight- 
digit number by adding zeros on the left, 
XI = 00913936, so that X 4 = 9139, and so 

*° Random numbers y n distributed uniformly 
within the interval from zero to unity are ob- 
tained from the numbers X n via the formula 
u n = XJ 10 4 , where n — 0, 1, 2, 3, . . ., so that 
y 0 = 0.8219, y, = 0.5519, y 2 = 0.4593, and so 

At first glance, the method looks attractive. 
However, a careful investigation demonstrated 
that this is not true. The main shortcoming of 
the method is that some starting numbers make 
the sequence “go into a cycle”. For instance, 
it was found that in the four-digit class of num- 
bers the sequences often terminate by a cycle 
6100, 2100, 4100, 8100, 6100. The period of this 
cycle is mere 4, and this is obviously unaccept- 

There is even a number that immediately 
reproduces itself. This is 3792 (3792 2 = 14379264). 
Zero also reproduces itself, and quite often se- 
quences generated by the mid-square method 
degenerate to zero. Consequently, nowadays 
the mid-square method is only of historical in- 
terest. 
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Exercises 

1. Compose a sequence of four-digit numbers, 
starting with 0085, 0067, 0032. Show that all i 
three sequences are monotone decreasing (each 
subsequent number is smaller than the preceding 
number) and rather rapidly degenerate to zero. 

2. Now prove that this constitutes the general 
shortcoming of the mid-square method: if 272 - 
digit numbers are used and a number b appear- 
ing in the sequence has zeros for n leading digits, 
the sequence becomes monotone decreasing and 
finally degenerates to zero. 

Linear Congruent Method 

At the present moment this method of generat- 
ing random numbers is regarded as the best. 

The idea is as follows. Four integers are chosen: 

(1) multiplier k, 

(2) shift c, 

(3) modulus /re, 

(4) the first number of the sequence X 0 . 

The sequence of random numbers is deter- 
mined by the formula 

X n+l = (kX n + c)modm ( 3 ) 

where the subscript n runs through 0, 1, 2, ... 

The symbol b mod m denotes the remainder after 
dividing b by m. For instance, 

b 25 6 30 3 147 

m 10 10 10 12 12 

b mod to 5 6 0 3 3 



Obviously, b mod m < m. Consequently, all 
numbers of the sequence, X n , satisfy the inequal- 
ity X n < 722 . The sequence of the numbers y n 
distributed uniformly within the interval from 
zero to unity is obtained by the formula 

y n — -jH' * n — 0> 1j . 2j . . . (4) 

It is not just any choice of the four starting 
numbers that leads to good results. Note first 
of all that the sequence X n must necessarily be 
periodic, and the period cannot be greater than 
772 . Indeed, all X n being integers, with X n < m, 
the number of different numbers cannot exceed m. 
For this reason, a number that has already oc- 
curred in the sequence will appear at least be- 
ginning with 72 = 7?2, and the sequence will re- 
peat itself. 

However, it is far from simple to generate a 
sequence with a maximum possible period L = m. 
If the starting numbers are not carefully selected, 
the generated sequences will have, as a rule, 
short periods. 

Exercises 

3. Write the sequence of numbers X n generated 
by means of formula (3) for 7c = 3,c=0,X 0 = 5, 
772 = 20. 

4. Write the sequence of numbers X n generated 
by means of formula (3) for k = 3, c = 1, X 0 — 
= 5, 772 = 20. 

5. Write the sequence of numbers X n generated 
by means of formula (3) for k — = 3, c = 2,X 0 = 5, 
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m == 20. Make sure that in all three cases the 
period of the sequences is essentially shorter than 
20. Analyze other examples. 

The following theorem is true. If a sequence is > 
generated by means of formula (3) for c ^ 0, its 
period equals m if and only if the following condi- 
tions are met : 

(i) c and m are coprime numbers (have no other 
common divisor than 1); 

(ii) b = k — 1 is a multiple of p for any prime p 
which is a divisor of m\ 

(iii) b is a multiple of 4 if m is a multiple of 4. 
Unfortunately, the proof of this theorem is too 
complicated to be given here. 

6. Make sure that the conditions imposed by 
the theorem given above did not hold in all the 
examples presented in Exercises 3-5. 

7. Make sure that the set of integers k = 11, 
c = 3, m = 5 meets the conditions of the above 
theorem and yields the period L = 5 for arbi- 
trary X 0 . 

Therefore, a generator with a maximum possi- 
ble period L will be obtained if we take for m the 
largest number with which a given computer can 
operate, and choose the other numbers in accord- 
ance with the theorem given above. 

However, the period is not the only character- 
istic of the quality of a random sequence. For 
instance, let us consider a sequence generated 
with k = c = 1. The sequence is 0, 1, 2, 3, . . ., 
m — 1, 0, 1, 2, 3, . . ., m — 1, 0, . . .. Its 
period is indeed m, but as a random sequence it 
is absolutely unacceptable. 

A complicated system of tests has been elabo- 
rated to determine the quality of a random- 
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number generator. Therefore, only reliable ge- 
nerators can be advised for applications. 

When a random-number generator is being 
chosen, the properties of the computer are im- 
portant not only from the standpoint of the choice 
of the period of maximum possible length. The 
rate at which random numbers are generated 
also depends on the choice of the starting numbers. 
And it is found that different generators prove 
to be optimal for different types of computers. 

Programs involving Monte Carlo computa- 
tions often have to address a random-number ge- 
nerator an enormous number of times (tens and 
hundreds of millions of times). Consequently, 
high speed is one of the most important character- 
istics of a generator. 

The generator recommended for a BESM-6 
computer has k = 5 17 , c = 0, m == 2 40 , and 
odd X 0 . This set of numbers does not meet the 
conditions of the above theorem (c =0), and 
the period of this generator is less than m. How- 
ever, another theorem has been proved for gene- 
rators with c = 0, and by virtue of this theorem, 
the period of the recommended generator is 
2 38 « 2.75- 10 11 . 

Determination of Percolation Threshold 
by Monte Carlo Simulation on a Computer. 
Distribution of Blocked 
and Nonblocked Sites 

Now we shall describe in detail a computer pro- 
gram that determines the percolation threshold 
by Monte Carlo techniques. Note that this pro- 
gram is not unique. In fact, each group of re- 
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searchers that deals with ‘^ problems prc- 
fors to work with their own program that lias some 
individual features. This 

features of different computers and to some ex 

tent from the experience accumulated by ma 

V1 The 1 problem that we mean is the site P roblei ^’ 
and for the sake of simplicity, we consider only 
a two-dimensional square lattice. Actually, 
till become clear later that a generalization of 
the method to an arbitrary lattice of arbitrary 
dimensionality is readily obtainable. 

Let us look at percolation in a square with the 
sidelontaimng L sites, so that the total number 
of sites is JT = L\ We assume the distance be- 
tween the sites to be unity, and 
sites by their coordinates X and Y. For ^ ^ 
a site with coordinates X — Y 
site in the ninth column from the left and m 

Tst'udy percolation, we must h* 
which°of e the° sites IJ hleehed and *ch are 



KWt-orl sites in order to exceed the percolation 
threshold. To achieve this, we first assign a spec,- 
fic number V to each site. Each site bei g 
terized bv two coordinates X and Y , this is q 
alent to introducing a function of two variables 
V (X Y) whose arguments X and Y do no 
thiough all possible values but are a 10^0 
assume only integral values within t ^e 
from 1 to L. Programmers can such a funcUo 

a two-dimensional array, and the values e^ 

on by this function are said to be ^ T7 /Q/| Q7N 



nnsK: :■ 



the element V (31, 97) 
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of an array is a number assigned to the site with 
coordinates X - 31, Y = 97. Altogether the 
Mray Fhas L X L = olomouts, and we have 
to reserve in the computer memory a space 
necessary to store j\T numbers. 

The work of the program begins with generat 
ing this array. Its elements are random numbers 
uniformly distributed between zero and unity. 

A random-number generator outputs a number 
1 / and this number is assigned to the array ele- 
ment V (1, 1). This means that this number is 
written into the corresponding memory cell ol 
the computer, and that from this moment on the 
computer “memorizes” that 7(1, 1) — y. The 
next number output by the generator is assigned 
to the element V (1, 2), and so forth. Thus all 
the elements of the array V are generated. 

Then a second two-dimensional array, that we 
denote by A, is formed. The elements of this 
array are zeros and unities, so that if, for exam- 
pie, K (25, 16) =0, this means that the site 
with coordinates X = 25, Y = 16 is b oc e , 
but if K (25, 16) = 1, this site is nonblocked 
The array K is generated by using the array V 
and a certain number t that falls within the in- 
terval from zero to unity. Varying t, we can 
change the number of blocked sites. 

The array K is generated by obeying the ioL 
lowing rule. Choose a site with coordinates A 
and Y. If V (X, Y)^t, then KJX, Y) = 1, and 
if V (X, Y) > t, then K (X, Y) = 0. The site 

with coordinates X and Y is treated as nonblocked 
in the former case and as blocked m the latter 
case. The quantities V being uniformly distri- 
buted within the interval from zero to unity, 
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w e can assume t being close to zero and obtain 
fliat almost all sites are blocked. Conversely, 
it t is close to unity, almost all sites are non- 
^locked. When t = 1/2, the number of blocked 
a *id nonblocked sites must be nearly 
e <lual. 

Making use of the distribution of random num- 
bers output by the generator, it is possible to 
r Hate the quantity t to the mean fraction x of 
n onblocked sites obtained as a result of the above- 
ascribed procedure. It can be shown (see Exer- 
cise 8) that t = x. However, this equality is val- 
H if we take a very large number JP of sites or 
if we generate many arrays K with the same t 
a hd then average the fractions of nonblocked 
s ltes found in each array. Actually, in each con- 
crete array, x may deviate from t to some extent, 
i*l both directions, but the equality holds the 
better, the larger JT is. 

The computer thus stores the array V , from 
vdiich we can generate an array K that describes 
vdiich site is blocked and which is not. The form 
pf the array K is dictated by the quantity t which 
i^ approximately equal to the fraction of non- 
blocked sites generated in this array. Smoothly 
Varying t , we can generate the distributions of 
blocked and nonblocked sites with smoothly 
Varying concentration of nonblocked sites, x . 

Exercise 

8. Prove that the mean fraction x of nonblocked 
s ites in an array K equals t. 
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Search for Percolation Path 

Let us assume that the array V is stored, the 
criterion t is fixed, and the array K with a certain 
fraction of nonblocked sites is found. 

Now the computer knows precisely which site is 
blocked and which is not, and we pass to the 
second stage of the program, namely, the search 
for percolation path. Assume that we look fo 
percolation from left to right. First we replace 
all unities in the leftmost column ( X = 1) by 
twos. The replacement consists in erasing a unity 
in the memory cell corresponding to a given ele- 
ment of the array K and in writing a two instead. 
A list of coordinates of the sites equipped with 
twos is compiled in the computer memory. Then 
the computer analyzes each site of this list. 
The computer calculates the nearest neighbors of 
the site under study and requests from the array 
K the information concerning these neighbors. 
If the nearest neighbor contains a unity, it is 
given a two, and its coordinates are added to 
the new list. After the first list has been analyz- 
ed, the computer memory contains a list of “sec- 
ond-generation” twos, that is, the list of unities 
replaced by twos because they were in contact 
with the first-generation twos. 

In order to economize on the computer memory, 
the first list is erased at this stage: it is no longer 
necessary; the appropriate memory cells are 
emptied. The computer starts an analysis of 
the second list and the formation of the third- 
generation twos. When this has been completed, 
the second list is erased and an analysis of the 
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third list is initiated. This analysis generates 
the fourth list, and so on. 

The number of twos in the array K grows in 
the course of this procedure. Twos are nonblocked 
sites linked by percolation to a nonblocked site 
of the leftmost column, that is, the twos label 
the percolation path. 

The search for percolation path terminates in 
two cases: 

1. A two appears on the right-hand side of a 
square. The computer recognizes the existence 
of percolation at the given value of t . 

2. There are no twos on the right-hand side of 
a square, and the analysis of the last list has not 
generated new sites marked with 2. This means 
that all paths broke down, and there is no per- 
colation at this t. 

Determination of the Threshold 

Assume that the computer recognized percola- 
tion at a given t. Then it diminishes t , and, mak- 
ing use of the same array F, finds a new array K 
with a reduced number of nonblocked sites. 
The search for percolation paths is then repeated. 
If percolation is found again, t is further reduced, 
until at a certain t percolation ceases. Then the 
interval between this value of t and the minimum 
value at which percolation was found is divided 
in two, and percolation paths are looked for at 
this intermediate value of t. If it is found that 
percolation paths are cut, the interval between 
this last value and the minimum value at which 
percolation was detected is again divided in two. 
If percolation is found, the interval to be divid- 
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Fie. 15. Pattern of the distribution of zeros, unities, and 
twos at the moment of percolation onset. The ^oken line 
shows the path by which twos percolated from the lett 
hand side to the right-hand side of the square In this 
case the computer did not stop calculations when the 
first. “2” aoneared on the right-hand side, but continued 



ed in two is the one between the last value of t 
and the value at which no percolation was found. 

Percolation threshold is thus bracketed 
with an arbitrary precision. If we find no percola- 
tion at the first chosen value of f, then t must be 
increased until percolation is detected, and this 
must be followed by “bracketing”. This method 
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enables us to find t corresponding to the perco- 
1 A a . tl ® n . threshold, with any prescribed accuracy 
At this t the computer finds the fraction of non- 
blocked sites, x, which is close to the value of t 
but is not necessarily equal to it. This value of 
x is taken for the percolation threshold found 
in th is . ex P e rmient. Figure 15 shows the first 

qh^qo 1011 Path from left t0 right found in a 
oU X dll square. 

.If en *i.9, a i experiments are then rerun many times 
with different sets of random numbers in the 
array V. This corresponds to changing the ran- 
dom sequence of blocked sites in the experiments 
with a wire mesh. The results of these experi- 
ments make it possible to determine the mean 
value x c {JT) of percolation threshold for a prescrib- 
ed number of sites, JT . (For this we have to add 
U fP u the obtained values of percolation thre- 
shold and divide the sum by the number of runs.) 

In order to find the true percolation threshold 
Xc ~ ** is necessar y to vary the numb- 

er of sites, ^ in the square and determine the func- 
tion x c (Jy)- This function must then be approxi- 
mated by an analytical expression of the type* 

z c (J r ) = x i .(oo)-f-iL 

jr K ' 

that is, we have to choose three quantities, 
f, c I'fA u ' and Y such that expression (5) fit in 
the best w ay the results obtained by means of 

There is no rigorous answer to the question whv this 

resuTte'of manvm ^ ^ P ower - la w manner, but the 

results ol many numerical experiments show that this 

is invariably the case in all percolation theory problems 



the computer. If this fit is achieved so that y > 0, 
we can say that it is the quantity x c (oo) that 
gives us the limiting value of x c . Indeed, by vir- 
tue of expression (5), 

lim * c (JT) = x c (oo) 

jr -+<x> 

The accuracy achieved by this procedure is the 
better, the greater the amount of data necessary 
to establish the dependence of x c ( JT ). In its turn, 
this amount of data is limited by the speed and 
memory of the computer employed. 

Exercise 

9. Look carefully at Fig. 15 and reconstruct the 
way by which individual groups of twos were 
formed. 
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Part II 

Various Problems 
of Percolation Theory 
and Their Applications 



Chapter 5 

Problems on Two-Dimensional Lattices 

We Are Planting an Orchard 
(the Bond Problem) 



Imagine that a vast orchard is being designed. 
The fruit trees in the orchard must be planted 
not arbitrarily but in a regular manner. They are 
to be located at the sites of a periodic lattice 
drawn on the ground. Many such lattices can be 
invented, but the following three will be suffic- 
ient for us here: a square, a triangular, and a 
hexagonal lattices (the latter is often referred 
to as the “honeycomb lattice”). The lattices are 
shown in Fig. 16. Land being expensive, it is 
natural to try and plant the trees as close to one 
another as possible, but this cannot be done for a 
number of reasons. One of the reasons is that 
the designers are afraid of infectious diseases of 
the trees. Let us assume that experts on tree 
infections supplied the following information*: 

* The author is not responsible for this information, and 
thus prays the reader not to be too serious about practical 
conclusions drawn from the solution of the problem formu- 
lated above. The problem is given merely to illustrate 
the potential inherent to percolation theory. 
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Fig. 16. Plane lattices: (a) square; (b) triangular; (c) ho- 
neycomb. 



| 



\ 



1. A tree blighted by the disease can infect 
only its nearest neighbors. 

2. Some nearest-neighbor trees inevitably in- 
fect one another when one of the trees has been 
infected. In other cases the disease is not trans- 
ferred (e.g. this may depend on the separation 
between the branches of a given pair of trees). 
A pair of trees in which the infection is inevit- 
able will be referred to as a connected pair. 

3. The experts also supply us with the form of 
the function x (a), mainly, the probability for 
a pair of nearest neighbors chosen at random to 
be connected. This probability depends on the 
distance a between the nearest neighbors in a 
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given lattice. Obviously, the function x (a) 
increases with decreasing argument a: two trees 
infect each other the easier, the smaller the se- 
paration between the trees. 

We are to answer the following question: What 
is the number of trees that a diseased tree can 
infect? This question can be given only a proba- 
bilistic answer. If a given tree forms a connected 
pair with one of its neighbors, this neighbor is 
infected. In their turn, the blighted trees infect 
their neighbors, and so forth. Therefore we can 
only ask: What is the probability for an infected 
tree to transfer the disease to a concrete number 
of trees in the orchard? 

At this juncture it is convenient to switch to 
the cluster language introduced in the preceding 
chapter. We assume that two neighbor sites 
with two trees forming a connected pair are 
linked by a bond represented by a piece of wire 
connecting the two sites. If two nearest-neighbor 
trees do not form a connected pair, the bond be- 
tween them (the wire) is broken (Fig. 17). 

Two sites will be regarded as connected if they 
are linked by an unbroken bond, or if they are 
linked by an unbroken chain of sites that are 
nearest neighbors and are linked by unbroken 
bonds (e.g. the sites A and 5, as well as the 
sites C and D , in Fig. 17 are connected). 

An ensemble of connected sites will be said 
to form a cluster. In the context of the given 
problem, the most important property of a cluster 
is that a diseased tree infects all the trees of 
its cluster and none outside this cluster. 

By definition, the fraction of unbroken bonds 
equals x. Further argument runs as in Chapter 3. 
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When x is small, unbroken bonds are mostly 
single, almost all clusters consist of two sites, 
three-site clusters are infrequent, and four-site 
clusters are even less so . When x is large, there 
exists an infinite cluster of connected sites. When 
x = 1 , this cluster comprises all the sites in the 
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Fig. 17. Fragment of a square lattice with broken bonds. 
Shown in the figure are three two-atom clusters (7, 2 , 3 ), 
one four-atom cluster (4), one six-atom cluster (5), and 
one ten-atom cluster (6), 

system. As x diminishes, this cluster loses some 
of the sites, and finally, at a certain critical 
value £ c , the infinite cluster ceases to exist. 

The infinite cluster is precisely the catastrophe 
from which the orchard must be saved. Let 
P h (x) be the probability for a randomly chosen 
site to belong to the infinite cluster. If x < x c , 
so that P b (x) =0, one blighted tree can infect 
only several other trees. And if x > x c , one 
blighted tree infects an infinite number of trees 
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in an infinite orchard with a probability P h (x). 
Therefore, if x < x c , the pocket of disease ac- 
cidentally brought to the orchard remains local- 
ized in the vicinity of the infected site, but it 
propagates all over the orchard if x > x c . 

In order to arrive at practical recommenda- 
tions, it is necessary to find x c and equate it to 
the function x (a) supplied by the experts. This 
yields the separation a c as the solution to the 
equation x ( a c ) = x c . The source of infection 
in the orchard remains localized if the separation 
between the trees exceeds a c , otherwise the in- 
fection propagates throughout the orchard. 

Percolation theory refers to the above-enun- 
ciated problem of determining x c as the bond 
problem in order to underline that here the ran- 
dom element lies in the bonds that can be un- 
broken or broken with a given probability. At 
first glance, the bond problem is quite similar 
to the site problem discussed in the preceding 
chapters. However, the two problems are not 
reducible to each other on a given lattice and 
have different answers. 

We shall have to somewhat complicate our 
notations in this and subsequent chapters. Let 
us denote the percolation threshold of the site 
problem by x s , and that of the bond problem by 
x h . These thresholds depend on the type of the 
lattice. Let us use the following abbreviations 
for the types of two-dimensional lattices: S for 
square, T for triangular, and H for hexagonal, 
or honeycomb. Then x s (H) denotes the percola- 
tion threshold of the site problem on a hexagonal 
lattice, x h (T) the percolation threshold of the 
bond problem on a triangular lattice, and so on. 
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The function P h (.z) introduced in this section 
for the bond problem should not be confused 
with the function P {x) defined earlier for the 
site problem. 

The bond problem can be formulated not only 
in cluster terms but also as a problem of percola- 
tion from one side of a square to the opposite 
side. Recall the experiment with a wire mesh 
with which we started the book. Some readers 
may have wanted to ask why it was necessary 
to cut simultaneously the four wires entering 
each site, instead of cutting randomly chosen 
individual wires (bonds). Now it can be clearly 
understood that by cutting the bonds the research- 
ers would find x ^ (S) instead of x s (S) which was 
actually determined in their experiment. Now 
we can explain why the site problem has been 
initially selected: it will be shown below that 
on a square lattice the bond problem has an 
exact analytical solution which yields that 
x b (S) = 0.5. Consequently, it would not be 
reasonable to carry out such a time- and labor-con- 
suming experiment for the sake of x b (S), with 
x s (S) being known only from approximate solu- 
tions. 

Exercise 

1. Find the function P h (x) for 1 — z < 1 for 
the three lattices shown in Fig. 16. 

Inequality Relating &b to x s 

When analyzing the site problem, the function 
P (x) is often replaced by a function P s {x) 
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related to P (x) by a relation 

P (x) = xP s (x) (1) 

By definition, P (x) is the probability for a ran- 
domly chosen site to belong to an infinite cluster. 
It can be written as the product of the probabili- 
ties of two independent events. In the language 
of the ferromagnetics problem, the first of these 
events is for a randomly selected site to be magnet- 
ic. The probability of this event equals x (see 
Exercise 1 to Chapter 1). The second event con- 
sists in this site being connected to an infinite 
cluster of magnetic sites. Therefore, the function 
P a (x) defined by formula (1) is the probability 
for a randomly selected magnetic site to be con- 
nected to an infinite cluster. In other words, 
P s (x) is the fraction of magnetic sites belonging 
to the infinite cluster, that is, the ratio of the 
number of sites belonging to the infinite cluster 
to the number of magnetic sites. We remind the 
reader that P (x) is the ratio of the number of 
sites belonging to an infinite cluster to the total 
number of sites. Naturally, the function P* (x) 
grows monotonically with increasing x, equals 
unity for x = 1, and vanishes for x^ x s . 

The English mathematician Hammersley, who 
was the first to enunciate percolation theory, 
proved a theorem that states that 

P s (x) P b (x) (2) 

Both functions P s and P b increase monotonically 
with increasing argument x. Therefore (Fig. 18), 
formula (2) implies that 

( 3 ) 



> 

i 



that is, the threshold for the bond problem is 
not greater than that for the site problem on any 
lattice (not necessarily two-dimensional). This 
result can be rewritten as a different inequality: 

1 x b 1 x s (4) 



that permits the following interpretation. Let 
us assume that we need to block electric current 




through a network of wires or block the flow of 
liquid through a network of pipes, and this can 
be done either by blocking network sites or by 
cutting network bonds (wires or pipes) connect- 
ing these sites. Inequality (4) means that the 
system is more readily blocked by cutting out 
its sites than by cutting the bonds. The fraction 
of blocked sites that terminates the flow is less 
than the fraction of broken bonds giving the 
same result. This conclusion appears quite natur- 
al because not one but all bonds entering a site 
are cut when this site is blocked. 
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Exercise 

2. Find the function P s ( x ) for 1 — 1 for 

the three lattices shown in Fig. 16. Compare the 
result with that obtained in Exercise 1 and check 
whether inequality (2) holds. 

Clarification. In Exercise 1 to Chapter 3 we 
recommended that you should find the function 
P ( x ) for the site problem, for 1 — a:<l. How- 
ever, there we meant only the principal term of 
the function, that is, P (x) = x. If this expression 
is substituted into formula (1), we obtain P s = 1. 
This is the correct result in the sense that 
lim P s (x) = 1. We suggest that the reader 

X-> 1 

find the small terms that make the function 
P s (x) differ from unity. Obviously, as x->i, 
these terms have zero for their limits. Conse- 
quently, the result may be written in the form 

ps (3) = 1 _ A (1 - x) n 

where A and n are positive numerical coefficients 
depending on the type of the lattice. 

Covering and Containing Lattices 

The site problem is more general than the bond 
problem. The bond problem is reducible to the 
site problem but on a different lattice said to 
cover the former. A covering lattice is constructed 
by using the following procedure: 

1. Place a site of the covering lattice in the 
middle of each bond of the initial lattice. 

2. Connect two sites of the covering lattice 
if and only if the bonds of the initial lattice on 
which these sites were placed meet at a site of 
the initial lattice. 
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a The result of such a construction is a new per- 
iodic lattice which is said to cover the initial 
lattice. 

Figure 19 shows the covering lattice in the 
case of an initial square lattice. Thin lines trace 
The initial square lattice. Semicircles show the 




Fig. 19. Covering lattice for the square lattice. 

places where the initial lattice has sites. The 
covering lattice consists of thin and solid lines, 
but there are no sites of the covering lattice 
where thin lines cross. The sites lie only at the 
crossings of solid lines and are marked by black 
circles. 

Each bond of the initial lattice joins three other 
bonds at one^end and three more bonds at the 
other end. Hence, each site of the covering lat- 
tice must be connected with six other sites. This 
is shown in Fig. 19. Each site is connected with 
four other sites by solid lines, and with two more 
sites by thin lines. 
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Now assume that a bond problem lias been 
formulated on the initial lattice, that is, a cer- 
tain fraction of randomly selected bonds has 
been cut. 

Assume now that if a bond of the initial lattice 
is broken, the site of the covering lattice placed 
at this bond is blocked. We thus obtain a site 
problem on the covering lattice. Its sites are 
randomly blocked, and the fraction of blocked 
sites equals the fraction of broken bonds in the 
initial lattice. 

Note that the existence of an infinite cluster 
of connected sites in the bond problem inevitably 
signifies the existence of an infinite cluster made 
of coupled unbroken bonds. Conversely, the ab- 
sence of an infinite cluster of sites signifies that 
bonds do not form an infinite cluster. 

As follows from the method of construction 
of the covering lattice, the existence of an infin- 
ite cluster formed by unbroken bonds in the 
initial lattice signifies that there exists an in- 
finite cluster formed by nonblocked sites in the 
covering lattice, and conversely, if there is no 
infinite cluster of bonds in the initial lattice, 
there is no infinite cluster of sites in the cover- 
ing lattice. Therefore, the percolation threshold 
of the bond problem in the initial lattice equals 
the percolation threshold of the site problem in 
the covering lattice. If the initial lattice is 
denoted by L, and the covering lattice by L cov , 
this statement can be written as formula 

x b (L) = x a (L eor ) (5) 

Introducing the concept of a lattice containing 
another lattice , it is possible to derive a number 
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of inequalities relating percolation thresholds on 
a number of lattices. Assume, for instance, that 
lattice L is obtained from lattice ^cont by cross- 
ing out a certain number of bonds. Then lattice 
Lcont is sa i<i t° contain lattice L. 

For instance, let us take a triangular lattice. 
If we erase all the bonds that are marked in 




Fig. 20. The triangular lattice contains the square lattice. 



Fig. 20 by two short lines, it will transform into 
the lattice shown on the right-hand side of the 
drawing. It can be easily seen that from the stand- 
point of the site or bond problems this new lattice 
is equivalent to a square lattice. Indeed, the fact 
that the angles between the bonds of the new 
lattice are not 90° is found to be immaterial if 
we analyze the relations between different sites 
(the new lattice can he just “straightened out”). 
The percolation threshold of the bond problem 
(as well as of the site problem!) on this lattice 
is exactly equal to the percolation threshold on 
the square lattice. For this reason, a triangular 
lattice is said to contain a square lattice. 

Now let us assume that a certain fraction of 
bonds in the containing lattice is broken. The 
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III k° nds of the containing lattice can be classified 
into the bonds that are common for the contain- 
ing and contained lattices, and the bonds that 
are specific for the containing lattice (these last 
are marked by two short lines in Fig. 20). The 
bonds being broken quite randomly, the frac- 
tion of broken bonds in one category of bonds is 
absolutely the same as in the other category, and 
equals the fraction of broken bonds in the whole 
lattice. Therefore, in order to obtain a contained 
lattice with the same fraction of broken bonds 
it is necessary to break additionally those bonds 
of the containing lattice that are intact but are 
specific for this lattice, that is, are completely 
absent in the contained lattice. 

This argument demonstrates that the number 
of unbroken bonds leaving each site of the con- 
taining lattice is not less than (is greater than 
or equal to) the number of unbroken bonds leav- 
ing the same site of the contained lattice. As 
a result, the probability for a randomly selected 
site to belong to an infinite cluster of the con- 
taining lattice is not less than that of the contain- 
ed lattice. This statement implies the inequality 

^(*)<^ cont (*) (6) 

The left-hand side of inequality (6) includes the 
function P b ( x ) calculated for the contained 
lattice, and the right-hand side includes this 
function for the containing lattice. By analogy 
to inequality (2) employing inequality (3), ine- 
quality (6) implies that 

I Zb (^cont)^^!) (A) 



(?) 
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Fig. 21. The square lattice contains the honeycomb lattice. 




Fi „ 22 Transformation of one unit cell of the lattice 
Swfon the rS”hand side of Fig. 21 into a honeycomb 
unit cell. 

that is, the containing lattice has a lower x b 
than the contained lattice. 

It has been already mentioned that a triangular 
lattice contains a square lattice. For this reason, 

* b (T) < z b (S) (8) 

Now assume that some bonds are deleted in 
a square lattice, as shown in Fig 21 This yields 
the lattice shown on the right-hand side of the 
drawing. Look attentively at this drawing. It 
is equivalent to a hexagonal lattice. Just pull 
it slightly upwards, somewhat warping the bonds 
and its cells (Fig. 22) will transform into the 
“honeycomb” lattice shown in Fig. 16. Conse- 
quently, a square lattice contains a hexagonal 
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lattice, and hence, 

x b (S) < x h (H) (9) 

Inequalities (8) and (9) imply that 

(T) < (H) (10) 

Let us turn now to the site problem. Assume 
that the same sites are blocked in the contain- 
ing and contained lattices (correspondingly, both 
have equal fractions of blocked sites). Let the 
contained lattice have an infinite cluster of intact 
sites. This means that it also exists in the con- 
taining lattice because additional bonds can 
only facilitate its formation. But if it is known 
that an infinite cluster exists in the containing 
lattice at a given fraction of sites, no conclusion 
can be drawn from this as to the existence of an 
infinite cluster in the contained lattice. The 
removal of some bonds between sites in the tran- 
sition from the containing to the contained lattice 
can be “lethal” for the infinite cluster. Conse- 
quently, the percolation threshold of the site 
problem in the contained lattice cannot be small- 
er than that in the containing lattice: 

x s (^cont)^^s (L) (11) 

Therefore, we can write a chain of inequalities 
for the triangular, square, and hexagonal lattices: 

(T) < z s (S) < * s (H) (12) 

quite identical to that for the bond problem. 

Let us turn again to Fig. 19 where the cover- 
ing lattice for the square lattice is shown. Imag- 
ine that the bonds traced by thin lines have been 
deleted. Obviously, this gives us a square lattice 
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merely rotated through 45°, which, of course, 
is quite immaterial for percolation theory prob- 
lems. 

The covering lattice for a square lattice thus 
contains a square lattice. By virtue of formula 

( 5 ), 

*s (£cov) = x b (S) ( 13 ) 

where L cov is interpreted as the covering lattice 
shown in Fig. 19. However, inequality (11) and 
the fact that this covering lattice contains a square 
lattice imply that 

*8 (£cov) < X S (S) (!4) 

Inequalities (13) and (14) imply that 

(S) < x s (S) (15) 

We have thus obtained that the bond problem 
threshold on a square lattice is lower than the 
site problem threshold. Inequality (15) is a parti- 
cular case of Hammersley’s general theorem that 
is embodied in formula (3) and was given earl- 
ier without proof. This theorem was not used in 
deriving (15), and we thus can say that the argu- 
ments given above prove this theorem in the case 
of a square lattice. 

“White” Percolation and 
“Black” Percolation 

Now let us look at the bond problem from a slight- 
ly different viewpoint. So far we normally have 
said that there are intact and broken bonds distri- 
buted randomly over a lattice and applied the 
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term ‘cluster” to an ensemble of sites connected 
by unbroken bonds. 

The problem can be reformulated in a more 
symmetric manner. Let us rename broken bonds 
“black”, and unbroken bonds “white”. An ensem- 
ble of sites connected by white bonds will be re- 
ferred to as a white cluster, and an ensemble of 
sites connected by black bonds as a black cluster 
(in our former terminology it was a “white” cluster 
that we called a’ cluster). Let the frac- 
tion'of white bonds be denoted, as before, by x. 
The fraction of black bonds will be denoted by 
q. Each bond has to be either black or white, 
so that q = 1 — x. 

After this reformulation we can speak of per- 
colation through both white bonds and black 
bonds. 

When the concentration x of white bonds is low, 
there is no infinite white cluster, but there is 
an infinite black cluster, that is, an infinite clust- 
er of sites connected through black bonds. Con- 
versely, when the concentration q of black bonds 
is low (i.e. at x nearly equal to unity), there is 
an infinite white cluster but not an infinite 
black cluster. 

When x varies from zero to unity, two events 
take place: the black infinite cluster vanishes, 
to be replaced by a white cluster, or, to say the 
same in different terms, percolation through black 
bonds vanishes, to be replaced by percolation 
through white bonds. But what is the sequence 
in which these events take place? 

The white and black bonds differ only in labels. 
It is obvious, therefore, that the critical concen- 
tration q i) at which percolation appears through 
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black bonds equals the concentration x h at 
which percolation appears through white bonds. 

Therefore, as x increases, percolation through 
white bonds appears when x = X\>, and percola- 
tion through black bonds stops when x = 
= 1 — q^ = 1 — Xft. The sequence in which 





I ( b) 

Fig. 23. 

these events take place depends on the sign of 
the difference x h — 0.5. 

■ If it is necessary for percolation through white 
bonds that the bonds be more than 50% white 
(this means that the number of white bonds must 
exceed that of black bonds at the percolation 
threshold), then at increasing x first percolation 
through black bonds ceases and then percolation 
through white bonds appears (Fig. 23a). In region 
/ in Fig. 23a there is percolation only through 
black bonds, in region III only through white 
bonds, and in region II there is no percolation 
through either white or black bonds. 
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If Xb < 0.5, first percolation through white 
bonds appears and only then percolation through 
black bonds ceases. In region / in Fig. 236 there 
is percolation through black bonds, in region III 
it proceeds through white bonds, and in region 
// through both black and white bonds. 

The same symmetric reformulation can be 
found for the site problem. We remind the reader 
that in this case all bonds are unbroken, but 
sites can be of two sorts. In the wire-mesh prob- 
lem we called them nonblocked and blocked. In 
the ferromagnetics problem these were magnetic 
and nonmagnetic sites. And now we introduce, as 
in the bond problem, a universal notation: 
nonblocked, or magnetic, sites are said to be 
white , while blocked, or nonmagnetic, sites are 
said to be black . White sites are said to be connect- 
ed if they are nearest neighbors or are connected 
by a string of white sites that are formed by pairs 
of nearest neighbors. Black sites can be connect- 
ed in the same symmetric manner. 

We are justified in speaking of percolation 
through white and black sites. If x s > 0.5, there 
is a region of x values with no percolation through 
either white or black sites (1 — x s < x < x s ). 
If x s < 0.5, the region x s < x < 1 — x s is 
the region in which percolation proceeds both 
through white and black sites. 

The symmetric approach proves constructive 
because sometimes it is possible to conclude from 
an analysis of the lattice that it cannot accommo- 
date either percolation through white or percola- 
tion through black, or that one of these must ne- 
cessarily be present. This information is suffic- 
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lent to make readily certain conclusions on per- 
colation threshold. 

For instance, let us consider the site problem 
on a triangular lattice. Assume that the lattice 
percolates through white sites. You will readily 
see that this excludes percolation through black 
sites. Let us assume that we study percolation 
of current from left to right on a very large net- 
work, as in the experiment whose description 
opened this book. But now the network is not 
a square but a triangular lattice. It is easy to 
understand that the structure of any triangular 
lattice is such that a rightward percolation path 
through white sites excludes any possibility for 
downward percolation through black sites. In- 
deed, black sites cannot “penetrate” through the 
broken line traversing the entire network from 
left to right and connecting white sites. 

As we explained in Chapter 3, an infinite cluster 
ensures percolation in any direction if the size 
of the system is sufficiently large. It then follows 
that a triangular lattice cannot contain, for the 
same value of x, two coexisting infinite clusters, 
one of white and another of black sites, that is, 
percolation through white sites and percolation 
through black sites are incompatible. Hence, 

(T)>0.5 

The same conclusion follows for square lattices: 
Zs (S) ^ 0.5 

A theorem has been proved for triangular lat- 
tices, stating that x s (T) = 0.5. The proof of 
this theorem cannot be given here, but the main 
points can be grasped more or less easily, making 

a 
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use of the concepts of percolation through white 
and percolation through black. Drawing various 
configurations of black and white sites, it^can 
be noticed that blocked rightward percolation 
through white sites necessarily implies downward 
percolation through black sites (square lattices 
do not possess this property!). 

A triangular lattice thus cannot simultaneously 
have percolation through white sites and percolation 
through black sites , but there must always be perco- 
lation through one sort of sites. This means that 
region II in Fig. 23 (# b in this figure must now 
be replaced by £ s ) degenerates to a point, that 
is, x s (T) = 0.5. 

In the cases of bond problems of this type an 
analysis can be facilitated by operating with dual 
lattices. 

Dual Lattices 

Only plane lattices can be dual. By definition, 
a plane lattice is a lattice that can be placed on 
a plane, with bonds intersecting only at points 
where lattice sites are located. For instance, all 
the lattices shown in Fig. 16 are plane, while 
the covering lattice in Fig. 19 is not plane be- 
cause its bonds intersect at points shown by arcs, 
but there are no sites at these crossing points. 
(These arcs act as “bridges” that uncouple right- 
ward and downward roads.) 

Each plane lattice divides the plane into cells. 
Lattice L d is said to be dual to lattice L if each 
bond of L d intersects one and only one bond be- 
longing to L, and vice versa, each bond of L 
intersects one and only one bond belonging to ZA 




Fig. 24. Construction demonstrating that square lattices 
are self-dual. The solid lines and black circles show the 
bonds and sites of the initial lattice; the dashed lines 
and open circles show the bonds and sites of the dual 
lattice. 




Fig. 25. Construction demonstrating that a triangular 
lattice is dual to a honeycomb lattice (and vice versa). 
Notations are the same as in Fig. 24. 
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‘ cannot at any point “break through it” and form 
a downward path. 

If a lattice is sufficiently large, the existence 
of a rightward path on the initial lattice signifies 
that the fraction of white bonds exceeds the thre- 
shold value: 



• z b ( L ) 



(16) 
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Furthermore, each cell of lattice L must con- 
tain only one site of lattice L d (and vice versa). 

You see from the definition that duality is a 
mutual property: if L d is dual to L, then L is dual 
to L d . Figures 24 and 25 show that a square lattice 
is dual to a square lattice, and a triangular and 
a hexagonal lattices are mutually dual. 

Let us return now to the bond problem. We 
assume that if a bond in the initial lattice is 
white (unbroken), the bond of the dual lattice 
intersecting the former is black (broken). There- 
fore, if the concentration of white bonds in the 
initial lattice is x , the concentration of white 
bonds in the dual lattice is q = 1 — x. 

In what follows it will be more convenient to 
interpret percolation threshold as that concen- 
tration of white bonds at which conductivity 
from left to right first appears (or disappears) 
in a very large wire mesh with electrodes soldered 
to it (see Fig. 1). This formulation of the problem 
is identical to that discussed at the very begin- 
ning of the book, but now the disrupted elements 
of the net are the bonds and not the sites of the 
lattice. 

Assume now that a rightward percolation path 
exists through white bonds of the initial lattice. 
It is easy to understand that this signifies no 
downward percolation through white bonds in 
the dual lattice. Indeed, by definition* a white 
bond of the initial lattice is intersected only 
by a black bond of the dual lattice. Hence, if 
the initial lattice contains a nowhere disconnected 
broken line composed of white bonds and tra- 
versing the entire network from left to right, this 
means that the white bonds of the dual lattice 



Conversely, the absence of a downward path 
through white bonds of the dual lattice indicates 
that the fraction of white bonds, q = 1 — #, 
in the dual lattice is less than the threshold 
fraction for the dual lattice, that is, 1 x < 
< z b (£ d )> or 

x > 1 - z b (L d ) (17) 

According to all this, all the values of x satis- 
fying inequality (16) also satisfy inequality (17). 
Hence, x h (L)^ 1 — x h (L d ), or 

Ib(i) +x b (Ld)>l ( 18 ) 

On a square lattice L = L d , so that inequal- 
ity (18) implies that 

z b (S)>0.5 ( 19 ) 

and making use of inequality (15), we addition- 
ally obtain 

z s (S)>0.5 ( 2 °) 

that is, we come to the same conclusion as in 
the preceding section. Recall that the experi- 
ment with a square lattice described in Chapter 
1 gave x s (S) = 0.59, which does not contradict 
inequality (20). 
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It has been rigorously proved for the square 
and triangular lattices that inequality sign in 
formula (18) must be replaced by equality sign , that 
is, 

(L) + * b (L d ) = 1 (21) 

This immediately leads to two new results: 
x b (S) - 0.5 (22) 

and 

z b (T) + * b (H) = 1 (23) 

We shall not give here the rigorous proof of 
formula (21): it would require that we introduce 
a number of new concepts that will not be nec- 
essary in the further exposition. However, this 
formula can be given a sufficiently clear interpre- 
tation. Tracing out various configurations of 
white and black bonds, it can be discerned that in 
the case of square and triangular lattices the 
absence of rightward percolation through white 
bonds in the initial lattice always signifies that 
there exists a downward percolation path through 
white bonds in the dual lattice. Let us take this 
statement for granted. Assume that the fraction 
of white bonds is such that there is no rightward 
percolation on the initial lattice. For a suffi- 
ciently large network this means that 

X < x b (L) (24) 

Correspondingly, the dual lattice contains a 
a downward percolation path. The fraction of 
white bonds in this lattice is 1 — x. Therefore, 

1 — x > x b (L d ), or 

x<l-x h m (25) 
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^he above arguments signify that any value of 
x satisfying inequality (24) must also satisfy 
inequality (25). This is possible if 1 — x b (L d )^ 
> # b (L), or 

z b (L) + z b (L d ) < 1 (26) 

Inequalities (18) and (26) yield equality (21). 

Formula (23) alone is not sufficient for find- 
ing separately x b (T) and x b (H). However, if 
we make use of the star-delta transformation 
known in the theory of electric circuits, we arrive 
at a relation between percolation thresholds of the 
bond problem on the triangular and hexagonal 
lattices. As a result, each of these thresholds is 
calculated: 

x b (T) = 2 sin (jt/18) « 0.347 296 

(27) 

x b (H) = 1 — x b (T) = 1 — 2 sin (ji/18) « 

» 0.652 704 

Exercise 

3. Look at Fig. 53 (see p. 252) illustrating 
the star-delta transformation. If you understand 
by looking at this drawing how to write the rela- 
tion mentioned above and arrive at formulas 
(27), you have done very well. If you do not go 
that far, do not despair because the problem is 
not that simple. The derivation of formulas (27) 
first obtained in 1963 by the English mathema- 
ticians Sykes and Essam was an important event 
in percolation theory. Now read attentively the 
text placed alongside the figure and you are likely 
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to feel pleased with the beautiful application of 
probability theory yielding such a nontrivial 
result. 

Results for Plane Lattices 

In conclusion we must give a summary table for 
percolation thresholds for plane lattices (Table 1). 



Table 1 

Percolation Threshold for Plane Lattices 



Lattice type 


*b 


*s 


Triangular 


0.3473 


0.5 


Square 


0.5 


0.59 


Hexagonal 


0.6527 


0.70 



Only two figures in this table, namely, x s (S) 
and x s (H), were obtained by approximate meth- 
ods. All the others represent exact solutions. In 
the next chapter we will show that the situation 
with three-dimensional lattices is much less 
satisfactory. In fact, not a single exact solution 
has been obtained. You need not regard this 
situation as strange. No universal methods exist 
for solving the problems of percolation theory 
analytically. Each exact solution mentioned 
above looks like a miracle. Rather, we should 
be surprised how many exact solutions have 
already been devised. 
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Exercise 

4. Let us revisit the orchard mentioned at the 
beginning of this chapter. Assume that the dis- 
tance separating the trees is chosen so as to 
satisfy the requirement that the fraction of con- 
nected pairs equals the threshold value. Let 
the function a (^) giving the distance between 
neighboring trees as a function of the fraction 
of connected pairs, x , be known. Naturally, the 
greater x is, the shorter a is because the trees 
infect one another the easier, the closer to one 
another they are planted. If x — se P a ~ 

ration between trees equals a (x^). Find the 
area per tree for the thus chosen separation, for 
three different lattices. The lattice with the 
least area per tree is the most profitable. Gould 
you say in advance which of the lattices gives 
the least area per tree, starting only with the 
information that the function a (x) monotoni- 
cally decreases with increasing x? 



Directed Percolation 

Assume now that the trees planted at the sites 
of a plane lattice form not an orchard but a for- 
est, and that this forest catches fire. Some near- 
est-neighbor trees have intertwined branches and 
readily pass on the fire. In accordance with our 
general terminology, we shall say that the sites 
containing such trees are connected by white 
bonds. Other nearest-neighbor trees do not light 
each other up. We shall say that the correspond- 
ing sites are connected by black bonds. The white 
and black bonds are distributed randomly over 
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the lattice, with the fraction of white bonds being 
equal to x. 

Our problem is to find the critical value x h 
such that at x < x h the pocket of fire remains 
localized, while at x > x h the fire spreads over 
the whole forest. 

Obviously, this is simply one more example of 
the bond problem. The value of x h can be found 
by looking it up in Table 1 (see p. 116). 

Assume now that the forest fire is accompanied 
by a strong wind, so that the flames propagate 
only along the wind. This results in a new and 
interesting problem, called the problem of di- 
rected percolation. 

We postulate that the forest was planted at the 
sites of a square lattice, and the wind is blowing 
along the diagonal of the square cells. In Fig. 26 
the wind direction is indicated by the arrow at 
the top of the figure, and the lattice has been 
rotated through 45° with respect to the standard 
arrangement of square lattices in this book. The 
problem is now reformulated as follows. Each 
white bond is converted into a vector whose 
arrow is placed so that the projection of the vec- 
tor onto the wind direction is positive. As before, 
black bonds do not let the fire pass whatever 
the direction, while white bonds let it pass only 
in the direction of the arrow. We want to find the 
critical fraction of white bonds beginning with 
which one tree in flames in an infinitely large 
forest can start a fire spreading to an infinitely 
large distance. 

The solid lines with arrows in Fig. 26 stand 
for white bonds; two percolation paths are shown: 
1 and 2. Path 1 traces the locus of the fire, while 



Wind 




& 



Fig. 26. Directed white bonds are shown me ^ 
lines with arrows; black bonds are shown by th n( j 

Fl»me propagates Irom Wit* stZnS rf pi* 2 

are marked by two dashes. 

path 2 does not correspond to the conditions of 
directed percolation: at two points marked by 
dashes the propagation goes counter to a white 
bond, that is, against the wind. 

We thus find that while m nondirected perco 
lation white bonds worked in both directions 
in directed percolation they work only in on 

direction. Hence, the critical frac V° n 

bonds in directed percolation, ^.cannot he 

smaller than in conventional percolation, th 

1S ' At^ the present time, approximate soiutions 
have been obtained for a number of problems 
Sh directed percolation. Thus, in the above- 
described bond problem on a square lattice 
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7T 0-63 or 0.64 (the results yielded by several 
methods slightly differ from one another). Recall 
that in the nondirected bond problem on a square 
lattice we had x h = 0.5. 

Quite a few physical problems are reducible 
o directed percolation. As an example, we can 
cite the motion of electrons in a strong electric 
held m a randomly inhomogeneous medium, 
these are media with properties that vary ran- 
domly from point to point, so that a moving 
electron runs into obstacles that it has to bypass 
At the same time, the uniform electric field acts as 
a wind that drives electrons in one direction only 
Directed percolation also realizes in the prob- 

rt e A ec A ric conduct ance of a wire net 
(see Chapter 1) if we assume that each bond be- 
ween the sites of the net contains a diode (a 
rectifying element) that lets the current flow 
in one direction only. Then the fraction of broken 
bonds at which the current through the net drops 
to zero corresponds to the threshold of directed 
percolation. A combined problem, with diodes 
present not in each bond of the net, has also been 



Chapter 6 

Three-Dimensional Lattices 
and Approximate Evaluation 
of Percolation Thresholds 



As was shown in Chapter 5, percolation threshold 
essentially depends on lattice type. The purpose 
of the present chapter is to explain qualita- 
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tively what exactly are the properties of lattices 
that are important for percolation thresholds. 
Having understood this aspect, it is possible to 
predict percolation thresholds (to within 10%) 
without solving the problem. The ability to come 
up with such predictions (be the prediction 
accurate or not) appears to be valuable because, 
first, the number of different lattices is quite 
substantial, and second, the computation of a 
threshold takes (on top of a high skill) about half 
an hour of computation time of the most advanced 
computers. 

Moreover, percolation problems are not re- 
stricted to lattice problems only. We shall be 
able to see that applications mostly require 
solving problems that are formulated not on 
lattices. It is found that the ideas of approximate 
evaluation presented in this chapter are quite 
fruitful for nonlattice problems. It was on the 
basis of these ideas that a number of percola- 
tion thresholds for nonlattice problems had been 
predicted with high accuracy long before these 
problems were solved with computers. 

In order to collect the necessary experience, 
we have to go beyond the framework of plane 
, lattices discussed in the preceding chapter and 
turn to three-dimensional lattices. 

Three-Dimensional Lattices 

The simplest of the three-dimensional (or simply 
3D) lattices is the simple cubic lattice . It is shown 
in Fig. 12. Its unit cell is a cube shown in Fig. 27. 
The vectors a x , a 2 , a 3 are called translation vec- 
tors. Elongating each of these translation vectors 
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an integral number of times (n lt n 2 , n 3 ) and then 
adding up the resultant vectors, we can obtain 
a vector R nt> „ 2> beginning at the origin of 
coordinates and ending at any site of the simple 
cubic lattice: 

n 2 , n, = Wi^i + n 2 a 2 + w 3 a 3 

In cubic lattices the three vectors a lt a 2 , a 3 are 
of identical length a, so that the numbers n 1 , 




»i 



Fig. 27. Unit cell of a simple cubic lattice. 

n 2 , n 3 are simply three Cartesian coordinates of 
lattice sites in units of a. The simple cubic lat- 
tice is said to be generated by translation (parallel 
displacement) of the cubic unit cell by vectors 
that are integral multiples of a lt a 2 , a 3 . 

The most important characteristic of a lattice 
is the number of nearest neighbors (also called 
the coordination number ), denoted by z. In the 
simple cubic lattice z — G. 

Such alkali halides as NaCl (common salt), 
KC1 (rock salt), LiF, Nal, and a number of 
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others crystallize into the simple cubic lattice. 
The alkali metal ions (e.g. Na + ) alternate in 
this lattice with the halogen ions (e.g. Cl“). 

Body-centered cubic lattice ( abbreviated to bcc). 
This lattice can be composed of two simple 




(«) *' (.b) 



Fig. 28. (a) Body-centered cubic lattice; (b) unit cell of 
the bcc lattice. 

cubic lattices by displacing them relative to 
each other by one half of the body diagonal of 
the unit cell. (By definition, the body diagonal 
of a cube is the line connecting the opposite 
vertices of the cube and passing through the cen- 
tral point.) 

A bcc lattice is shown in Fig. 28a. Its unit 
cell is shown in Fig. 28 b. The entire bcc lattice 
can be generated by translating this unit cell 
by vectors a l5 a 2 , a 3 ; however, it will not be suffici- 
ent to place the vector R ni> „„ at the origin if 
we want to generate all sites of the lattice. We 
also need to place the origin of the vector at 
the central atom. Only then will the two simple 
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cubic sublattices of the bcc lattice be translated 
together. 

The nearest neighbors of each site of a bcc lat- 
tice lie in the direction of the body diagonal of 
the unit cell (see Fig. 28 b). Each site has eight 
such neighbors, so that z = 8. The distance to 
a nearest neighbor equals one half of the body 
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Fig. 29. Unit cell of a face-centered cubic lattice. 

diagonal, that is, |Y3a/2, where the letter a 
denotes, as before, the length of the edge of the 
unit cell. 

Univalent alkali metals, such as lithium, sod- 
ium, potassium, rubidium, cesium, bivalent 
barium, and a number of other substances 
crystallize into bcc lattices. 

Face-centered cubic lattice (abbreviated to fee). 
A unit cell of an fee lattice is shown in Fig. 29. 
In comparison with a simple cubic lattice, it 
contains additional sites placed at the centers 
of each face. In order to translate this cube, the 
origin of the vector R ni> n . it V9 must be placed 
at the origin of coordinates and also at the centers 
of the three nonopposing faces of the cubic unit 



cell. The nearest neighbors of each site of an 
fee lattice lie in the direction of face diagonals. 
The distance to a nearest neighbor is a/]/~ 2. Each 
of the three mutually perpendicular plane in- 
tersecting at a given site contains four nearest 




neighbors of this site, so that the number of 
nearest neighbors, z, of a site is 12. 

Substances that crystallize into fee lattices 
are such metals as copper, silver, gold, alumin- 
ium, lead. 

The last lattice that we introduce into our ana- 
lysis is the diamond lattice. It is shown in Fig. 30a. 
You obtain it if you imagine two face-centered 
cubic lattices displaced with respect to each other 
along the body diagonal of the cubic unit cell 
by one fourth of its length. 

Substances that crystallize into diamond-type 
lattices are elements of the fourth group of the 
periodic table: carbon (diamond), as well as the 
two most important semiconductor elements: 
germanium and silicon. All these elements are 
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tetravalent, and the atoms in their lattices are 
bonded by covalent forces. To simplify, we can 
imagine that each atom has four “aims” corre- 
sponding to four valence electrons. An atom in 
a lattice then holds hands with its four nearest 
neighbors. The diamond lattice is perfectly 
suited to this type of bonding. Each site of this 
lattice lies at the center of a regular tetrahedron 
formed by other sites (Fig. 30b). The number of 
nearest neighbors is z = 4. 

Percolation Thresholds for 3D Lattices 

The bond and site problems are formulated for 
3D lattices exactly as for plane lattices. As there, 
we again assume that bonds connect only near- 
est-neighbor sites. 

Table 2 summarizes percolation thresholds 
of the site and bond problems in the 3D lattices 
described above. As we have already mentioned, 
no exact solution was obtained in the 3D case. 
All the results listed in Table 2 were obtained 
by various approximate methods that, as a rule, 



Table 2 

Percolation Threshold for 3D Lattices 



Lattice type 


*b 


X s 


Simple cubic 


0.25 


0.31 


Body-centered cubic 


0.18 


0.25 


Face-centered cubic 


0.12 


0.20 


Diamond 


0.39 


0.43 



employ computers. Naturally, certain slight 
discrepancies are found between the results 
published in the scientific literature. Table 2 
is a selection of results that are the most reliable, 
from our point of view. 

Now the problem is to try and understand, by 
looking at this table and Table 1 (see p. 116), 
which summarizes the results for plane lattices, 
why some lattices have relatively high percola- 
tion thresholds, and why other lattices have 
relatively low thresholds. We begin with the 
bond problem. 

Factors Determining Percolation Threshold 
in the Bond Problem 

If all bonds are intact, each site is connected with 
z other sites, where the number z of nearest neigh- 
bors varies considerably among lattice types. 
At a fixed fraction x of intact bonds, each site 
is connected, on the average, with zx other sites. 
Let us try and test the following hypothesis: 
Is it realistic for the quantity zx that gives the 
average number of sites with which each site 
is connected to carry sufficient information for 
deciding whether percolation is present or absent 
in the lattice? Can it be that we need no other 
information on the properties of a lattice but its 
coordination number z, and that percolation 
sets in in all lattices at the same value of zx? 
It is rather clear that this hypothesis cannot be 
very accurate. But is it approximately valid? 

This is very easy to prove or disprove. Let 
us find the product zx b for all lattices with known 
percolation thresholds of the bond problem. If 
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this product is universal, that is, identical or 
at least nearly identical for all lattices, the pro- 
posed hypothesis holds or holds approximately. 



Table 3 

The Product zx b for Different Lattices 



Lattice type 


z 


*b 


ZX b 


Plane lattices 








Square 


4 


0.5 


2.0 


Triangular 


6 


0.35 


2.1 


Honeycomb 


3 


0.65 


2.0 


3D lattices 








Simple cubic 


6 


0.25 


1.5 


Body-centered cubic 


8 


0.18 


1.4 


Face-centered cubic 


12 


0.12 


1.4 


Diamond 


4 


0.39 


1.6 



The required data are summarized in Table 3. 
You see that in plane lattices we have 

zx h = 2 (1) 

with an error less than 10%, and in 3D lattices 
zx b = 1.5 (2) 

The hypothesis of the universal mean number of 
bonds per site required for the onset of percola- 
tion is thus not accurate, but it does hold approx- 
imately. If we remember that each of the quan- 
tities z and x b varies by at least a factor of two 
both in the group of plane lattices and in that 
of 3D lattices, the accuracy with which the pro- 
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duct zx b is constant within each group must be 
considered remarkably high. 

We conclude that percolation threshold of the 
bond problem can be approximately evaluated 
if we know the number of nearest neighbors, and 
use formula (1) for plane lattices and formula 
(2) for 3D lattices. Percolation threshold of the 
bond problem is most sensitive to the number of 
nearest neighbors, but is much less sensitive to 
all other properties of lattices (e.g. to the number 
of next-to-nearest neighbors, i.e. the sites 
that are second in their closeness to a given site 
after its nearest neighbors). 

We have thus arrived at a very simple and at 
the same time a relatively accurate method of 
evaluating percolation thresholds of the bond 
problem applicable to arbitrary lattices. 

How to Evaluate Percolation Threshold 
in the Site Problem 

Let us analyze now a similar scheme for the site 
problem. It will be natural to test first the fami- 
liar version, that is, check the variation of the 
quantity zx s among different lattices. It will be 
readily found that the product varies almost as 
each of the quantities z and x s does independent- 
ly. One should not be surprised: the product 
zx b has a well-defined physical meaning for the 
bond problem, namely, that of the mean number 
of intact bonds per site. In the case of the site 
problem, a bond is effective if it connects two 
white sites, but is ineffective in all other cases. 
Consequently, the product zx s is virtually 
meaningless. 
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Fig. 31. Construction of tangent circles in a honeycomb 
lattice. The lattice itself is shown in Fig. 16c. The circles 
have radii equal to one half of the distance to the nearest 
neighbor. The white sites correspond to white circles, and 
the black sites to black circles. The percolation paths 
through white circles are traced with the solid lines. 

In 1970 the American physicists H. Scher and 
R. Zallen proposed a different method for evalu- 
ating percolation threshold in the site problem. 
Their idea was to put each site in correspondence 
with a specific part of space. Then they went on 
to say that percolation through white sites appears 
when the fraction of space occupied by these 
sites exceeds a certain critical value slightly 
dependent on the lattice type. 

Imagine around each lattice site a sphere (or a 
circle if we deal with a plane lattice) of radius 
equal to half the distance to the nearest neigh- 
bor. The spheres (circles) constructed around 
neighboring sites are then tangent to one another 
(Fig. 31). A white site is assigned a white sphere, 
and a black site a black sphere. If two white sites 
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are connected, there is a path between them 
through tangent white spheres (see Fig. 31). There- 
fore, the onset of percolation means the formation 
of infinitely long paths through tangent white 
spheres. 

Now we assume that percolation sets in when 
the fraction of the total volume (or of the total 
area) occupied by white spheres (or circles in a 
plane lattice) exceeds a certain critical value in- 
dependent of lattice type. In order to test this 
hypothesis, we need to calculate the fraction of 
volume occupied by white spheres for x — x s 
for different lattices with known values of x s and 
then compare the results. 

First we need to calculate the fraction of volume 
filled up by white spheres for x = 1, that is, 
in the case when all spheres are white. This 
quantity is denoted by letter / and referred to as 
filling jactor. The filling factor equals the fraction 
of volume occupied by the spheres constructed 
around each lattice site and having a radius equal 
to half the distance to the nearest neighbor. Fil- 
ling factor essentially depends on lattice type 
and must be calculated for each specific lattice. 

The fraction of volume occupied by white spheres 
for x < 1 is found by multiplying the filling 
factor by the fraction of white spheres, that is, 
by x. The fraction of volume filled up by white 
spheres is thus fx. At percolation threshold it 
equals fx s . If the hypothesis on the universality 
of the fraction of volume at which percolation 
sets in is correct, the product fx s must be identi- 
cal for all lattices. 

The filling factors for a number of lattices are 
listed in the second column of Table 4. In order 
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Table 4 

The Product fx s for Different Lattices 



Lattice type 


/ 


*s 


fx s 


Plane lattices 








Square 


0.79 


0.59 


0.47 


Triangular 


0.91 


0.5 


0.46 


Honeycomb 


0.61 


0.7 


0.43 


3D lattices 








Simple cubic 


0.52 


0.31 


0.16 


Body-centered 


0.68 


0.25 


0.17 


Face-centered 


0.74 


0.20 


0.15 


Diamond 


0.34 


0.43 


0.15 



to illustrate how they were calculated, let us de- 
termine / for the honeycomb lattice given in 
Fig. 31. It was shown in Exercise 4 to Chapter 5 
that the area per site of the honeycomb lattice is 
(3 / 3/4) a 2 , where a is the length of the side of 
the unit cell. This result has the following mean- 
ing: let us draw on a plane, on which the lattice 
was constructed, a square, a rectangle, a circle, 
or any other geometric figure, but necessarily 
such that its size be much greater than the distance 
between the neighboring sites of the lattice. 
Divide its area by the number of lattice sites that 
fell within this figure. The area per site is the 
limit to which this ratio tends when the figures 
become infinitely large. 

The fraction of area occupied by circles equals 
the limit of the ratio of the area occupied by the 



circles to the area of the large figure. The area 
occupied by the circles equals the product of the 
number of sites belonging to the large figure by 
the area of a single circle. In other words, the 
quantity / equals the ratio of the area of one cir- 
cle to the area per site. 

The circles drawn in Fig. 31 have the radius 
a/2, and hence, their area is na 2 /4. This gives 

no*/ 3/3a » = w Q.gQg 
7 4 / 4 3/3 

The filling factors for other lattices are calculated 
in a similar manner; you can see from Table 4 
that / varies within a wide range. 

The last column of Table 4 lists the products 
fx s . You can see that the hypothesis that fx s is 
independent of lattice type does not hold too 
well. However, its variations are small both 
within the group of plane lattices and within the 
group of 3D lattices. Hence, to within 10 to 15% 
the following formulas hold: 

fx s = 0.5 (^) 

for plane lattices, and 

fx s = 0.16 (4) 



for 3D lattices. 

The calculation of the filling factor / being rel- 
atively simple, formulas (3) and (4) make it 
possible to evaluate percolation threshold of the 
site problem for arbitrary lattices. 

It is readily understood that the critical fraction 
of volume occupied by white spheres, at which 
percolation appears, decreases monotonically 
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with increasing dimensionality of space. In one- 
dimensional space, that is, in a linear string of 
sites, percolation through white sites vanishes 
at no matter how small concentration of black 
sites. Even a single black site blocks the perco- 
lation path because no bypass is possible. The 
possibility to go around black sites appears in 
plane (two-dimensional) lattices, and such pos- 
sibilities become more numerous in 3D lattices 
because detours are not restricted to a single 
plane. 

The critical volume concept proves fruitful not 
only in lattice problems. In Chapter 9 we shall 
encounter a problem in which white and black 
balls are not located at lattice sites but are ran- 
domly poured into a jar. We shall be interested 
in percolation through white balls in contact with 
one another. This percolation is also found to arise 
when the volume occupied by white balls comes 
to about 0.16 of the total volume. This result 
changes only slightly if the balls have unequal 
radii. 

Chapter 10 treats the problem of space painted 
in a random manner by white and black paints. 
It is found that percolation through regions of one 
of these colors appears in the 2D case when the 
fraction of surface area painted white (or black) 
is exactly equal to 0.5, and in the 3D case when 
the fraction of volume painted white (or black) 
is approximately equal to 0.16. 

Exercise 

1. Check whether the filling factors / listed 
in Table 4 were correctly calculated. 
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Chapter 7 

Ferromagnetics 

with Long-Range Interaction. 

The Sphere Problem 

So far, when considering site and bond problems, 
we have always assumed that each site can be di- 
rectly connected only with its nearest neighbors, 
and that the connections between farther removed 
sites are established through strings of sites each 
of which is connected with its nearest neighbors. 
In this chapter the site problem is generalized to 
the case in which the sites that are not nearest 
neighbors are directly connected with each other. 
This problem may be of practical importance, and 
so it may prove useful to know that it has been 
sufficiently well elaborated. 

If the number of sites with which a given site 
is connected becomes very large, the site problem 
transforms to a quite new problem referred to as 
the sphere problem. In fact, this problem plays an 
especially important role in percolation theory: 
it is used to understand the transition to metallic 
conductivity taking place in semiconductors as 
the impurity concentration in them increases, its 
solution is at the basis of the theory of hopping 
conduction in semiconductors, which is an impor- 
tant and interesting phenomenon observed at very 
low temperatures. Consequently, numerous scien- 
tists worked with the sphere problem, and intere- 
sting results have been obtained for this and 
similar problems. 

The sphere problem is also interesting because 
this is the first nonlattice problem that we en- 
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counter in this book. The random elements with 
which it operates are defined not at sites of a 
periodic lattice. 

Ferromagnetics with Long-Range Interaction 

Let us return to the problem of a ferromagnetic 
with nonmagnetic impurity atoms, discussed in 
Chapter 3. There it was assumed that magnetic 
atoms orient their magnetic moments parallel to 
each other only if they are nearest neighbors or 
if they are connected through a string of magnet- 
ic atoms that form successive pairs of nearest 
neighbors. But if all nearest neighbors of a mag- 
netic atom were nonmagnetic atoms, the magnet- 
ic moment of this atom was assumed to be arbi- 
trarily oriented. 

This model was based on the fact that the inter- 
action between magnetic moments resulting in 
the parallel arrangement decreases very rapidly 
with distance; in fact, the decrease is so steep 
that the magnetic moments that are not nearest 
neighbors simply “ignore” one another, that is, do 
not interact at all. 

In crystallography, the group of all nearest- 
neighbor atoms of a given atom is usually called 
its first coordination group , and it has already 
been mentioned that the number z of nearest 
neighbors is called the coordination number. 
The ensemble of equivalent atoms that are 
next-to-nearest neighbors of the given atom 
is called the second coordination group , and 
so on. As an example, let us consider the 
simple cubic lattice (see Fig. 12). The first coor- 
dination group in this lattice consists of six atoms 
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located at the edges of the cube projecting from 
the initial atom. The second coordination group is 
formed by 12 atoms located at the diagonals of the 
cube’s faces passing through the initial atom. 
And finally, the third coordination group consists 
of eight atoms located at the body diagonals of 
the cube passing through the initial atom. 

In Chapter 3 we assumed that the interaction 
between atomic magnetic moments vanishes be- 
yond the first coordination group. Under this 
assumption, the calculation of the critical con- 
centration of magnetic atoms at which spontane- 
ous magnetization appeared (or vanished) re- 
duced to the site problem with bonds only between 
nearest neighbors. 

The assumption of short-range interaction does 
not always hold, so that it is justifiable to discuss 
a problem in which the interaction between mag- 
netic moments covers several coordination groups, 
and to determine the critical fraction of magnetic 
atoms resulting in spontaneous magnetization. 

This problem reduces to the site problem in 
which bonds are formed not only between nea- 
rest neighbors. In fact, its formulation contains 
nothing new. Sites may be white and black (mag- 
netic and nonmagnetic). Two white sites are re- 
garded as connected if the bonds connect the 
coordination groups in which they are located. 
If site A is connected with site B , and site B with 
site C, then A is connected with C. The ensemble 
of connected sites forms a cluster. The percola- 
tion threshold is defined as the fraction of white 
sites at which an infinite cluster is born. 

When bonds are made to extend to farther coor- 
dination groups, the percolation threshold x s must 
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Table 5 



Percolation Threshold in the Site Problem with 
Bonded Nonnearest Neighbors* 



Lattice type 


z 


*s 


Zx s 


Plane lattices 






Honeycomb, 1 


3 


0.700 


2.10 


Square, 1 


4 


0.590 


2.36 


Triangular, 1 


6 


0.500 


3.00 


Square, 1, 2 


8 


0.410 


3.28 


Triangular, 1, 2 


12 


0.295 


3.54 


Honeycomb, 1, 2, 3 


12 


0.300 


3.60 


Square, 1, 2, 3 


12 


0.292 


3.50 


Triangular, 1, 2, 3 


18 


0.225 


4.05 



3D lattices 



Diamond 


4 


0.425 


1.70 


Simple cubic, 1 


6 


0.307 


1.84 


Body-centered, 1 


8 


0.243 


1.94 


Face-centered, 1 


12 


0.195 


2.34 


Body-centered, 1, 2 


14 


0.175 


2.45 


Simple cubic, 1, 2 


18 


0.137 


2.47 


Face-centered, 1, 2 


18 


0.136 


2.45 


Simple cubic, 1, 2, 3 


26 


0.097 


2.52 


Body-centered, 1, 2, 3 


26 


0.095 


2.47 


Face-centered, 1, 2, 3 


42 


0.061 


2.56 



* The listed thresholds for a single coordination group are 
not always exactly equal to the data given in the preceding 
tables. The reason for this is that those tables gave the data 
that we considered the most reliable, while the present table 
compares the data calculated by the same method. 
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obviously diminish. The greater the number of 
bonds originating at a given white site, the great- 
er the probability that at least one of these 
bonds will lead to another white site, and cor- 
respondingly, the smaller the number of white 
sites necessary to ensure percolation. 

Clearly, this problem is not a bit simpler (rath- 
er, it will be more complex) than the standard 
site problem. However, a number of problems of 
this type was solved by various approximate 
techniques, and the results of one of them are 
listed in Table 5. The first column of the table 
gives the lattice type and the number of coordi- 
nation groups to which the bonds were extended. 
The second column shows the number of sites, 
Z, with which each site is connected, that is, the 
total number of sites located in the considered 
coordination groups (in the case of a single coor- 
dination group it coincides with the coordination 
number z). 

The last column of the table gives the product 
Zx s . This product, as we indicated in the preced- 
ing section, is strongly dependent on lattice 
type in the case of the site problem with bonds 
only in the first coordination group. However, as 
we find in Table 5, at high values of Z it changes 
the slower, the greater Z is. This is especially 
clear for three-dimensional lattices where large 
values of Z are used. Obviously, the product Zx s 
tends to a number of the order of 2. 6-2.7 inde- 
pendent of lattice type. 

In the 3D case the number 2.7 is now consid- 
ered the most likely value (to within ±0.1) of 
the quantity B c defined as the limit of Zx B for 
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very large Z: 



B c =limZx 8 (1) 

Z-+00 

It will be explained in the subsequent sections 
why this limit exists and why it is independent 
of lattice type, being only a function of the di- 
mensionality of space, that is, of whether the lat- 
tice under consideration is plane or three-dimen- 
sional. 

In order to understand this feature, it is neces- 
sary to study the so-called sphere problem. 

Exercise 

1. Find the location of 42 sites belonging to the 
first three coordination groups of an fee lattice. 

The Sphere (Circle) Problem 

Now we shall consider a different problem or, 
rather, a problem that seems to be different at 
first glance. Assume that a plane is covered with 
circles of identical radius R whose centers are 
spread on the plane quite randomly and, on the 
average, uniformly. This means that both coor- 
dinates of centers of circles are random numbers 
distributed uniformly in the range from zero to 
L, where L is a very large length (in comparison 
with R ) characterizing the size of the system under 
discussion. The important distinctive feature of 
this problem is that there is no limitation on 
the extent to which the circles overlap. The 
mean number of centers of circles per unit area is 
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I N. In other words, N is the concentration of 
J centers of circles. 

Two circles are said to be mutually connected 
if the center of one of them lies within the other; 
such circles are sometimes said to be inclusive. 
If circle A is connected with circle 5, and B is 




Fig. 32. Percolation paths through inclusive circles are 
traced by the broken lines. The dots mark the centers of 
circles. 



connected with (7, then A is connected with C. 
Therefore, too far removed circles can be connect- 
ed via a string of inclusive circles (Fig. 32). 

The problem now is to find the critical concen- 
tration N at which percolation sets in through 
inclusive circles, that is, paths are formed tra- 
versing the whole system and composed of inclu- 
sive circles. (In other words, an infinite cluster is 
formed out of mutually connected circles.) 

The problem involves two parameters: the con- 
centration N and radius R. (We have left out 
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the size of the system, L, but it is clear that if 
the system is sufficiently large, the critical value 
N is almost independent of L .) It is not difficult, 
however, to ascertain that the presence or absence 
of percolation depends not on two parameters but 
only on one, namely, on the dimensionless pro- 
duct NR 2 . (The dimensionality of concentration 
in a plane problem is cm -2 .) Let us choose for 
this parameter the mean number of centers of 
circles within a circle. It equals 

B = nNR 2 

The following arguments will immediately 
show that percolation sets in at a certain value 
of the parameter Z?, regardless of the specific val- 
ues of N and R making up B. Imagine looking at 
a plane with the circles drawn on it. Let us mag- 
nify this picture severalfold, for instance, through 
a projector. This will be a transformation vary- 
ing N and R but leaving B invariant because the 
mean number of centers of circles within a circle 
will not be changed by this magnification. 

It is also easy to understand that this transfor- 
mation does not affect percolation. If there was 
no percolation in the initial picture, there will 
not be any in the magnified picture, and con- 
versely, if there was percolation through inclusive 
circles in the initial picture, this percolation will 
not vanish because of magnification. 

The transformation changing N and R but not 
B thus does not affect percolation. Consequently, 
it depends only on the value of the parameter B 
whether there will be or not be percolation in a 
system. Percolation is present when B is large, 
and absent when it is small. 
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This new problem that has just been formulat 
ed above is referred to as the circle problem. Its 
three-dimensional analogue is called the sphere 
problem. The sphere problem is formulated as 
follows. The coordinates of centers of spheres of 
radius R are generated in the three-dimensional 
space by a suitable random-number generator. 
Two spheres are said to be connected (or inclu- 
sive) if the center of one sphere lies within another 
sphere. We need to determine the critical con- 
centration of centers at which percolation sets in 
through inclusive spheres. 

It will be readily understood that, as in the 
plane problem, the presence or absence of perco- 
lation is determined exclusively by the value 
of the parameter B which is the mean number of 
centers of spheres within a sphere: 

B = (4/3) nNR 3 

where N is the mean number of centers of spheres 
per unit volume (the dimensionality of volume 
concentration is cm" 3 ). 

We have already mentioned that the sphere 
problem has important applications to the theory 
of electric conduction in semiconductors at low 
temperatures. For this reason, it was analyzed 
by a number of authors resorting to most various 
techniques. By today’s data, the critical value B c 
at which percolation through spheres sets in is 
2.7 -I- 0.1. The circle problem was not studied as 
thoroughly, so that a substantial spread is found 
among the results published by different authors. 
It appears that B c = 4.1 ± 0.4. 
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The Circle (Sphere) Problem 

Is the Limiting Case of the Site Problem 

Let us return to the site problem in which not only 
nearest neighbors are connected, and explain 
why the limit on the right-hand side of formu- 
la (1) exists, why it is independent of lattice 
type, and finally, why we denote it by the same 
symbol B c as the threshold values of the circle 
and sphere problems. 

We begin with slightly reformulating the site 
problem. To be specific, we first discuss plane 
lattices. Generalization to 3D lattices proves to 
be quite simple. 

Let us surround each white site with a circle of 
radius R that was chosen to be greater than the 
distance from this site to the sites of the farthest 
coordination group with which the given site is 
connected, but smaller than the distance to the 
sites of the next coordination group. We assume 
that two circles are connected if the white sites 
around which they were drawn are connected. 
This means that two circles are connected if the 
center of one of them lies within the other, that is, 
if the circles are inclusive. 

The formation of an infinite cluster of connect- 
ed white sites is equivalent to the onset of per- 
colation through inclusive circles. In other words, 
percolation through the circles of radius R drawn 
around white sites sets in at a critical fraction 
of white sites, x s (Fig. 33). 

The meaning of the quantity Zx is easily un- 
derstood. The quantity Z is the number of sites 
(both black and white) that lie within a circle. 
The product Zx is the mean number of centers of 
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other circles that lie within a circle (or the mean 
number of white sites within a circle). The quan- 
tity Zx s is the mean number of centers of circles 
within a circle at which percolation sets in, 
that is, at which infinite paths can be traced 
through inclusive circles. 

It is now obvious that the product Zx s has the 
same meaning as the quantity R c in the circle 




Fig. 33. Percolation path through inclusive circles on a 
square lattice. The interaction is taken into account at 
a distance of up to three distances to the nearest neigh- 
bors. The percolation path is traced by the broken line. 

problem. Quite likely, now it may be more dif- 
ficult for the reader to point out the difference 
between the site and circle problems than to no- 
tice the similarity. But the difference is there, 
and it is an important one. The fact is that in 
the circle problem the points that can serve as the 
centers of circles can lie anywhere on the plane, 
while in the site problem these are only the sites 
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of the lattice in question (see Fig. 33). If the 
number of sites within a circle is small, the differ- 
ence between the two problems becomes very 
important. It is then natural that the critical 
value of Zx s is a function of lattice type. The to- 
tal number of sites lying within a circle is Z. As 
we see from Table 5, Zx s indeed varies among 
lattices when Z is not too large. 

However, if Z is large, the difference between 
the two problems fades away. Imagine that we 
have started with the circle problem, but then 
shifted the center of each circle to the lattice 
site nearest to it. This is already a site problem. 
If the number of lattice sites within a circle is 
very large, this shift will not, with a high prob- 
ability, connect the circles that were not con- 
nected before the shift, and vice versa. This argu- 
ment shows that the site and circle problems 
become equivalent for Z oo. The value of B c 
determined by formula (1) is independent of lat- 
tice type and coincides with the quantity B c de- 
fined in the circle problem. 

This chain of arguments can be carried over 
entirely to the 3D case. The quantity B c defined 
by formula (1) coincides for 3D lattices with B c 
defined in the sphere problem. 

To recapitulate, the site problem on any plane 
lattice is reducible at large Z to the circle prob- 
lem, and the site problem on a 3D lattice is re- 
ducible to the sphere problem. The limit in formu- 
la (1) is, therefore, indeed independent of lattice 
type but depends on the dimensionality of the 
space on which the problem has been formulated. 
(The values of B c are different for the circle and 
sphere problems!) 
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Chapter 8 

Electric Conduction 
of Impurity Semiconductors. 

The Sphere Problem 

Percolation theory proved extremely useful for 
the understanding of the processes in doped se- 
miconductors (called extrinsic semiconductors to 
distinguish them from pure ones, called intrinsic 
semiconductors ) . It can rightly be said that at the 
present time impurity semiconductors represent 
one of the main fields to which percolation theory 
is applied. Several chapters of this book deal 
with impurity semiconductors. The present chap- 
ter begins with an exposition of current notions 
of semiconductor science, briefly interrupting the 
story of percolation theory. 

Intrinsic Semiconductors 

Let us start to discuss pure semiconductors, choos- 
ing as a basis semiconducting elements of the 
fourth group of the periodic table, such as ger- 
manium and silicon. These elements crystallize 
into diamond lattice (see Fig. 30a, b ). The four 
electrons making up the outer shell of each atom 
form bonds with the four nearest neighbors. The 
electrons move in such a way that their density 
is distributed nonuniformly around the atoms, 
being concentrated in four “strands stretching 
from the centers of tetrahedrons, in which the 
atoms are located, to the vertices of these tetra- 
hedrons. It is these strands that bind up the atoms 
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| of the crystal and do not let them fly apart. 

It can be imagined, by way of illustration, that 
each atom has four hands and holds hands with 
its four neighbors. Each electron is then strongly 
| hound to a neighbor, so that a moderately strong 

electric field applied to the semiconductor can- 
not generate an electric current: the number of 
i electrons is exactly equal to that necessary to 

form the bonds, while the energy required to re- 
lease a bonding electron is quite high. This 
energy, or rather, the minimum work that has to 
be done to transfer an electron from the bound 
state to the state in which it can move freely 
through the crystal is called the energy gap width 
(or simply gap width) and is denoted by E g . 

The gap width is an extremely important char- 
acteristic of a semiconductor, determining to a 
large extent all its electrical properties. For in- 
stance, let us consider electric conduction. The 
electric conduction in pure (intrinsic) semiconduc- 
tors exists because a certain fraction of bonds 
are broken. A bond breakdown releases an elec- 
tron ready to carry electric current and a “crip- 
pled”, “three-armed” atom. This “cripple” is also 
a charge carrier, only its sign is opposite to that 
of the electron. Indeed, an electric field can 
transfer an electron from a neighbor to the “three- 
armed” atom, changing its energy only slightly. 
As a result, another atom becomes “three-armed”. 
This process will continue, and you easily un- 
derstand that in contrast to electrons moving in 
the field, say, from left to right, a “three-armed” 
atom moves from right to left. (It is important 
i to realize that the atoms themselves do not 
move. In fact, a moving electron turns the atom 
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which it abandons into a “three-armed” atom.) For 
this reason, a “three-armed” atom can also be 
regarded as a charge carrier but, in contrast to the 
electron, it must be assigned the positive, not 
negative, charge. 

The microscopic picture is somewhat different 
in other semiconductors consisting of atoms with 
a different valence, but in all cases the breakdown 
of a bond produces two charge carriers with oppo- 
site signs. One of them (the negative) is called 
the electron, and the other (the positive) is cal- 
led the hole. A “three-armed” atom is thus a par- 
ticular case of hole. 

The energy required to break down a bond is 
supplied from the energy of thermal motion of 
the atoms. As it is proved in statistical physics, 
at high temperatures the mean energy of thermal 
motion per atom that performs small-amplitude 
vibrations around the equilibrium position is 
3 kT, where T is temperature in K, and k is the 
Boltzmann constant ( k = 1.38- 10" 16 erg/K). 

As a rule, the gap width E g is measured in elec- 
tron volts (eV). One electron volt is the work 
done by an electron traversing a potential differ- 
ence of IV: leV-1.6.10- 19 J - 1.6-10" 12 erg. 
In germanium E g = 0.7 eV, and in silicon E g = 
= 1.1 eV. At room temperature (300 K) the 
energy 3 kT is only 0.08 eV. This is much less 
than necessary for a bond breakdown both in 
germanium and silicon. 

However, thermal motion is chaotic, and at 
some random moments the energy of motion of an 
atom may be very high. This happens extremely 
rarely to each individual atom, but the total 
number of atoms is very large. There are roughly 
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10 22 atoms in 1 cm 3 . (The reader may remember 
that the number of atoms per mole equals Avo- 
gadro’s number, ^6-10 23 .) 

For this reason, the concentration of electrons 
and holes released by broken bonds is much smal- 
ler than the concentration of atoms; nevertheless, 
it is so high that it can result in an appreciable 
electric conductivity in a certain temperature 
range. Calculations show that in germanium the 
concentration of electrons is about 10 13 cm -3 at 
T — 300 K. In silicon, where the gap is wider, 
the concentration at the same temperature is 
much lower (10 10 cm -3 ). As temperature is reduced, 
the concentration of electrons decreases sharply, 
and the resistance of an intrinsic semiconductor 
undergoes a corresponding sharp enhancement. 

The number of charge carriers at a given tem- 
perature is the smaller, the wider the energy gap 
is. The difference between a semiconductor and 
a dielectric is merely the width of the gap. As a 
rule, materials with Eg of about 5 eV or more are 
referred to as dielectrics. They contain practical- 
ly no charge carriers at room temperature. 

While the difference between dielectrics and 
semiconductors is more quantitative than quali- 
tative, the difference between metals and dielec- 
trics is a principal one. The energy gap of a me- 
tal is zero, and the concentration of charge car- 
riers is high even in the immediate vicinity of 
the absolute zero of the temperature scale. 

Impurity Semiconductors 

Let us assume that germanium or silicon is doped 
with an impurity element from the fifth group of 
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the periodic table, such as phosphorus, antimony, 
or arsenic. An atom of these elements has five 
electrons on the outer shell. If such an atom sub- 
stitutes, for example, a germanium atom located 
at the center of a tetrahedron, four of its elec- 
trons go to form the bonds with the^four neigh- 
bors and the fifth electron remains dangling. 
This electron is bound to its atom because if the 
electron moves away, the atom becomes positive- 
ly charged and attracts the electron back. There- 
fore, the electron can be removed to a considerable 
distance away from the atom only by doing a work 
against attractive forces. This work is what we call 
the energy of bonding of the electron to an atom. 

As we will show below, the bonding energy of 
this extra electron is relatively low, so that the 
distance between this electron and the impurity 
atom is rather large (in comparison with the period 
of the lattice). As a result, the structure of the 
impurity atom resembles that of the simplest 
among atoms, namely, the hydrogen atom. We 
remind the reader that the hydrogen atom con- 
sists of a positively charged heavy nucleus and a 
light negative electron, the size of the nucleus 
being negligibly small in comparison with the 
distance from the nucleus to the electron. 

An impurity atom has a similar structure. The 
unit that acts as the nucleus of a hydrogen atom 
is not just the nucleus of the impurity atom but 
includes the inner-shell electrons and the four 
electrons that form the bonds. We shall be able 
to show now that the size of this complex is small 
compared with the distance to the extra elec- 
tron; the charge of the complex is positive and 
equals in magnitude the electron charge. 
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The electron in a hydrogen atom is known to be 
at a distance of the Bohr radius a B from the 
nucleus: 

= -^2 « 0.53- 10 -8 cm (1) 

Here h = 1.05- 10 -27 erg-s is Planck’s constant 
(divided by 2n), m = 9.8- 10~ 28 g is the electron 
mass, and e = 4.8 -lO -10 GGSE units is the elec- 
tron charge. 

The fundamental constant h was introduced by 
the German physicist Max Planck in 1901 in 
connection with the hypothesis on the quantum 
nature of electromagnetic waves. This constant 
enters into the quantum-mechanical equation 
describing the motion of an electron around the 
nucleus. The only other constants in this equation 
are e and m. The Bohr radius a B is the only pos- 
sible quantity with the dimensionality of length 
that can be composed of e, m, and H. 

The bonding energy of the electron in the hyd- 
rogen atom equals 

*b=-£~ 13.6 eV (2) 

The following interpretation of this formula by 
the order of magnitude is possible. The positively 
charged nucleus creates a potential e/r at a dis- 
tance r. At a distance of order a B it equals ela B , 
and at infinity it is zero. Consequently, the work 
that has to be done to transfer an electron located 
at a distance of order a B from the nucleus to an 
infinitely large distance equals, by the order of 
magnitude, e 2 /a B (the potential difference times 
the electron charge). 
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f Let us return to the impurity atom. The for- 
| mulas written for the hydrogen atom have to be 
5 modified in order to take into account that the 
extra electron moves not in a vacuum but in a 
semiconductor crystal. This modifies the form of 
Coulomb’s law. The force acting on an electron 
separated by a distance r from the nucleus is 
now e 2 /er 2 (in CGSE units), where e is the dielec- 
tric constant, or permittivity, of the crystal, 
and not e 2 /r 2 . This complication can be taken 
into account by replacing the quantity e 2 by 
e 2 /e in the expressions for a B and As- Further- 
more, it should be borne in mind that the mass 
describing the motion of an electron through the 
crystal does not equal that of a free electron, rn. 

The point is that the nuclei of the semiconduc- 
tor’s atoms and the inner-shell electrons produce 
in the crystal a periodic electric potential. One 
of the most interesting conclusions of the quantum 
theory of solid state is that the electron very 
nearly “overlooks” this potential if it is exactly 
periodic. The qualification “very nearly” consists 
in that the neglection of the periodic potential 
in the equations of motion of the electron must be 
simultaneously accompanied by replacing the elec- 
tron mass m by a mass m* which is a function 
of crystal properties. The quantity m* is called 
the effective mass. 

An impurity atom produces a nonperiodic po- 
tential that can by no means be omitted. But, 
when the motion of the electron in this potential 
is described, we can neglect the periodic potential 
of the crystal provided we replace m by m*. 

The equation describing the motion of the extra 
electron in a crystal around a charged impurity 
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I atom thus differs from the equation of motion of 
an electron in the hydrogen atom by the re- 
placements e 2 -> e 2 /e and m m*. 

Let us denote the characteristic distance to 
which the extra electron is removed from the im- 
purity atom by ag, and the bonding energy of 
this electron by E%. Making use of formulas (1) 
and (2), we obtain 

| <ib = 0.53-10- 8 e [cm] (3) 

I (4) 

As a rule, the effective masses in semiconduc- 
tors are substantially smaller than the mass of 
the free electron, and the dielectric constant 
equals 10-20 (e.g. in germanium m* ^0.1 m 
and s = 16). Consequently, the characteristic 
distance at which the extra electron is found 
equals, in typical semiconductors, from several 
tens to several hundreds of angstroms (1 A = 
= 10 ~ 8 cm), which is much greater than the in- 
teratomic distance (e.g. in germanium this spa- 
cing is 2.45 A). 

We have thus obtained that the extra electron 
in an impurity atom is removed to a considerable 
distance from the atom and is kept at this distance 
by the attractive forces caused by the nucleus. 
The fact that the distance is large is already 
a sufficient indication of smallness of the bonding 
energy E g. Indeed, the data cited above show 
that in typical semiconductors the bonding energy 
E% is from several hundredths to several thous- 
andths of an electron volt; hence, it is much 
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I smaller than the gap width E g (e.g. in germanium 

__ o 7 eV). 

^This is natural because it is easier to break 
loose the extra electron off the impurity atom 
than the electron bonding the host atoms. The 
thermal energy kT reaches the level of 0.01 eV 
at a temperature of the order of 100 K. As 
a rule, at this temperature a considerable frac- 
tion of extra electrons separate from their impu- 
rity atoms and take part in the transfer of elec- 
tric current. Therefore, impurity atoms belong- 
ing to the fifth group of the periodic table 
let go off their extra electrons rather easily. 
For this reason, they were given the name 
donors. 

Suppose that the impurity atoms belong to an 
element of the third group of the periodic table, 
for example, boron, aluminium, gallium, or others. 
These atoms have three electrons on the outer 
shell, so that they lack one electron for forming 
the bonds with their four neighbors. This elec- 
tron is readily borrowed from the neighboring 
atoms of the semiconductor, but then one of the 
neighbors becomes “three-armed” or, in other 
words, a hole appears in the neighborhood of the 
impurity atom. The impurity atom captures a 
fourth electron and thereby becomes negatively 
charged. The hole is thus attracted to this atom 
by electric forces, and a work has to be done to 
pull the hole away. This work is called the bond- 
ing energy between the hole and the impurity 
atom. 

The calculation of the bonding energy of the 
hole again leads to the problem of a hydrogen-like 
atom, but in this last case it is a positively charged 
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hole that moves around a negatively charged 
nucleus. The role of the nucleus is now played by 
an impurity atom together with the captured ad- 
ditional electron. This additional electron turns 
into a strand forming an interatomic bond, so 
that the region in which it moves does not exceed 
one interatomic distance. At the same time, a 
hole is bonded much less strongly to an impu- 
rity atom. A distance from the atom to the hole, 
and the bonding energy, are given by formulas (3) 
and (4). We only have to take into account that 
the effective mass of the hole that must be used 
in this case differs, generally speaking, from the 
electron effective mass. As a rule, this mass is 
also considerably smaller than the free electron 
mass, so that a hydrogen-like atom with a hole 
also has a size of the order of tens of angstroms, 
and its bonding energy is of the order of several 
hundredths of one electron volt. 

At temperatures of about 100 K the thermal 
motion breaks the bonds between holes and 
impurity atoms, after which the holes begin an 
“independent way of life” and, when an electric 
field is applied, participate in the electric cur- 
rent. 

Impurity atoms from the third group of the 
periodic table thus readily accept an electron 
and form a hole; hence, the term acceptors for 
such dopants. 

Let us summarize this section. Mobile charge 
carriers are formed in semiconductors only at the 
expense of the thermal motion energy. They can 
be formed when bonds of the lattice are broken. 
This requires that a work be done equal to the 
energy gap width E g . The outcome is the simul- 
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taneous appearance of an electron and a hole. An 
electron or a hole can also be born individually 
when their bond to an impurity atom is broken. 
The bonding energy between an electron or a hole 
and an impurity atom being much smaller than 
the gap width, the probability for the electron 
to break away from the impurity atom is much 
higher than the probability of breaking a lattice 
bond. On the other hand, the concentration of 
impurity atoms is usually lower by many orders 
of magnitude than the concentration of lattice 
atoms. Consequently, when temperature is in- 
creased, first of all charge carriers on impurity 
atoms are separated, the concentration of intrin- 
sic carriers (electrons and holes born of broken 
lattice bonds) being negligibly small. This tem- 
perature range is called the range of extrinsic 
conduction. However, as the temperature is 
further increased, the concentration of intrinsic 
charge carriers becomes comparable with that of 
impurity atoms and then grows larger than it. 
This range is called the range of intrinsic conduc- 
tion. 

The conclusion which is the most important for 
further exposition is that at a very low tempera- 
ture at which the thermal motion energy kT is 
small in comparison with the bonding energy of 
electrons to impurity atoms, the semiconductor 
has neither extrinsic nor intrinsic charge carriers. 
The bonds are intact where they ought to be, and 
impurity electrons and holes are localized in the 
vicinity of their atoms. In this range, electric 
conduction of the semiconductor drops very 
sharply to zero as its temperature is lowered. 
This behavior distinguishes a semiconductor 
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from a metal in which the concentration of mobile 
charge carriers remains high no matter how 
small the temperature is. 

Transition to Metallic Electric Conduction 
at Increased Impurity Concentrations 

This principal difference between semiconductors 
and metals suddenly disappears when impurity 
concentration is increased. If concentration goes 
beyond a certain critical value N c , electric 
conductivity remains relatively high and weakly 
dependent on temperature no matter how low the 
temperature is. This electric conduction is refer- 
red to as metallic. This does not mean that the 
electric conductivity of a semiconductor becomes 
comparable with that of good metals. Far from 
it: the electric conductivity of a semiconductor 
always stays many orders of magnitude lower. 
The term only reflects the behavior of electric 
conductivity at low temperatures. The transition 
to metallic electric conduction that takes place at 
an increased impurity concentration is called the 
metal-insulator transition, or the Mott transition 
(after the famous British physicist Sir Nevill 
Mott). 

Experiments have demonstrated that the criti- 
cal impurity concentration N c at which the Mott 
transition takes place varies quite considerably 
among different semiconductors. A reliable esti- 
mate of the quantity N c is given by a relation 

N c a B ~ 0.02 (5) 

For instance, in germanium N c = 10 17 , and in 
silicon N c = 3-10 18 cm -3 . 
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Even now we lack a satisfactory mathematical 
theory of the Mott transition in semiconductors. 
In fact, this is one of the most complicated prob- 
lems in solid state theory. The qualitative picture 
looks as follows. 

The structure of impurity atoms resembles that 
of univalent elements of the first group of the 
periodic table (H, Li, Na, K). Like these ele- 
ments, an impurity atom has a single electron on 
its outer shell. When crystallized, elements of 
the first group form good metals. Hydrogen forms 
a molecular crystal which is an insulator. How- 
ever, there is every reason to believe that this 
insulator also converts to a metal at a sufficient- 
ly high pressure. (The research in this field is 
conducted nowadays with such vigor that it could 
make the topic for a separate book.) 

It appears as nearly obvious that if hydrogen- 
like impurity atoms are distributed in a semi- 
conductor with a sufficiently high density, they 
will also form a metallic system. 

It might seem strange at first glance that this 
argument has a bearing on the Mott transition. 
Indeed, the concentration of atoms in a metallic 
sodium is of the order of 10 2 “ cm 3 , that is, great- 
er than the concentration of impurities in ger- 
manium, at which the Mott transition occurs, 
by a factor of about 100 000. Obviously, sodium 
cannot be a metal at such low concentra- 
tions. 

However, it is not just the concentration but 
the degree to which the electron shells of neigh- 
boring atoms overlap that determines whether a 
system is metallic or not. If the atoms are at 
such a large distance that the regions of space in 
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which their outer-shell electrons move are far 
from one another, these are simply isolated atoms. 

But if the nuclei of the atoms are sufficiently 
close for these regions to overlap, the atoms lose 
their individuality. Their outer-shell electrons 
“cannot know” to which of the nuclei they be- 
long. A system of a large number of atoms appro- 
priates, or collectivizes, electrons. These electrons 
form a separate system capable of conducting an 
electric current. This material is a metal. 

The degree of overlapping is determined for \ 
hydrogen-like atoms by a dimensionless parame- 
ter Na%\ (The quantity (4/3) jt Na% is the mean 
number of atomic nuclei within a sphere of radius 
equal to the effective Bohr radius.) It has already 
been mentioned that since an impurity atom is 
located inside a semiconductor, its radius is ano- 
malously large. For this reason, condition (5) is 
met already at N = 10 17 cm -3 (for germanium). 

It must be clear by now from the arguments 
given above why the Mott transition occurs at ap- 
proximately the same value of the parameter 
Na% in all conductors, although critical con- 
centrations N c may vary quite considerably. 
Indeed, it is this parameter that determines the 
overlapping of neighboring atoms! 

The following question can be posed: What 
should be the concentration of atoms in a hydrogen 
crystal for the overlapping of its atoms to be the 
same as the overlapping of impurity atoms in a 
semiconductor at the concentration corresponding 
to the Mott transition? The overlapping is given ) 
by formula (5), but now it is natural to substitute 
into it the Bohr radius found from formula (1). 

This will give an enormously high concentration: 
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10 23 cm" 3 . Hence, the overlapping corresponding 
to the Mott transition is extremely large. 

There is no ground, therefore, to be surprised 
that at a concentration greater than N c impu- 
rity atoms form a metallic system. This fact 
lies at the foundation of the modern theory of 
semiconductors with a high impurity concen- 
tration. 

The Mott Transition and Sphere Problem 

Percolation theory offers a simplified description 
of the Mott transition which takes into account 
that impurities are located in a semiconductor 
randomly and can form tighter or looser ensem- 
bles. Assume that impurity atoms in a certain re- 
gion are so close to one another that their electron 
shells strongly overlap and the outer-shell elec- 
trons become “collective property”. This region 
is a piece of metal: a potential difference applied 
to this region would cause an electric cur- 
rent. 

The presence of such regions is not sufficient, 
however, for a large sample to behave as a me- 
tal. If metallic regions are infrequent, they are 
not in contact, forming only isolated metallic 
islands in the material that at low temperatures 
behaves as a dielectric. Taken as a whole, the 
combination is equally dielectric. 

As impurity concentration increases, the frac- 
tion of space occupied by metallic regions grows, 
and at a certain critical concentration N c the 
metallic regions form a connected system of “lakes 
and canals” that permeates the whole semiconduc- 
tor crystal. Beginning with this concentration, 
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the electric conduction of a large sample be- 
comes metallic. 

Obviously, the concepts presented above must 
be formulated mathematically in terms of perco- 
lation theory, although the formulation will not 
come easily. The main difficulty is that we do 
not know the impurity concentration at which a 
region can be regarded as metallic. 

The simplest model of the Mott transition 
enunciated in the sixties reads as follows. Ima- 
gine that each impurity atom is a metallic ball 





Fig. 34. (a) Overlapping spheres representing atoms with 
collectivized electrons; (b) a string of overlapping spheres 
forming a metallic channel through which electric current 
can flow. 

of radius r 0 . The balls can penetrate into one 
another (Fig. 34a), that is, the regions of space 
that they occupy can overlap. The balls form 
strings (Fig. 34 b) and regions of more complex 
shape that are regarded as metallic by definition. 
We want to find the concentration N c of balls 
beginning with which the metal regions make the 
whole large sample electrically conductive. 

But how to choose the radius r 0 ? A hydrogen- 
like impurity atom does not have well-defined 
boundaries. The probability of finding an electron 
at a distance r from the nucleus diminishes with 
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Increasing r by the law exp (— 2r/ag), where 
: e = 2.718 is the base of natural logarithm. The 
probability of finding an electron at r = ag/2 is 
less by a factor of 2.7 than that at r = 0; at r = 
«= ag it is less by a factor of 7.4 than that at r = 
= 0; at r = 1.5 ag it is smaller by a factor of 20. 
The effective radius of the atom must definitely 
be proportional to the length ag: 

r 0 =qa% (6) 

This is an important statement. It signifies 
that the numerical coefficient q must be fairly 
universal: it varies rather weakly in going from 
one semiconductor to another, while the length 
ag varies quite considerably. 

It is very difficult to calculate q on the basis of 
physical arguments. A simpler approach is to find 
the critical concentration N c corresponding to 
the Mott transition in terms of the model for- 
mulated above. This will give us N c as a func- 
tion of the ball of radius r 0 . Then, making use 
of formula (6), we must express N c in terms of q 
and ag. According to experimental data, N c is 
given by formula (5). Comparing the theoretical 
expression with formula (5), we can find g, that 
is, determine it “as if” from experimental data. 

Let us start the realization of this program. 
We have to solve the following problem of perco- 
lation theory. Spheres of radius r 0 , whose cen- 
ters are distributed in space randomly and on the 
average uniformly, are drawn in the three-di- 
mensional space. The mean number of centers of 
spheres per unit volume is N. Two spheres are 
regarded as connected if they overlap (see 
Fig. 34 b). We need to find the critical concentra- 
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tion N c at which percolation sets in through 
overlapping spheres, that is, paths appear that 
traverse the whole system and are composed of 
overlapping spheres (see Fig. 34 b). 

This problem differs from the sphere problem 
formulated in the preceding chapter in that the 
connected spheres in the sphere problem were not 




Fig. 35. String of inclusive circles of radius 2 r 0 is shown 
by the dashed lines. The centers of circles are marked by 
1 the crosses. The percolation path through these centers 
is shown by the broken line. The percolation path through 
overlapping circles of radius r 0 follows the same 
centers. 

merely overlapping but inclusive, that is, their 
centers were at a distance smaller than r 0 , not 
than 2r 0 . However, this difference is not too 
significant, and the results of one problem are 
easily transferred to the other. Indeed, if an in- 
finite path exists through inclusive spheres of 
radius 2r 0 at a certain concentration of centers 
i of spheres, there also exists a path through over- 
lapping spheres of radius r 0 . This is just the 
same path, that is, the path through the same 
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centers. This argument is illustrated in Fig. 35 
for the plane problem. (The reader will easily 
visualize the corresponding drawing for the 3D 
case.) If there is no path through inclusive spheres 
of radius 2 r 0 , there is no path through overlap- 
ping spheres of radius r 0 . Hence, the critical con- 
centration corresponding to the percolation thre- 
shold through inclusive sphere of radius 2r 0 equals 
the critical concentration corresponding to the 
percolation threshold through overlapping spheres 
of radius r 0 . 

According to the results given in the preceding 
chapter, the critical concentration for percolation 
through inclusive spheres of radius 2r 0 is deter- 
mined by the condition 

(4/3) ji N 0 (2r 0 ) 3 = B c « 2.7 (7) 



Substituting (6) into (7), we obtain 




(8) 



The coefficient q can now be found by comparing 
(8) with expression (5) obtained from experimen- 
tal data. In order for these expressions to com- 
cide, we have to set q = 1.6. Therefore, the effec- 
tive radius r 0 equals 1.6ag. The probability to 
find an electron at the point at a distance ot 
r = r 0 from the nucleus is less by a factor of 24 
than the probability to find it at the point r == 0. 

The most important feature of this distribution 
of the metal-insulator transition lies not so much 
in calculating the quantity q that determines the 
effective atomic radius but, of course, in gaining 
an insight, through percolation theory, into the 
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internal organization of semiconductors at impu- 
rity concentrations close to N c . 

If the suggested description is correct, a semi- 
conductor contains a system of metallic channels 
that run across the matrix. Electric current flows 
through these channels as through wires. In 
Part III of this book it will be shown that the 
fraction of volume of a semiconductor occupied by 
these channels is very small if impurity concen- 
tration is close to N c . This entails specific prop- 
erties of the electric conduction and other im- 
portant characteristics of semiconductors. 

Exercise 

1. Indium antimonide (InSb) is a semiconduc- 
tor with a very narrow gap (0.18 eV at room 
temperature). Effective masses are also small in 
such semiconductors. Assuming that the effective 
electron mass m* is 0.015 m and the dielectric 
constant is e = 18, calculate the effective Bohr 
radius a$ and the critical concentration N 0 cor- 
responding to the Mott metal-insulator transition. 



Chapter 9 

V arious Generalizations 
of the Sphere Problem 

Inclusive Figures of Arbitrary Shape 

We have already mentioned that the motion of 
electrons through a crystal is described by the 
effective mass m *; this m* may considerably 



f Ch. 9. Generalizations of the Sphere Problem 167 

differ from the free electron mass, that is, the 
mass of an electron in an empty space. It was 
also found that the effective mass may differ from 
the free electron mass not only in magnitude. 
The point is that not all directions in a crystal 
are equivalent. For instance, the motion of an 
electron along the edges of a cube may differ 
from the motion along body diagonals and from 
the motion along face diagonals of the cube. 
Consequently, the effective mass is not necessari- 
ly identical in different directions. As a result, 
the region in which an electron moves in the 
neighborhood of a donor atom will not be spheri- 
cal. It may be ellipsoidal or even more complex. 

These arguments explain why the sphere prob- 
lem was generalized to the case of figures of arbi- 
trary shape. The new problem is formulated as 
follows: the sites (centers) are distributed in 
space randomly and on the average uniformly. 
The site concentration is N. The sites are sur- 
rounded by identical closed surfaces of arbitrary 
shape. 

The surfaces surrounding different sites are iden- 
tical not only in shape but also in surface orien- 
tation in space. If, for instance, the surfaces are 
fish-shaped, the tails of all the fishes must face in 
the same direction. 

The volume within one surface is V. Two sites 
are said to be connected if one of them lies with- 
in the surface drawn around the other site (inclu- 
sive surfaces). We have to find the critical con- 
centration N c at which percolation through con- 
nected sites sets in. 

As in the sphere problem, the presence or ab- 
sence of percolation is determined only by the 
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value of the parameter B which equals the pro- 
duct VN. This is the mean number of sites within 
the volume bounded by one surface. This number 
will not change if we enlarge or contract the 
scales along all the three directions, that is, if we 
multiply or divide the coordinates of all sites in 
the system and of all points on the surfaces by 
the same factor. It is also clear that this rescal- 
ing will not suppress percolation if it existed; 
neither will it bring it to existence. Therefore, 
percolation “does not respond” to those changes 
in N andF that leave B unaltered. Consequently, 
it will be more convenient, as in the sphere prob- 
lem, to argue not in terms of the critical concen- 
tration N c but in terms of the critical parameter 



If we keep the shape of the surface drawn around 
each site constant but increase the volume V 
bounded by this surface, say, twice, the critical 
concentration N c will be half as large, while the 
parameter B c will not be effected. It depends 
only on the surface shape . 

Formula (1) generalizes the formula 



B c = (4/3) jt N C R 3 



that we used earlier for the sphere problem. 

In the general case, the critical value of B c at 
which percolation sets in must not be equal to 
the value of B c = 2.7 obtained for the sphere 
problem. At present the value of B c for different 
figures is a subject of intensive study. 
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The Ellipsoid Problem 

The Soviet scientists Ya. G. Sinai and B. I. 
Shklovskii proved that certain distinct surfaces 
have the same value of B c . Thus, the ellipsoid an 
sphere (and the ellipse and circle m the plane 
paspI are such surfaces. ■ _ . 

First we remind the reader of the definitions of 
ellipse and ellipsoid. The ellipse is a closed curve 




Fig. 36. Ellipse E obtained from circle 0 by extending 
it along the z-axis. 

drawn on a plane obtained from a circle by extend- 
ing (or contracting) it along one of its axes 
(Fig. 36). In order to realize this extension or 
contraction, it is necessary to transform each 
point M of the circle with x-, {/-coordinates ^nto 
the point M' with coordinates x' = kx, y — y, 
where k is the extension factor (k > 1 corresponds 
to extension, and k < 1 to contraction). 

The surface obtained from a sphere by exten- 
sion (or contraction) along one of its axes is 
called the ellipsoid of revolution. In order to realize 
this extension or contraction along the z-axis 
it is necessary to transform each point on the 
sphere with x-, {/-, z-coordinates into the point 
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with coordinates x' = x, y f = y, 2 ' = kz. The 
figure thus obtained is shown in Fig. 37. It is 
called the ellipsoid of revolution because a rota- 
tion by any angle around the 2 -axis (called the 




Fig. 37. Ellipsoid of revolution. 

axis of revolution) transforms this solid into 
itself. If it is cut by a plane drawn through the 
2 -axis the section is an ellipse. The section by a 
plane drawn perpendicularly to the 2 -axis is a 
circle. 

The general-type ellipsoid is obtained from an 
ellipsoid of revolution by extending (or contract- 
ing) it along one of its axes perpendicular to 
the axis of revolution. The new extension factor 
is not necessarily equal to k. The section of the 
general-type ellipsoid by any plane parallel to 
the xOy- } xOz- } y 02 -planes is an ellipse. 
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Now we shall prove that ellipsoids and spheres 
have identical B c . Let the sites be randomly 
distributed in space with the mean concentration 
iV. Draw a sphere of radius R around each site. 
If B = (4/3) nNR 3 > 5 s c ph there is percolation, 
otherwise, there is no percolation through con- 
tacting spheres. Here B s c ph is the critical value 
of B for the sphere problem (B S <F* = 2.7 ± 0.1). 

Let us apply the extension y' = y, x = k x x , 
z’ = k 2 z along the 2 - and £-axes in order to 
transform the coordinates of both the sites and 
points on the spheres. (If we dealt with a plane 
problem, the picture could be modelled by mark- 
ing the sites and circles on a rubber band which 
is then extended in one direction. Likewise, we 
could imagine a three-dimensional “rubber space” 
in which we have marked the sites and spherical 
surfaces. Then this “rubber space” is extended 
in two directions.) 

If the ^-coordinates of sites are random numbers 
uniformly distributed within the interval from 
0 to L, where L is the size of the system, the new 
^'-coordinates obtained by multiplying x by k x 
are also random numbers uniformly distributed 
within the range from 0 to k x L. The 2 '-coordinates 
are uniformly distributed within the interval 
from 0 to k 2 L, while the ^-coordinates remain 
unaltered. The concentration N ' of sites is not 
equal to N. 

All spheres now transformed into ellipsoids 
with volume V'. (It can be shown that V' = 
= k,k 9 -(i/3)nR s . This is irrelevant for the 
derivation to follow.) As a result, the sphere 
problem converted to the ellipsoid problem. Here 
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we can introduce a quantity B f = N'V' and 
find the critical value of B ' at which percolation 
through ellipsoids sets in. Let us denote it by 5® 11 . 

The rest of the proof separates into the following 
steps: 

1. B = B '. All sites that lay within a certain 
sphere prior to transformation are found inside 
the ellipsoid obtained from this sphere after 
the transformation. Indeed, when the sphere is 
extended (or contracted), the points internal to 
the sphere remain internal, and those external to 
it remain external at each step of extension. 
Hence, the mean number of sites, B , which fell 
inside one sphere prior to transformation equals 
the mean number of sites, S', which fell inside 
one ellipsoid after the transformation. 

2. If B > £c ph , then B>B e c 11 (3) 

If B<Bl ph , then B<Bl n (4) 

Indeed, if two spheres were connected prior to 
transformation, that is, the center of one of them 
was inside another, the two ellipsoids formed 
out of these spheres are also connected because 
the transformation will leave the internal points 
within, and the external points without, these 
surfaces. If two spheres were not connected, the 
ellipsoids obtained from them are not connected 
either. This yields that if B > 2?c Ph , that 
if infinite percolation paths exist through con- 
nected spheres, infinite percolation paths also 
exist through connected ellipsoids, and this sig- 
nifies, in turn, that B f is greater than the threshold 
value 5° 11 . This gives us condition (3) because 
B* = B. If B < 5c Ph , that is, no percolation 
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exists through connected spheres, there is no 
percolation through connected ellipsoids, that is, 
B’ < 5® 11 . This gives us condition (4). 

3. Since conditions (3) and (4) must hold for 
any value of 5, they imply that 5c Ph = B c , 
which was to be proved. 

Other Surfaces 

The ellipsoid problem is very important in semi- 
conductor physics, but in some cases the ellip- 
soid is too simple a surface. Physicists were there- 
fore much interested in how the critical value of 
B c depends on the shape of the surface when exact 
relations cannot be obtained. To gain this knowl- 
edge, Monte Carlo computer simulations were 
carried out, and the following surprising and 
extremely important fact was found out: B c is 
only very weakly dependent on the surface 
shape. 

The shapes that were analyzed were a cube and 
a tetrahedron. It was found that the critical 
values for these solids do not differ, within the 
error of calculations (±0.1), either from each 
other or from the value of B c for the sphere prob- 
lem. Attempts to find a solid with least resem- 
blance to a sphere led to a “three-dimensional 
cross”, that is, a figure formed by three elongated 
parallelepipeds intersecting at the origin of 
coordinates (Fig. 38). It was found that B c for 
this surface is only 20% less than that for the 
sphere. 

The same conclusion was obtained for plane 
figures. An analysis demonstrated that B c for 
squares differs from B c for circles by mere 2%. 
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Consequently, B c is, to a satisfactory accuracy, 
universal within a class of figures of identical 
dimensionality (i.e. either plane or three-dimen- 
sional). The reason for this observation may be 
a certain relation similar to that proved in the 
preceding section, so that exact equalities in 
fact hold where computer simulation discovers 
only a slight deviation or no deviation at all, 




Fig. 38. “Three-dimensional cross”. 

within computational errors. (Small deviations 
found by computer calculations must be looked 
at critically because it is not always easy to carry 
out a correct evaluation of results.) Unfortunately, 
this is all we know at present. 

Another Experiment at the House Kitchen. 

The Hard-Sphere Problem 

In 1974 three students of the Harvard University 
in the USA, J.B. Fitzpatrick, R.B. Malt, and 
F. Spaepen, conducted the following simple 
experiment. Five thousand small balls, part 
of them made of aluminium and another part of 
plastic, were placed in a jar. Before starting the 
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m. 

I measurements, the balls were thoroughly stirred 

I and then the jar was intensely shaken to make 
the packing of balls maximally tight. A metal 
foil electrode was placed at the bottom of the jar, 
and another electrode was placed on top of the 
balls, after which the system was compressed by 
a load of 15 kg. 

Aluminium is a highly conductive metal, and 
the plastic was an insulator. The aim of the 
experiment was to determine the critical fraction 
x c of aluminium balls at which electric current 
flows between the electrodes, that is, paths appear 
through the aluminium balls in contact. It was 
found that x c « 0.25. 

In addition, the experiment made it possible 
to study the electric conductivity of the system 
as a function of x at x > x c . 

Here we encounter another problem of per- 
colation theory, namely, the hard-sphere problem . 
The closest relation to this new problem is per- 
haps the site problem. Let us recall our approach 
to an approximate evaluation of percolation 
threshold in the site problem used in Chapter 6. 
A sphere is constructed around each lattice site, 
its radius being equal to half the distance to the 
nearest neighbor. The spheres constructed around 
white sites are said to be white, and those around 
black sites are said to be black. Percolation 
through white sites is equivalent to the existence 
of percolation paths through tangent white 
spheres (see Fig. 31). 

The difference between the site problem and 
the new problem lies in that the centers of spheres 
in the site problem are the sites of a regular lat- 
tice, while in the new problem the centers can be 
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located at any point. It will be shown below 
that this difference is not very important. 

The difference between the new problem and 
the sphere problem discussed in the preceding 
chapters and the preceding sections of this chap- 
ter lies in that the spheres that we operate with 
in the new problem are assumed to be hard and 
thus cannot overlap. This is an essential distinc- 



tion. 

It was shown in Chapter 6 that in the site 
problem percolation through white sites sets in 
when the fraction of space filled by the spheres 
constructed around these sites is approximately 
0.16. It was found that this number is almost 
independent of lattice type. It is then natural 
to conjecture that if it is nearly independent of 
lattice type, it should not strongly depend on 
whether the lattice is or is not there . If this hypo- 
thesis is correct, the fraction of volume filled 
with metallic balls at which percolation appears 
through these balls must be approximately 0.16. 

Let us denote by /, as in Chapter 6, the filling 
factor, that is, the fraction of volume occupied 
jointly by both aluminium and plastic balls. 
By definition, the quantity x is the ratio of the 
number of aluminium balls to the total number 
of balls. Hence, the fraction of volume occupied 
by aluminium balls equals fx. If this fraction 
of volume equals 0.16 at percolation threshold, 
the critical value x c can be found from the con- 
dition fx c = 0.16. 

The filling factor of a system of tightly packed 
but randomly distributed balls is well known. 
Such a system is familiar to mankind from the 
times of antiquity. If it was necessary to measure 
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a certain amount of grain or some other free- 
flowing material, it was poured into a special 
vessel (a measure), shaken, and compacted. In 
modern science this system is a model of atomic 
arrangement in amorphous metals. Amorphous 
metals, also referred to as metallic glasses, are 
materials with metallic electric conductance but 
having no crystal structure. It was found that 
the atomic arrangement of amorphous meta 
strongly resembles the arrangement of tightly 
packed incompressible balls. For this reason, 
the properties of tightly packed random y 
tributed balls were studied very thoroughly 
(mostly in computer models). ^^as found tha 
the fraction of volume filled with the halls 

1 Le°t' 6 us 7 ‘return to the problem of percolation 
through metallic bails. Having determined x c 
£ means of the formula ,x c = 0.16// we obtain 
Xo — 0.25, in complete agreement with the 

result of the three students. 

Experiments on determining the percolation 
threshold were often repeated, each 
better instruments. Experiments were conducted 
in which balls in the mixture were of different 
radii. The radii of both metallic and dielectric 
balls were varied in the mixture m a wide range. 
It was found that in this case the critical volume 
fraction fx c is approximately 0.17, bei »g al “ 10 ® 
the same, to within experimental accuracy, as 
in the case of identical balls. 

Percolation threshold in the hard-sphere prob- 
lem can thus be evaluated rather easily after 
we assume that the critical volume fraction filled 
with metallic balls is approximately 0.16. Note 



178 Part II. Various Problems of Percolation Theory 




that this fraction is much greater in the problem 
of overlapping spheres. 

The hard-sphere problem proved very impor- 
tant for applications as well. It accounts for the 
theory of heterogeneous materials manufactured 
from an insulator with tiny metallic inclusions. 
Such materials are now an important object of 
study. They are prepared and utilized both as 
thin films and as bulk samples. The electrical 
properties of such materials in the vicinity of 
percolation threshold are fascinating. For in- 
stance, the capacitance of a capacitor filled with 
such a material grows to infinity when the frac- 
tion of volume filled with the metal tends to per- 
colation threshold. This phenomenon stems from 
an enormously high mutual capacitance of large 
metallic clusters. The description of electrical 
properties of heterogeneous materials is growing 
nowadays into an independent domain of per- 
colation theory. 
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Chapter 10 

Percolation Level 



“The Flood” 



Six days and [six] nights 
Blows the flood wind, as the south-storm 
sweeps the land. 

When the seventh day arrived, 

The flood (-carrying) south-storm 
subsided in the battle, 

Which it had fought like an army. 

The sea grew quiet, the tempest 
was still, the flood ceased. 

I looked at the weather: stillness 
had set in, 

And all of mankind had returned to clay. 
The landscape was as level 
as a flat roof. 

I opened a hatch, and light fell upon 
my face. 

In: Ancient Near Eastern Texts 
Relating to the Old Testament » 
Ed. by J. B. Pritchard, 
Princeton Univ. Press, 1955, 
p. 94 

(Accadian Myths and Epics , 

the Epic of 

Gilgamesh, 

Gilgamesh, Tablet XI, trans- 
lator 

E. A. Spaiser) 

This description of the Flood is found in the 
Babylonian Epic of Gilgamesh, that most an- 
cient work of literature which dates back roughly 
to the second millennium BC. 



181 



180 Part II. Various Problems of Percolation Theory 

When the waters began to subside, the sum- 
mits of the highest mountains appeared above 
the surface. The water subsided lower and lower 
until it reached its normal level. Just imagine 
this vast picture: a huge system of mountain 




Fig. 39. Voyage during the Flood. 



ridges such as the Himalayas, gradually emerged 
irom the water. First the highest peaks broke 
the surface and formed islands, then the zone of 
alpine meadows was liberated, and finally, water 
dropped to foothills. 

Suppose that we want to answer the question: 
To what level must the water drop in order for 
the last waterway across the whole system of 
ridges to disappear (Fig. 39)? Clearly, such a 
way exists as long as a certain fraction of moun- 
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Fig. 40. Map of a mountain system flooded with water. 
The submerged regions are shaded. 

tainous passes remain under water, but then it 
dries up. 

Figure 40 illustrates a geographic map of a 
mountain system on which the thin lines trace 
the contour lines (curves of identical height). 
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The thicker lines mark the contours corresponding 
to water level. They separate water from dry 
land. The submerged regions are shaded in the 
drawing. Figure 40a shows the very beginning of 
the Flood: the lakes do not communicate with 
one another. Then the water level climbed, and 
Fig. 406 now shows the waterways across the 
system. And in Fig. 40c only individual summits 
are above the water. 

The level of water at which this waterway 
appears (or disappears) is called percolation level. 
The problem of finding this waterway is a plane 
problem of percolation theory. It can be refor- 
mulated. Let us assume that a plane is randomly 
painted by white and black paints. Let the frac- 
tion of area painted white be x. At small values 
of x, white spots form isolated islands, while at 
x nearly equal to unity, the isolated spots are 
black. We want to find the critical value of x 
at which a noninterrupted path appears (or dis- 
appears) across the whole system, going only 
through white regions. 

Likewise, we can also formulate the correspond- 
ing 3D problem, filling a large volume with a 
white and a black substances. Then we have to 
vary the fraction of volume occupied by one of 
these substances until percolation sets in. 

How to Construct a Random Function ** 

This new problem is not a lattice problem. Let 
us analyze its mathematical formulation in the 
plane case. Wo have to define on the whole plane 
a random function V ( X , Y), where X and Y 
are coordinates. In the Flood problem this func- 



tion is the altitude above sea level of the point 
on the earth’s surface with the given coordinates. 
We shall confine ourselves only to the case of 
Gaussian random functions (named after the great 
mathematician Karl Friedrich Gauss). 

The simplest method of constructing a Gaus- 
sian random function is to ascribe to each point 
in space a random number not related in any 
way to a neighboring random number. ^ Such 
random functions were given the name “white 
noise”. The values it takes at neighboring points 
differ abruptly, that is, this function is discon- 
tinuous. To obtain a continuous function, the 
white noise has to be “smoothed”. This proce- 
dure consists in ascribing to each point of space 
a quantity equal to the average of the values that 
the “white noise” function assumes in a certain 
region around this point. The continuous random 
function is then formed by these values. The 
function varies only slightly from this point to 
the neighboring point because of the small size 
of the spatial region over which the white noise 
is averaged. Let us denote by r 0 the size of the 
region over which averaging is carried out. This 
size is referred to as the correlation radius of a 
random function. Its main property is that the 
function changes only slightly when the argu- 
ment of the function changes by an increment 
small compared with r 0 . ,. 

The quantities V can be described by a dis- 
tribution function that we denote by / (T)- By 
definition, the probability for the function 
V (X Y) to have at a point taken randomly in 
space’ a value within a small interval from F x 
to V x + AF equals / (F x ) AF. 
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The function V (X, Y) can always be con- 
structed in such a way that, averaged over all 
points of a plane, it gives zero. To achieve this, 
we need to use in the construction of the “white 
noise” function the random numbers distributed 
symmetrically with respect to zero. It can be 
shown that the function V ( X , Y) obtained from 
such “white noise” has the Gaussian distribution 
(see Chapter 2) of the type 

/<v) = exp (~w) <<) 

Note, that the probability density (1) is sym- 
metric relative to the positive and negative val- 
ues of V. In “mountain lingo” this means that 
peaks and valleys are encountered with equal 
probability. 

Now we are ready to enunciate the percolation 
problem. The contours shown in Fig. 40 are de- 
termined by the condition V ( X , Y) = const. In 
order to prescribe the water level, we need to 
introduce a quantity t that can vary from — oo 
to +oo. The regions on the plane where V (X, Y) < 
< ^ will be referred to as white (covered with 
water), and those where V (X, Y) > t as black 
(protruding above water). The thick contours in 
Fig. 40 correspond to V (X, Y) = t. 

The percolation level is defined as the critical 
value t c at which white regions form paths 
going uninterruptedly across the whole system 
It is also possible to speak of the critical fraction 
x a of space filled with white regions at the moment 
when percolation sets in. This fraction of space 
equals the probability for a continuous random 
variable V to assume a value in the range — oo 
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< V < t c . By virtue of the definition of distribu- 
tion functions, 

*c 

* c = j f(V)dV (2) 

— 00 

Formula (2) relates the critical fraction x 0 of 
space to the percolation level f c . 

Analogy to the Site Problem ** 

Imagine that a plane lattice (its type is immate- 
rial) with small period is superposed over the 
geographic map shown in Fig. 40. Let the lattice 
period be substantially smaller than the charac- 
teristic extensions of land and water regions. 
Let us “paint” the lattice sites in water-covered 
areas white, and those in land areas black. As in 
the familiar site problem, we again consider 
white sites connected if they are nearest neigh- 
bors. The critical fraction of area under water 
equals the critical fraction of white sites at which 
percolation through white sets in. 

However, the new problem is not identical to 
the site problem. Remember Chapter 4 where 
the construction of a system of white and black 
sites has been described in detail. (Sites were 
labelled in that chapter as blocked and non- 
bloeked.) It was an important aspect of that 
construction that each site became black or 
while as prescribed by a random-number gene- 
rator, in no relation to the color of the sites in 
the nearest neighborhood. As a result, the white 
and black sites were thoroughly mixed. The 
situation is quite different in the new problem. 
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The lattice period being very small, black and 
white sites form large groups. Neighbors of a 
white site are almost certainly white, and those 
of a black site are most probably black. (The 
size of blocks is determined by the correlation 
radius r 0 .) 

It must be clear that to introduce a lattice into 
the percolation level problem is a purely formal 
operation. This is in fact a nonlattice problem, 
and its solutions must be independent of both 
the lattice period (provided it is sufficiently 
small) and lattice type. However, this method 
allows us to make use of a well-developed appa- 
ratus of lattice problems. 

In Chapter 4 we gave an algorithm of the Monte 
Carlo solution of a site problem with a computer. 
This algorithm is entirely transferable to the 
percolation level problem. An array K ( X , Y) 
consisting of zeros and unities is generated 
by addressing a random function V (X, Y) by 
the method described in Chapter 4. White sites 
correspond to unities, and black sites to zeros. 
Then percolation paths are searched for, and the 
critical fraction of space at which percolation 
sets in is determined by the same techniques. 

Percolation Levels 

in Plane and Three-Dimensional Problems ** 

A plane problem has an exact solution if the 
properties of a random function V (X, Y) are 
on the average symmetric with respect to V = 0. 
The Gaussian functions described above are 
among the functions with such properties. 

In order to arrive at the exact solution, we 
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must use the symmetric formulation of the 
percolation problem given in Chapter 5. As be- 
fore, we refer to the regions in which V (X ,Y) <C t 
as white, and to those in which V (X, Y) > t 
as black. In addition to percolation level through 
white regions, t c , we can introduce percolation 
level through black regions, f 0 '. The further proof 
goes by the following steps: 

1. By virtue of the symmetry of the function 
V (X, Y), we have f c = — C Indeed, replace V 
by —V at each point of the plane. This will 
give a function V' = —V that on the average 
has the same properties, so that percolation 
levels calculated by V' must be equal to those 
calculated by the original function V. The 
inequality V < t c yields the inequality — V > 

> - t c , that is, V' > -t 0 . At t = t c the regions 
that are white with respect to the function V 
(P < t 0 ) form an infinite cluster. These same 
regions of space are black with respect to the 
function V' and t = — t 0 . Indeed, here V' > 

> — t c . Therefore, the value t = — 1 0 is per- 
colation level through black regions for the 
function V r . But it has already been said that 
the functions V and V' should have identical 
percolation levels. Consequently, percolation level 
through black regions for the function V is 

t' — f 

C 2. If t 0 < 0, then, as t increases, first percola- 
tion through white regions appears (at t = 
= t c < 0) and then (at t = — f c > 0) percolation 
through black regions disappears. In the range 
t — t c there is percolation both through 

white and black regions. If t c > 0, first per- 
colation through black regions disappears, and 
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only then percolation through white regions 
sets in. Ariel in the range — t c < t < t c per- 
colation is impossible. 

3. In plane problems, percolation through 
black regions excludes percolation through white 
regions, and vice versa. Indeed, if we can sail 
across a mountain ridge from west to east, it 
means that we cannot cross it by land from 
north to south. This excludes the case t c < 0. 
On the other hand, the absence of percolation 
through white regions necessarily means the 
presence of southward percolation through black 
regions. (You can make sure of that by analyzing 
pictures like Fig. 40.) Therefore, the case t e > 0 
is also excluded. This leaves us with the single 
possibility: t c = 0. This is the sought result: 
percolation level stands at zero. 

Formula (2) makes it possible to calculate the 
critical fraction x c of area. As follows from the 
normalization condition for distribution func- 
tion (see formula (1) of Chapter 2), at t c = oo 
the right-hand side of formula (2) equals unity. 
The symmetric nature of function / (y) implies 
that at t c = 0 the fraction of area x c = 0.5. 

In three-dimensional cases, eastward percola- 
tion through white regions does not preclude 
southward percolation through black regions 
because percolation channels can be easily un- 
coupled. (Remember highway bypasses arranged 
at different levels.) For this reason, in the 3D 
case, t c < 0 and, by virtue of formula (2), x c < 
< 0.5. Monte Carlo calculations by the above- 
described procedure for Gaussian random func- 
tions demonstrated that the 3D critical fraction 
x c = 0.16 ± 0.01. 
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The quantity x c can be approximately evalu- 
ated by the method we used to evaluate per- 
colation thresholds of the site problem in Chap- 
ter 6. The method consists in surrounding each 
site by a sphere (or circle in the plane case) of 
radius equal to the distance to the nearest neigh- 
bor. A sphere around a white site is regarded as 
white, and a sphere around a black site is regar- 
ded as black. It was found that percolation 
through white spheres in contact appears when 
the fraction of space filled with these spheres is 
approximately the same for all lattices. It is 
natural to assume that this fraction must not be 
very different from the value x c that we en- 
counter in the percolation level problem. Accord- 
ing to formulas (3) and (4) of Chapter 6, in the 
2D case the fraction of area occupied by white 
circles is approximately 0.5, and in the 3D case 
the fraction of volume occupied by white spheres 
is approximately 0.16. The two estimates thus 
coincide strikingly well with the solution of the 
percolation level problem. We can expect that 
x c = 0.16 is a fairly good estimate for non- 
Gaussian random functions. 

Impurity Compensation in Semiconductors 

Percolation level plays a veryamportantYole in 
the theory of extrinsic semiconductors. Assume 
that a semiconductor was doped by both donor 
and acceptor impurities in identical amounts. 
Donor impurity atoms have an extra electron 
on the outer electron shell, while, contrary to 
this, acceptor impurity atoms lack one electron. 
As a result, donors eagerly donate their extra 
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I! electrons to acceptors and thereby become posi- 

I tively charged. Acceptors accept electrons and 

become negative. (This phenomenon is called 
impurity compensation.) Impurities in semicon- 
ductors being distributed chaotically, a random 
system of positively and negatively located char- 
ges is formed. Each charge produces an electric 
potentm 1 ±e/er, where e is the absolute magnitude 
of the electron charge, e is the dielectric constant 
and r is the distance from the charge. The sign of 
the potential depends on the sign of the charge. 

the potential of any point in space is the 
sum of the potentials produced by all impurity 
atoms. The impurities being distributed ran- 
domly through the semiconductor, their poten- 
tial is also a random function. 

, ^ concen ^ ra tio n of donors is slightly higher 
than that of acceptors, some electrons remain on 
the donors. If the donor-electron bonding enerev 
is relatively low, the thermal motion of atoms 
easily ionizes the donors. In principle, the liber- 
ated electrons can participate in current transfer, 
but this transfer is inhibited by an electric poten- 
tial produced by charged impurity atoms within 
the semiconductor. The product of electric poten- 
tial by the electron charge is the potential energy 
of interaction between an electron and the electric 
held of impurities. Potential energy essentially 
changes the character of motion of electrons. 

Motion of a Particle 
with Nonzero Potential Energy 

The total energy E of a particle consists of its 
kinetic energy mv 2 / 2, where m is the mass of the 
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particle and u is its velocity, and potential energy 
H V (r) which is a function of coordinates of the 
point at which the particle is located: 



where r is a vector from the origin of coordinates 
to the location of the particle. 

The trajectory of a particle is described by a 
function r ( t ), where t is time. The fundamental 
principle of mechanics is that the total energy E 
of a particle does not change in the course of 
motion. This means that a change in coordinates 
of the particle entails a change in its velocity 
so as to compensate for the change in the poten- 
tial energy V (r). The law of energy conservation 
and the fact that the kinetic energy mv 2 l 2 must be 



a positive quantity impose important restrictions 
on the motion of the particle. 

Let us assume for the sake of simplicity that V 
is a function of only one coordinate, namely, X, 
and that velocity is also directed only along the 
X-axis. Let V (X) have the form shown in Fig. 41. 
The type of motion of the particle is determined 
by its total energy. If the total energy is E x , 
the kinetic energy is positive only in the range 
X x < X < X'. Points X x and X' are the bounds 
of motion. Velocity vanishes at these points, hence 
the term turning points . A particle with energy 
E x is trapped into the potential well between the 
turning points and cannot leave this well. Its 
range of motion is bounded on both sides. A par- 
ticle with energy E 2 cannot penetrate deeper 
on the left than point X 2 , but its range of motion 
is not bounded on the right. 
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The motion of a particle can thus be bounded or 
unbounded depending on the relative values of 
the total and potential energies. 

Note that quantum mechanics allows a particle 
to penetrate into regions with a negative kinetic 




Fig. 41. Potential energy V of an electron as a function 
of the X-coordinate (the solid curve). 

energy. However, if V (r) is a sufficiently smooth 
function (and in semiconductors with high con- 
centration of compensated impurities it is always 
smooth), this penetration plays only a minor 
role. 

Motion of an Electron 
in the Field of Impurities 

The potential energy of an electron, stemming 
from the potentials of randomly distributed im- 
purity atoms, is a random function of coordinates. 
If the total energy E of the electron is low, it 
can move only in a bounded region of space 
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everywhere surrounded by turning points. If the 
energy is sufficiently high, the electron is allowed 
to move throughout the whole space. This is the 
necessary condition for participating in the 
transfer of electric charge (in electric current). 

The reader must understand the difference 
between the one-dimensional problem discussed 
in the preceding section and two-dimensional or 
three-dimensional problems. If a particle can 
move only along a straight line, then the motion 
can be unbounded in both directions only if the 
energy of the particle is greater than all the values 
assumed by the potential energy along this straight 
line. Hence, a stringent condition E >> V max 
must be satisfied (V max is the maximum value of 
potential energy). 

No such condition is required in the two- and 
three-dimensional cases. A particle can bypass 
the regions in which its motion is forbidden. It is 
only necessary that the regions where the motion 
is allowed form a system of “lakes and canals” 
through which the particle could travel to infinity 
in an infinite system. 

Now it is obvious that the problem that emerges 
is one of finding percolation level. Let V (r) be 
a random function describing the electron poten- 
tial energy. We fix the total energy E and refer 
to spatial regions where E > V (r) (the positive 
kinetic energy) as white, and to those where 
E < V (r) as black. We want to find percolation 
level, that is, the critical value E c at which 
percolation through white regions sets in. 

Only those electrons whose energy exceeds E c 
are free and take part in current transfer (some- 
times this energy is called the mobility edge). 
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At low temperatures, E c is substantially greater 
than the energy of thermal atomic motion kT. 
Consequently, the probability for an electron 
to reach the level of E e is low. Correspondingly, 
the concentration of electrons capable of charge 
transfer is small. As temperature increases, this 
concentration grows steeply, thereby increasing 
the electric conductivity of the system. The per- 
colation threshold E c thus determines the tem- 
perature dependence of the electric conductivity 
of impurity-compensated semiconductors. 

At very low temperatures, electrons accumulate 
in the deepest potential wells and cannot take 
part in charge transfer. Consequently, at low 
temperatures, impurity-compensated semicon- 
ductors turn into insulators. 

It was mentioned in Chapter 8 that at a suf- 
ficiently high concentration of impurities of 
one sort (e.g. donors), the semiconductor’s con- 
duction becomes metallic, only marginally de- 
pending on temperature down to absolute zero 
(the Mott transition). Impurity compensation 
(e.g. by adding acceptors) enhances the random 
potential energy and switches off the metallic 
conduction. A detailed theory of this phenomenon 
has been developed on the basis of percolation 
level concepts. 



(Critical Behavior 
of Various Quantities 
Near Percolation Threshold. 
Infinite Cluster Geometry. 

$ , : 

l 

In this part we discuss those aspects of percolation 
theory that appear as the most interesting from 
physics’ standpoint: the behavior of various 
quantities in the immediate neighborhood of 
percolation threshold. It was said in the previous 
parts that such physical quantities as the spon- 
taneous magnetization of a doped ferromagnetic 
or the electric conductivity of a network with 
blocked sites vanish at the threshold point. In 
the present part we discuss the laws that describe 
their behavior in the vicinity of percolation 
threshold. The derivation of these laws requires 
that the geometrical properties of infinite clusters 
be understood. 



Chapter 11** 

The Bethe Lattice 

It was shown in Chapter 5 that exact solutions 
can be found in some plane percolation threshold 
problems. However, it was never said that it 
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was also possible to find the function P (x), 
that is, the probability for a site to belong to 
an infinite cluster. No exact expressions for this 
function (or for the electric conductivity of a 
network) are known at present either for plane 
or for three-dimensional problems. The only 
exception is the Bethe lattice which, as we shall 
show below, must be classified as a lattice in the 
space of infinitely high dimensionality. In what 
follows we will give the formulation and solution 
of the site problem on the Bethe lattice. 

Rumors 

“What nonsensical stories they do spread about 
the town. What are things coming to when you 
can hardly turn round before there is some 
scandal going about you, and not a word of sense 
in it either...”* — in these words Nikolai Gogol 
described how a preposterous gossip spread by 
two ladies ruined Ghichikov’s promising fraud. 
“This enterprise they contrived to carry out in 
just a trifle over half an hour. The town was posi- 
tively stirred up; everything was in a ferment; 
and if there were but anybody that could make 
out anything!”** 

And indeed, rumors do spread unbelievably 
fast. But this speed ceases to be surprising if 
we analyze the mathematical model suggested 
below. 

* Dead Souls , a poem by N. Gogol. Ghatto & Windus, 
London, 1922; translated from the Russian by Constance 
Garnett. 

** Dead Souls , a poem by N. Gogol. The Heritage Press, 
New York, 1944; translated from the Russian by B. G. Gu- 
erney. 
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Fig. 42. The Bethe lattice for q = 3 The open circles 
represent people of category 0, and the black circles re- 



The model is shown in Fig. 42a. Assume that 
a “lady agreeable in every respect”, denoted m the 
figure 'by circle A, passed the news to three of her 

acquaintances, B u B ir # 3 . Fach tlies * r . e ® 
passed the message on to three of her (or his) 
rumor-mongering acquaintances, so that this 
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I “second-hand” information was received by nine 
people marked by circles C. Each of these nine 
people again transmitted information to their 
three contacts, making the news known to another 
27 people. It is readily found that the “tenth- 
hand” information reaches 3 10 = 59 049 people! 
Assuming that it takes each gossip twenty min- 
utes to pass the information to his (her) three 
listeners, we find that the whole chain takes 
200 min = 3 h 20 min. 

Of course, this model strongly simplifies the real 
process. It assumes that all people have an iden- 
tical number of acquaintances. Furthermore, it 
assumes that each person receives information 
from only one informant. This means that only 
one line enters each circle (see Fig. 42a). By vir- 
tue of this property, the model resembles a tree 
that branches out infinitely in all directions. 
Each circle can be regarded as a trunk of its 
own tree, and the trees that grew from, for in- 
stance, circles B x , B 2 , B 3 have no common circles. 
The same can be said about the trees stemming 
from circles C , and so on. 

In scientific literature, this model is called 
just that: a tree . It is also called the Bethe lattice , 
after the famous physicist Hans Bethe. The 
circles shown in Fig. 42a are the sites of this 
lattice. The number of lines originating from 
each site of the Bethe lattice can be arbitrary 
(but identical for all sites). Let us denote this 
number by q . Figure 42a shows the lattice for 
q = 3. 

Now recall that most people have their own 
point of view and do not participate in spreading 
rumors. Let us divide all people into two cate- 
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gories: category 0 shown by the open circles 
consists of people who transmit the received in- 
formation to the next trio (Fig. 42h). The se circles 
give rise to q arrows. Category T shown by the 
black circles consists of people who do not par- 
ticipate in spreading the rumors. No arrows 
originate from the black circles (see Fig. 42b). 

The introduction of the black circles strong y 
affects the propagation of a rumor. Consider the 
configuration shown in Fig. 42 b. Among 
three circles B u B v B s , only one is white and 
transmits the rumor further on. Circles C x , L t , 

C 3 would he glad to gossip but received nothing 

from 5 X and B 3 . Among C t , C t , C 6 > i°?{L 

to 0 Therefore, instead of nine recipients of th 
second-hand information, only three were reached 
by the rumor, and only C 6 will transmit it along 

^Imagine that the system under discussion is 
not bounded and thus has an infinite number of 
circles. Then the following question can be 
nosed: Will a rumor originating from point A 
die out after a finite number of transfers, or 
will it stray to an infinite distance from A and 
become known to an infinite number of persons 
in an infinite system? We see from Fig. 426 that 

this depends on the relative number of open aud 
black circles and on configurations that arise 

ar °In n fact, S1 we speak of a site problem of per- 
colation theory, only formulated on the Beth 
lattice. Let the fraction of people belong g 
be x. This means that a'person selected at random 
will be found to belong to category 0 with 
probability x and to category T with probabili y 
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1 — x. The question we have to answer is as 
follows: What is the probability P (x) for a rumor 
transmitted to a randomly selected person to 
become known to an infinite number of people? 
Obviously, this probability is zero if x is small, 
but becomes greater than zero beginning with a 
certain critical value x = x c . 

Solution of the Site Problem 
on the Bethe Lattice 

It will be convenient to introduce, instead of 
P ( x ), the probability for a rumor transmitted 
to a randomly chosen person not to become known 
to an infinite number of persons. Let us denote 
this probability by Q (x). Obviously, 

Q(x) = i-P (x) (1) 

because these two events form a complete system. 

An algebraic equation can be written for 
Q {x). The following line of reasoning must be 
used. The spreading of a rumor can be stopped 
for two incompatible reasons. The first is that the 
person chosen at random belongs to category T. 
The second reason is that although this person 
belongs to category 0 and passes the rumor to q 
people, all channels leading from these people 
will be cut at different stages. The probability Q 
is thus the sum of the probabilities of two in- 
compatible events. The probability for a random- 
ly selected person to belong to category T is 
1 — x. The probability that the person belongs 
to category 0 but that the rumor propagation 
will be interrupted at later stages is denoted, for 
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the time being, by W'. Then 
Q = 1 — X + W' 

Now let us take a look at W'. The outcome 
that it represents is the result of two simulta- 
neously realized events: (i) the person chosen 
randomly belongs to 0 (the probability of this 
event is x ), (ii) all q channels leading from the 
acquaintances of the person chosen at random 
interrupt at some step. Obviously, these two 
events are independent. Hence, the probability 
W ' equals the product of the probabilities: W = 
= xW (#), so that 

Q = 1 — X + X W ( x ) ( 2 ) 

where W (x) is the probability for all q channels 
to be interrupted at some step (of course, this 
can happen at different steps in different chan- 
nels). 

Let us consider one of the q channels that ori- 
ginate from one of the acquaintances of the 
person chosen at random. The event consisting 
in this channel being interrupted at some step 
is equivalent to the statement that the rumor 
communicated to this acquaintance does not 
reach an infinitely large number of people. By 
definition, the] probability of this event is Q (#). 

It is very important for what follows that the 
trees rooted at q acquaintances of a randomly 
selected person have no common circles. Hence, 
a certain configuration of open and black circles 
found in one tree does not affect at all the pro- 
bability of any possible configuration of circles 
in other trees.*" (Obviously, this last statement 
would be incorrect if the trees had common cir- 
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cles.) Consequently, the events consisting in the 
interruption of the rumor in one or another chan- 
nel are independent. 

The probability for all q channels to be inter- 
rupted is therefore equal to the product of the 
probabilities for each of the q channels to be 
interrupted: 

W (x) = [<? (x)l 3 (8) 

Substituting formula (3) into (2), we obtain 
an equation for Q ( x ): 

Q (x) = 1 — x + x (x)] q (4) 

Note that the decisive factor that enables us 
to reduce the problem to algebraic equation (4) 
is the independence of channels. This property 
is an exclusive property of the Bethe lattice, 
so that the method we use is not successful when 
applied to usual lattices, although it is often 
employed to obtain an approximate solution. 

Let us analyze equation (4). It is meaningful 
for all x within the interval 0 ^ x ^ 1. Re- 
writing it in terms of P (x) — 1 — Q (x), we find 

H — P(x)] 9 x-i-P(x)—x = 0 (5) 

One of the solutions of (5) is P (x) = 0 for all x , 
although for q > 1 equation (5) is nonlinear and 
has other solutions. Thus, P (1) = 1 is also a 
solution for x = 1, and this second solution is 
physically meaningful because the probability P 
must be unity and not zero if all circles are 
open. 

For q = 2 the solutions of equation (5) are 
easily found. There are two solutions in this 
case: P (x) — 0 and P (x) = 2 — 1 lx. If x > 
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> 1/2, the second solution has a physical mean- 
ing. If x < 1/2, the solution becomes negative 
and thus meaningless. 

Therefore, for q = 2 we have the following 
solution: 

( 0 for0<z<l/2 

P{x)= | 2 — 1/* for 1/2<*<1 

In this case the percolation threshold x 0 equals 

1 / 2 . , W1 
A similar solution exists for all q > 1, although 
the percolation threshold x c is a function of q. 
In the general case it is possible to find x c and 
the form of P (x) for x close to x c , assuming 
from the start that P (x) < 1 (this is always 
true in the neighborhood of percolation thre- 
shold). The term (1 — P) q in equation (5) can be 
expanded via the binomial theorem: 

( l_P). = l- g p + ’-fc! (7) 

Each subsequent term here is much less than the 
preceding term because P 1 • Therefore, we 
substitute (7) into equation (5), assuming that 
only the three terms written above are significant 
on the right-hand side of (7). This gives 

- ^ xP z = qPx — P 

z 

Assuming P 0, we divide both sides of this 
equality by P and obtain 



p (s— l/g)-2 

x(q-i) 



( 8 ) 
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Formula (8) coincides with (6) for q = 2. This 
solution vanishes at x = l/q, which implies that 
x c — l/q. Solution (8) is meaningful if q > 1 , 
x > l/q , and only if x is very close to l/q. There- 
fore, we can set x = l/q in the denominator of 
(8). Finally, we obtain 

P(x) = (x ~!^}- 2 g - (9) 

Formula (9) describes the probability P (x) in 
the vicinity of percolation threshold. 

Discussion 

If q = 1, the function P (x) = 0 for all x within 
the range 0 ^ x < 1 , and equation (5) has no 
other solutions within this range. If q — 1 and 
x — 1, equation (5) holds for any value of P. 

If q = 1, the Bethe lattice transforms into a 
linear chain of lattice sites. And no matter how 
small number of black sites interrupts percola- 
tion through white in such a chain. It is therefore 
natural that whatever the value of x within the 
range 0 ^ x < 1, the interruption of the rumor 
propagation is unavoidable, that is, P (x) — 0. 
In a linear chain x c = 1. 

We have shown in the preceding section that 
in the general case x c = l/q. This result could 
be predicted in advance. Each person in the 
problem under discussion transmits the rumor 
to q of her (his) acquaintances. The mean number 
of people of category 0 among these acquaintances 
equals qx. Hence, each transfer of information on 
the average creates qx sources instead of one. 
The quantity qx is therefore its branching ratio. 
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In order for the process to continue, the branching 
ratio must be greater than unity. Hence, the 
critical concentration x c is obtained from the 
condition qx c = 1, that is, x c = l/q. 

Recall that the condition necessary to sustain 
the uranium fission chain reaction is written 
in the same form. In fact, the process of rumor 
spreading is also a chain reaction and is des- 
cribed in the same terms as a nuclear explosion. 

It is of interest to compare the obtained value 
of x c with the results derived for lattice prob- 
lems in spaces with a greater dimensionality. 
Percolation threshold of the site problem was 
approximately computed for the so-called hypei- 
lattices. These are lattices of the same type as 
the square and simple cubic ones but in a space 
with a greater number of dimensions. The coor- 
dination number z (the number of nearest neigh- 
bors) is found in such lattices from the formula 
z = 2d, where d is the dimensionality of space 
(it gives z = 4 for d = 2, and z = 6 for d — • )• 
The percolation threshold £ s was computed tor 
d = 4, 5, 6. It was found that the formula 

_/ 1 ,6J\_ L_ (10) 

x s (J 2 / z — 1 

provides a good fit to the obtained results. 

The second term in parentheses can be neglect- 
ed if d is sufficiently large, so that * 3 = 
= l/( z — 1). But recall that z — 1 = q for the 
Bethe lattice. Indeed, only one bond enters 
each site of this lattice, and q bonds originate 

We see that percolation threshold of the Bethe 
lattice (x e = 1/?) is the same as that of a hyper- 
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lattice with a very large dimensionality. Hence, 
the Bethe lattice as if represents an infinite- 
dimensional space. 

The Bethe lattice is the only system for whicli 
it proved possible to find the exact form of P (x) 
at percolation threshold. It was found (see (9)) 
that in this case P (x) tends to zero linearly as 
x -> x c : P OC (x — x c ). We shall see later that 
this is a specific property of the Bethe lattice as 
well as of all lattices of high dimensionality. 

Exercise 

1. Analyze the bond problem on the Bethe lat- 
tice. Consider all sites to be identical, and the 
bonds to be intact or broken. Let the fraction 
of intact bonds be x. Find the function P (x) 
defined as in the text. 



Chapter 12 

Structure of Infinite Clusters 

The Shklovskii-de Gennes Model 

Now we shall consider, for the sake of definiteness, 
the site problem and assume that the concentra- 
tion of nonblocked sites is slightly above the 
threshold value, so that an infinite cluster exists. 
It consists of infinite strings (chains) of inter- 
connected sites. If all connected sites of an 
infinite cluster are joined by straight segments, 
we obtain a set of intersecting broken lines (see 
Fig. 15 showing one such line). 
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By definition, the structure of an infinite 
cluster is its geometry on a scale much greater 
than the lattice period. On this scale, the sharp 
knees occurring at individual lattice sites are 
not resolved to the eye, and strings appear as 
smoothly bent curves. 

Figure 43 shows a small fragment of an infinite 
cluster. The cluster does not end at points A 




Fig. 43. Fragment of an infinite cluster with dead-ends 

and B but stretches leftwards and rightwards to 
infinity. Now let us introduce the following clas- 
sification of points and lines in an infinite cluster: 
the elements of an infinite cluster belong to either 
its backbone or dead-ends. 

A point is said to belong to the backbone of an 
infinite cluster if there are at least two paths that 
emerge from this point in opposite directions 
and lead to infinity. Point C in Fig. 43 is one of 
such points. Going to the left or to the right of 
this point, we can move infinitely far from it. 
If only one path leaving a point leads to infinity, 
this point belongs to a dead-end. For instance, 
only the motion upwards from point D in Fig. 43 
leads to infinity. The downward motion ends at 
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a cul-de-sac. Point D is thus said to lie on a dead- 
end. 

Let us imagine that all dead-ends are deleted; 
now we can try to discern the structure of the 
backbone of an infinite cluster. The simplest 
backbone model was suggested independently by 
the Soviet physicist B.I. Shklovskii and the 




Fig. 44. Backbone of an infinite cluster. 

French physicist P.G. de Gennes. In a plane 
problem, this model resembles a very large, old, 
and fairly worn-out fishnet. The net lost its 
regular periodicity, its cords have slack, some 
nots are broken, while others slipped off to wrong 
locations, but nevertheless it’s still a net on 
the average” (Fig. 44). 

The characteristic linear size R of a cell ol 
this net is called the correlation radius of an in- 
finite cluster. It grows drastically as we approach 
percolation threshold: 



where l is a length that in order of magnitude 
equals the lattice period, and v is a positive 
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quantity called the exponent of correlation radius. 
The network thus becomes less and less dense as 
percolation threshold is approached. 

A correlation radius approaching infinity is 
the general property of all critical phenomena. 




Fie. 45. Correlation radius as a function of ^- The graph 
shows the width 5 of the critical region for the L X L 
square (see the next section). 

The power-law behavior (1) is not a rigorously 
proved property; nevertheless, it forms the basis 
for the modern theories of critical phenomena 
and appears to be well supported by experimental 

da A correlation radius remains meaningful at 
x<x c , that is, below the threshold. In this 
region it describes the maximum size of finite 
clusters. If x ->■ x c on the side of smaller values 
(x C xA a correlation radius also tends to 
infinity as in (1). This means that when percola- 
tion threshold is approached from below, finite 
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clusters grow infinitely and merge into an infinite 
cluster at x = x c . The dependence of R on x 
thus has a shape schematically traced in Fig. 45. 

For three-dimensional problems the Shklovskii- 
de Gennes model is formulated similarly. It 
resembles a badly damaged wire skeleton of a 
three-dimensional lattice whose mean cell size 
is given by formula (1). But it must be borne 
in mind that the numerical values of the expo- 
nents of correlation radius are different for plane 
and 3D problems. 

Let us consider now the corollaries derived from 
the network-like structure of infinite clusters. 

Role of the System’s Size 

We emphasized in Chapters 1, 2, and 3 that the 
concept of percolation threshold is truly meaning- 
ful only in an infinite system. In a finite system, 
percolation threshold varies among samples, that 
is, it is a random variable. However, the values 
assumed by this random variable fall, with very 
high probability, into a certain range of width 
6 (jfP) called the critical region. As the number of 
sites in the system grows, the width of this region 
undergoes a power-law decrease (see formula (8) 
of Chapter 1), so that percolation threshold 
acquires a clear-cut sense as^-> oo and con- 
verts from a random variable into a certain 



quantity. 

This information was practically quoted, with 
no attempt to derive it, in the opening chapters 
of the book. The concept of correlation radius 
makes it possible to understand these concepts 
and to derive formula (8) of Chapter 1. 
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For definiteness, let us discuss the wire-mesh 
experiment on a grid L X Lin size, the schema- 
tic of this experiment being shown in Fig. 1. 



Assume that numerous experiments using various 
random sequences of blocked sites have been 




Fig. 46. Square superposed on an infinite network, with 
R <C L and x > x c . 



conducted, and the result is a set of percolation 
threshold values. Recall that the configurations 
of blocked sites are very dissimilar in different 
experimental runs. 

It is convenient to pursue the following line of 
reasoning. Imagine an infinite network (a wire 
mesh), with a prescribed bed fraction x of non- 
blocked sites. Imagine now that an L X L square 
is superposed on the mesh at different places, and 
percolation is analyzed from the left- to the 
right-hand side of this square through the non- 
blocked sites that fell within the square (Fig. 46). 
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Placing the square over the different parts of 
the infinite wire mesh, it is possible to sample 
one by one the results of successive experiments 
with a finite wire mesh. 

Percolation sets in in an infinite network exact- 
ly at x = ar c , but as we shall see below, this does 
not mean at all that percolation is always realized 
in an L X L square when x > x c . 

If x > x Ci an infinite cluster exists in the 
infinite system. Let us picture its backbone as 
a fishnet (Fig. 44). The parameter extremely 
important hereafter is the ratio of correlation 
radius R to the square size L. First we assume 
that L is much greater than R. Then (see Fig. 46) 
the square covers a large number of unit cells 
of the infinite cluster netting that ensures per- 
colation between the opposite sides of the square. 
These cells may be different in size, and the net- 
ting of the infinite cluster may contain large 
holes, but the probability for the cluster to con- 
tain a square-size hole is negligibly small if on 
the average the square contains very many cells. 
The following conclusion is therefore drawn: 



If x > £ c , percolation threshold of a square 
cannot lie' within the range of x satisfying 
the strong inequality L > R (z). This region ' ' 

must lie above the threshold. 

According to formula (1), a correlation radius 
grows infinitely when x tends to x c , so that at 
some value of x the radius inevitably becomes 
comparable with L. Now we cannot say anything 
definite about percolation across this square. 
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Everything depends on the specificities of the 
configuration of blocked sites within the square. 

Now let x < x c , with the correlation radius 
much smaller than L. At x <C x c the correlation 
radius is the maximum length of finite clusters. 

If R -C L, no cluster exists long enough to con- 
nect the opposite sides of the square. Hence, 
another definite conclusion is made: 

If x < Xf,, percolation threshold of a square 
also cannot lie within the range of x satisfy- (I1) 
ing the strong inequality L >R (x). this 
region must lie below the threshold. 

If x < X e but X is very close to x c , the cor- 
relation radius exceeds L. In this case we cannot 
say anything definite about percolation across 
the square. An infinitely large system contains 
finite clusters whose size is greater than L, but 
these clusters have holes of a similar size, so 
that everything depends on a specific configuration 
of blocked sites within the square. 

Now we can evaluate the size of the critical 
region in which the values of percolation thresh- 
old for an L X L square can lie. According to 
the corollaries (I) and (II), this region must obey 
the condition L ^ R. Figure 45 shows that this 
region is the narrower and squeezes up to per- 
colation threshold for an infinite system the 
closer, the larger L is. The width 8 of the region 
obeys the condition R (6) = L. From formula 
(1) we find l! S v = or 

Mir 



(2) 
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Percolation thresholds of L X L squares are 
distributed uniformly within the critical region, 
that is, for \ x — x c | <C 8 (see Fig. 5 where the 
distribution function of percolation thresholds 
is plotted). Nothing singles out the point x — x G 
within this region. Indeed, this is the point at 
which percolation sets in in an infinite system. 
But it is impossible to establish whether per 
eolation is achieved or not if we work with finite- 
size squares. If L < /?, it is impossible to estab- 
lish, by placing the square over different regions 
of an infinite network, whether the network con- 
tains only finite clusters or whether these clusters 
have already merged and formed an infinite clus- 
ter. A study of percolation across a finite-size 
square only allows us to determine the width of the 
critical region. 

In this section we discussed only plane prob- 
lems. In fact, the arguments completely transfer 
to three-dimensional problems. Formula (2) de- 
termines the width of the critical region in 3D 
problems equally well. A slight difference appears 
if the width 5 is written not in terms of the 
system’s size L but of the total number of sites, jfT . 
The point is that JfT = (. L/a ) d , where a is the 
lattice period, and d is the dimensionality of 
space. Consequently, we have by (2) 

( 3 ) 

where C is a numerical coefficient that cannot be 
found by equally simple reasoning. In plane cases 
(d = 2), formula (3) coincides with formula (8) 
of Chapter 1. It was by means of this formula 
that the exponent of correlation radius of the 
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plane problem was established for the first time 
from the dependence 8 (jfT) determined numeri- 
cally on a computer. It was found that v 2 = 1.33. 
(Here and below the subscript 2 indicates that 
we deal with the exponent of a two-dimensional 
system ) The exponent v of three-dimensional 
problems is different: v 3 = 0.8-0.9. (The sub- 
script 3 denotes that the exponent refers to oD 
problems.) 



Electric Conduction 
Near Percolation Threshold 

Let us consider specifically two- or three-dimen- 
sional networks with blocked sites. As we said 
in the opening chapters, such networks are 
electrically conductive at x !> x C7 but at the 
percolation threshold x c their conductivity va- 
nishes. Experimental data and the data obtained 
by computer simulation show that the specific 
electric conductivity of networks tends to zero 
in a power-law fashion, 

a (x) = a 0 (x — x c Y (4) 

where the factor o 0 equals, by the order of ® a S" 
nitude, the specific electric conductivity of the 
network without blocked sites. The quantity t 
is called the critical exponent of electric conduc- 
tion; it constitutes the subject of a very careful 
analysis, mostly by numerical computer simula- 
tion (e.g. one of the more recent computations 
involved a square lattice with 800 X 800 sites). 
It was established that t a = 1.3 in two-dimen- 
sional networks, and t 3 = 1.6-1, 7 in three-di- 
mensional networks. 
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The network model of an infinite cluster makes 
it possible to derive formula (4) and to relate 
its exponent t to the exponent of correlation 
radius. Electric current flows only through an 
infinite cluster, namely, through its backbone. 
No current flows in dead-ends that couple to the 
backbone only at one end. If electric current 
were sufficiently strong, so as to heat the wire to 
a glow, the backbone of an infinite cluster would 
be visually discernable in the dark, as a rate of 
i uminated channels against a dark background, 
bar from threshold the whole network glows more 
or less uniformly, with the distance between glow- 
,? n ^ s ' increasing as we approach the 
threshold; finally, the glow fades out when the 
threshold is reached: the current ceases to flow. 

Let us calculate the specific electric conducti- 
Vi * l , ba ? kb °ne of an infinite cluster. It 
must be borne in mind that we cannot expect to 
ca culate correctly the numerical factors. Our 
calculation will only give how a depends on 
x ~ x This dependence will not be altered if 
we replace an irregular, distorted network bv 
an ideal network with period R. 

First consider a plane case (see Fig. 46) The 
resistivity equals the resistance of a square with 
unit length side. The number of wires intersect- 
ing this square is 1 /R, where R is the spacing 
between wires, given by formula (1). Let us denote 
the resistance of one wire with unit length by p ft . 
The circuit is the parallel connection of all wires 
therefore, the resistivity is 

p = (1/i?) = (5) 
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and the electric conductivity is 

a=p-‘ = p- 1 i?-‘ (6) 

Substituting (1), we obtain 

a = a 2 (x — £ c ) v U) 

where <r 2 = Po'J -1 - , , , 

In the three-dimensional cases we need to 
calculate the resistivity of a wire skeleton that 
models, for instance, the simple cubic lattice 
with period R (lattice type affects only the 
numerical coefficient). The resistivity equals the 
resistance of a cubic cell with unit-length edge. 
The number of wires connected in parallel pass- 
ing through a face of such a cube equals \IR 2 . 
The resistivity therefore equals 

p= iiw r=pofli (8) 

and the electric conductivity is 

a = po 1 ^ -2 = a 3 (^ — x c) 2V ( 9 ) 

where a s = po^ -2 - , , 

In order to avoid confusion, note that tne 
electric conductivity cr has different dimension- 
alities in the two- and three-dimensional cases: 
it is measured in ohm”*' in the 2D case, and in 
ohm” 1 -cm” 1 in the 3D case. 

The coefficients o 2 and a 3 correspond, by the 
order of magnitude, to the electric conductivities 
of two- and three-dimensional networks with no 
blocked sites. Indeed, as we see from formulas 
6) and (9), the electric conductivity a (x) reduces 
to a 2 or a 3 , respectively, if R = l, that is, if the 
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network of an infinite cluster coincides with the 
initial network on which the problem has been 
formulated. The factor <r 0 in formula (4) is thus 
equal to a 2 in the 2D case, and to o 3 in the 3D 
case. 

A comparison of formulas (6) and (7) with 
formula (4) shows that t = v in the 2D case, and 
t — 2v in the 3D case. Using v 2 — 1.3 and 




Fig. 47. Backbone of an infinite cluster, with “sinuosity” 
taken into account. 

v 3 = 0.8-0. 9, we find t 2 = 1.3 and t 3 = 1.6-1. 8, 
quite close to the data given above. This agree- 
ment strongly supports the Shklovskii-de Gennes 
model. 

The backbone model of an infinite cluster can 
be generalized as follows. Imagine that the wires 
forming the backbone are very sinuous (Fig. 47). 
The distance between their intersection points 
is, as before, R (. x ) and is given by formula (1). 
However, if the wire segment between two inter- 
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section points is straightened out, its length 
will prove considerably greater than R. We 
denote this length by X and write it in the form 



# = 



l 

(x — x c )^ 



(x > X c ) 



where £ > v. As x -> x G1 the ratio XI R = 
= (x — # c )-(s-v) tends to infinity. 

The introduction of the length X is justified 
only if X > R, that is, £ > v. However, it has 
been proved that £ = 1. Hence, the suggested 
generalization of the Shklovskii-de Gennes model 
(in fact, this generalization has been carried by 
the authors of the model) is meaningful only if 
v < 1. Since v = 1.3 in the 2D case, no “sinuosi- 
ty” is expected there. But in the 3D case v < 1, 
and there is a good reason to believe that “sinuosi- 
ty” indeed characterizes the backbone of an infi- 
nite cluster. 



Exercise 

1. Express the critical exponent t of electric 
conductivity in terms of the exponents £ and v 
in the 3D case, making use of the generalized 



model. 



Function P(x) Near Percolation Threshold. 

Role Played by Dead-Ends 

The function P (x), which is the fraction of sites 
belonging to an infinite cluster, vanishes at x = 
— x , together with electric conduction. An 
analysis demonstrated that in the immediate 
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vicinity of the threshold this function varies as 

P (x) = D (x - x c )t (10) 

where D is a numerical coefficient of the order 
of unity, and (} is another critical exponent. It 
has been found that p 2 = 0.14 (2D problems), 
and p 3 = 0.4 (3D problems). These results were 

I mostly obtained in computer simulation. 

All sites of an infinite cluster, both those be- 
longing to the backbone and to dead-ends, con- 
tribute to P (x). The infinite cluster model makes 
it possible to find which of these subsets has a 
greater number of sites. First assume that dead- 
ends are completely absent and calculate the 
I contribution of the backbone of the infinite clus- 
ter to P (x). 

The number of sites belonging to the backbone 
per cell of an infinite cluster in the 2D case is of 
the order of R/a, where a is the lattice period 
(as in the preceding section, this is also an 
evaluation that is not meant to find numerical 
coefficients). The cell area is about R 2 , and 
hence, the total number of all sites in a cell is 
about R 2 /a 2 . Therefore, the fraction of sites 
belonging to the backbone of an infinite cluster is 

^bb ( x ) ~ -Jf ~ (* — X c ) v (H) 

Here the symbol ~ stands for equality to within 
an order of magnitude (neglecting numerical 
coefficients of the order of unity). 

[ In the 3D case the number of sites of the back- 

| bone per each cell of an infinite cluster is also 
R/a 1 but the total number of sites per cell is 
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about (R/a) 2 . Consequently, in the 3D case 

(x) 2 ~ (a: ' - * c)2v (12) 

Comparing formulas (11) and (12) "frith for- 
mulas (7) and (9), we obtain that the fraction 
of sites belonging to the backbone of an infinite 
cluster coincides, to within an order of magni- 
tude, with the function o ( x)la 0 = (x — x 0 ). 
Comparing (11) and (12) with formula (10), we 

notice that 



_ , x _ x 
P (x) V c ' 

in the 2D case, and 

P bb(*) _ ( r — r \ 2v >-P* 

__T - (x ar 0 ) 



in the 3D case. - 

Recall that v 2 = 1.3 and v 3 « 0.9. Hence, 
v, - 6, = 1.2 and 2v s - p„ « 1.4. Therefore, 
the ratio P bb {x)IP (x) rapidly tends to zero 
as x X(S both in the 2D and 3D cases. This 
means that the sites forming the backbone of an 
infinite cluster comprise an infinitesimal fraction 
of the total number of sites belonging to the 
infinite cluster. The “ mass ” of an infinite cluster 
predominantly resides in its dead-ends and is com- 
pletely useless from the standpoint of electric 
conduction. Consequently, in the vicinity of 
percolation threshold a (x)/a 0 « P (x) (see 
Fig. 10). However, the spontaneous magnetiza- 
tion of a doped ferromagnetic in the neighbor- 
hood of percolation threshold is determined 
precisely by these dead-ends (see Chapter 3). 
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Universality of Critical Exponents 

We have encountered three critical exponents 
that describe the behavior of various quantities 
in the neighborhood of percolation threshold, 
i.e. v, t, and p. This behavior is also referred to 
as critical because the functions R (x), a (x) and 
P ( x ) have singularities at x = x c . For instance, 
R (x) tends to infinity, while P (x) has a dis- 
continuous first derivative (at x c ). The derivative 
is zero on the left of x c and is infinitely large on 
its right. In the case of g(x), it is a second deriv- 
ative that is discontinuous. Many other quan- 
tities with critical behavior are known in per- 
colation theory, and correspondingly, many other 
critical exponents. 

For each of the above-discussed critical expo- 
nents we gave two values: one for 2D and another 
for 3D problems. However, the list of relevant 
2D and 3D problems would be extremely long. 
For instance, we know of three-dimensional site-, 
bond-, and sphere problems, the problem of per- 
colation level in a random potential, and many 
others. What, then, were the exponents we dis- 
cussed so far? Now we come to what seems to be 
the most exciting feature of percolation theory. 
According to modern concepts, the critical expo- 
nents are identical for all problems in the space of 
a given dimensionality. (The only exception among 
the problems mentioned in this book is the prob- 
lem of directed percolation.) The statement of 
universality of critical exponents is rather a 
convention than a rigorously proved proposition. 
Nevertheless, numerous tests of this statement, 
carried out on computers, failed to refute it. 
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What are the physical reasons for the univer- 
sality of exponents? Presumably, the exponents 
are determined by the structure of clusters in the 
vicinity of percolation threshold. The properties 
that play the principal role here are the geometri- 
cal properties of clusters that are felt at large 
distances (of the order of correlation radius). 
In the neighborhood of the threshold these distan- 
ces become much greater than the lattice period (in 
lattice problems) or the radius of spheres (in sphere 
problems). Consequently, the geometry of clusters 
is independent of the type of lattice on which 
the problem is formulated. The structure of large 
clusters will remain unaltered even if a problem 
is given not on a lattice but on sites randomly 
distributed in space. But, of course, the dimen- 
sionality of space strongly affects the geometry 
of clusters because, for instance, it is much easier 
to ensure “bypassing” of curves in a three-dimen- 
sional space than in two dimensions. For these 
reasons, critical exponents do not depend on 
problem type, but do depend on space dimen- 
sionality. 

It is interesting to note that critical exponents 
change with increasing dimensionality of space 
only up to dimensionality 6. Beginning with 
d = 6, exponents remain constant, and the expo- 
nent p == 1 as in the Bethe lattice. When d ^ 6, 
the critical exponent problem considerably sim- 
plifies and allows an exact solution. 

Critical exponents thus possess a certain uni- 
versality in contrast to percolation thresholds 
which strongly depend on the type of problem. 
This leads to a simple conclusion. If the results 
of a physical experiment are treated in terms 
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of percolation theory, and the microscopic struc- 
ture of the investigated system is not quite clear, 
the first characteristics to be compared with the 
theory must be critical exponents because they 
are almost invariant. This is the strategy one 
uses in analyzing experimental data on electric 
conduction in heterogeneous materials (see Chap- 
ter 9). 

Percolation theory borrowed the idea of uni- 
versality of critical exponents from the theory of 
second-order phase transitions (e.g. among second- 
order transitions are the transition of a metal 
from the superconducting to the normal state and 
the transition of a ferromagnetic to the nonferro- 
magnetic state, both brought about by increasing 
the temperature of a sample). As in the vicinity 
of percolation threshold, large regions differing 
in properties from one another form close to the 
point of a second-order phase transition. The 
difference lies in that the boundaries of these 
regions are not “frozen” as they are in percolation 
theory, but change with time owing to thermal 
motion. The size of such regions is also called 
the correlation radius and also given by for- 
mula (1). 

The theory of phase transition supplied another 
important idea as well: the scaling hypothesis. 
We shall formulate it in terms of percolation 
theory. Assume that different parts of a net- 
work with some blocked sites have been photo- 
graphed at two values of x , namely, x x and x 2 . 
Both these values are on one side of x c and are 
close to it. Let, for instance, x x — x c > x 2 — 
— x c > 0. At x = x x the correlation radius is 
R x = l (x x — £ c )" v , and at x — x 2 it is R 2 — 
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__ / t x x )~ v . These conditions imply that 

R 1 < J? 2 - The scaling states: if the photographs 
taken at x = x x are magnified by a factor RJRii 
they will not differ, “on the average , from the 
photographs taken at x = x 2 . The photographs 
are assumed to be so coarse-grained that only 
large blocks show up, and individual sites and 
bonds are not resolved. The scaling hypothesis 
thus states that the large-scale geometry of the 
system is transformed on approaching percolation 
threshold in a self-similar manner, with all linear 
dimensions growing in proportion to the corre- 
lation radius. 

Note that the Shklovskii-de Gennes model 
satisfies the scaling hypothesis, but that this 
hypothesis goes far beyond the model. It is neither 
restricted to the backbone of an infinite cluster 
nor does it assume that the network is subdivided 
into its backbone and dead-ends. 

The mathematical formulation of scaling hypo- 
thesis makes it possible to find a relation between 
the critical exponents v, p and a third exponent 
that we chose not to introduce here. Calculations 
show that this relation holds quite well. 

Scaling concepts that were first introduced by 
the Soviet physicists A. Z. Patashinsky and 
V. L. Pokrovsky and by the American physicist 
L. Kadanoff form the basis of the modern 
theory of phase transitions and of percolation 

theory. , f 

Scaling concepts brought about a number or 
new mathematical techniques for calculating 
critical exponents. A dramatic progress has been 
achieved in these methods in recent d_ eca “ e > so 
that certain perfection has been reached. At pre- 
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sent, these methods help to calculate almost all 
important exponents. However, computations 
are so complicated that they could not be dis- 
cussed in this book. 



Chapter 13 

Hopping Electric Conduction 

The preceding chapters gave a detailed picture 
of how percolation theory is applied to calculate 
the electric conductivity of systems consisting 
of randomly coupled elements. Each element 
was either conducting or insulating, but the 
resistances of all conducting elements were as- 
sumed equal. Examples of such systems are net- 
works with broken bonds or blocked sites, mix- 
tures of metallic and dielectric spheres, and so 
forth. 

Let us turn now to a different class of systems 
which are also composed of different elements, 
but the resistances of these elements can take 
on any values from extremely small To enorm- 
ously large. It was found that the 'resistance 
of systems composed of a very large number of 
such elements can also be calculated in terms 
of percolation theory. The theory of hopping 
electric conduction of semiconductors, based on 
percolation ideology, was constructed in 1971. 
In this chapter we make acquaintance with hop- 
ping electric conduction and present its mathe- 
matical description. 
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Mechanism of Hopping Conduction 

Consider an extrinsic semiconductor doped with, 
say, donor impurities. The dopant concentration 
will be considered small in comparison with the 
critical concentration N c at which the Mott 
transition (metal-insulator transition, see Chap- 
ter 8) takes place. In these conditions, outer 
electron shells of neighboring impurity atoms 
overlap very slightly. Each donor exists therefore 
as an isolated hydrogen-like atom whose outer 
electron is at a distance about a% from the nucleus 
and whose bonding energy is about E% (see Chap- 
ter 8). Assume that the temperature of the semi- 
conductor is so low that the thermal energy o 
vibrating atoms is not sufficient to ionize a donor. 
What could be the mechanism of electric con- 
duction in this situation? 

Imagine that some donors lost their outer elec- 
tron Usually this results from impurity com- 
pensation (see Chapter 10). If the semiconductor 
contains both donor and acceptor impurities, 
each acceptor accommodates one electron from 
a donor. If the number of acceptors is smaller 
than that of donors, a fraction of donors will have 
retained an outer electron, while the other frac- 
tion will have lost their electrons and become 
positively charged. 

The mechanism of hopping electric conduction 
consists in the “hops’ 7 of an electron from one 
donor to another which before this hop had no 
outer electron. 

Now we are going to consider the case ot an 
outer donor electron whose potential energy 
depends only slightly on the spatial position 
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of the donor, that is, on the specificities of the 
configuration of charged impurities surrounding 
this donor. This case corresponds to a low con- 
centration of impurities. In conditions typical 
tor hopping conduction, the spread in the ener- 
gies of outer electrons belonging to different 
donors is roughly 0.1 of the bonding energy E%. 

In this situation, the attraction of an electron 
to the donor to which it belongs at the initial 
moment is the main obstacle to a “jump” from 
one donor to another. From the standpoint of 
classical mechanics, an electron can be trans- 
erred from the outer orbit of one donor to the 
outer orbit of another if some work is done 
against attractive forces because it is necessary 
to remove the electron to one half of the distance 
between the donors. After this point the electron 
is attracted to the second donor. If donors are 
spread very thinly, this work is almost equal to 
the work required to remove the electron belong- 
ing to an isolated donor to infinity, that is, equal 
to the bonding energy E%. 

What then is the advantage that hopping elec- 
tric conduction holds over the electric conduc- 
tion by free electrons? The point is that hopping 
electric conduction is a quantum phenomenon. 
Quantum mechanics allows an electron to pass 
from one donor to another without rising to a 
tree state and without borrowing energy from the 
thermal motion of atoms. This is the tunnelling 
transition. Tunnelling transitions obey the law 
of energy conservation. The law imposes the con- 
straint that the electron energy be equal in the 
initial and final states. Consequently, if the 
energies of the electrons at the first and second 
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donors differ owing to the potentials of the sur- 
rounding impurities, a deficit of energy has cer- 
tainly to be borrowed, and a surplus of energy 
has to be dissipated. But this energy is ten 
times less than E%. As a result, hopping electric 
conduction defeats the electric conduction by free 
electrons at very low temperatures. 

Resistor Network 

It must not be overlooked, however, that a tun- 
nelling transition is a very low-probability event 
in conditions we discuss. As was said in Chapter 8, 
the probability of finding an outer electron of 
a donor at a distance r from the donor’s nuc- 
leus decreases with r as exp (—2 r/ag). If two 
donors 1 and 2 are at a distance r 12 ^ from 
each other, the probability of finding the outer 
electron of donor 1 near the nucleus of donor 2 
equals exp (—2 r 12 /ag). This exponential enters 
the probability of tunnelling transition. As fol- 
lows from Chapter 8, at a donor concentration 
much less than the critical concentration N c 
of the metal-insulator transition, the mean dis- 
tance between donors is much greater than ag, 
so that the quantity exp (—2 r 12 /a^) = 

= 1/exp (2r 12 /aB) is, as a rule, very small. 

Nevertheless, tunnelling transitions do take 
place time and again between neighboring donors. 
If an electric field is applied to a semiconductor, 
these transitions will be more frequent in the 
direction against the field (along the applied 
force) than along the field. The result is an elec- 
tric current proportional to the electric field 
strength. This is the essence of the phenomenon 
we call hopping conduction* 
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The model used to calculate the resistivity of 
a semiconductor is the so-called “resistor net- 
work”. This model is formulated not in terms 
of atoms and tunnelling transitions but of con- 
ventional electric circuit resistors. Imagine that 
a resistor is connected between each pair of 
donors. The donors themselves can be pictured 
as tiny metallic balls to which the wires from 
numerous resistors are soldered. The second end 
of each of these resistors is soldered to another 
ball. The result should resemble an irregular 
three-dimensional resistor network. Obviously, 
an intention to assemble an actual network model 
would require that the scale of the system be 
substantially magnified. Indeed, the mean dis- 
tance between donors is of the order of 10 ~ 5 cm. 

The resistors connecting two donors must be 
found by calculating the tunnel current that 
flows between these donors at a given electric 
field. This means solving a quantum-mechanical 
problem that will not be given here. Note only 
that, in view of the above description, the elec- 
tric current produced by tunnelling transitions 
between donors in a given field is the smaller, 
the greater the separation between these donors is. 

Accordingly, the resistance connecting do- 
nors at a distance r from each other can be written 
in the form 

M ( r ) = exp (2r/ag) (1) 

where is identical for all resistances (it can be 
of the order of 1 ohm). The problem is to find 
the resistivity of a system composed of a tre- 
mendously large number of donors (10 16 -10 19 ) 
with ran dom patial distribution. 



Properties of Resistor Network 



The main feature of the model we discuss is that 
the resistances described by formula ^(1) are 
spread over a fantastically large range. The 
mean distance r don between donors is related to 
the donor concentration iV d on by (4/3) itr don /V don = 
= 1, which means that on the average one donor 
is nested in a sphere of radius r don . Typically, 
hopping electric conduction is observed in con- 
ditions in which r don exceeds a B by a factor of 
6 to 12. 

Assume, for instance, that r don = 10a B . Then 
the resistance connected between donors at a 
distance 1.5r don is greater than that between 
donors at a distance r don by a factor 
exp(r don /ag) = * 10 = 2.2-10 4 . 

Donors at a distance r don occur almost as fre- 
quently as those at a-distance 1.5r don . Conse- 
quently, a negligible change in separation be- 
tween donors results in a tremendous variation 



of resistance connected between them. 

In principle, the resistor network model stip- 
ulates that each pair of donors is connected 
through a resistor. However, the resistances 
between remote donors are so large that they 
can undoubtedly be neglected. As a rule, the 
same two donors are connected through a chain 
of resistors connecting nearest-neighbor donors. 
Although the length of this chain is greater than 
the shortest distance between remote donors, the 
chain resistance is much less than one resistance 
connecting these two donors. Such are the prop- 
erties of the exponential function: 1 

M d x , » (, then 
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Consequently, it is amply sufficient to retain 
in the resistor network only those resistors that 
connect each donor with four or five of its nearest 
neighbors. 

The Sphere Problem Revisited 

Our next step is to calculate the resistances. We 
suggest the following line of reasoning. Let us 
disconnect all resistors supposed to join the 
balls that stand for donors and start connecting 
them in a predetermined sequence. 

First we solder in the resistances that connect 
the donors separated by a distance less than a 
certain length r\ This means that we include 
the lowest resistances smaller than M f = 
— ex P (2 r'/a§). If r' is much smaller than 
the mean separation r don between donors, then 
at this stage we connect only the rare donors 
separated by anomalously short distance. As a 
rule, these resistances are not interconnected and 
cannot sustain the flow of current through the 
system. 

Let us gradually increase r', each time adding 
new groups of resistances. Beginning with a 
certain value r' — r c , the resistors form an in- 
finite cluster. At r' > r c the system is electri- 
cally conductive and has a finite resistivity. 

The problem of finding r c is nothing less than 
the sphere problem (see Chapter 7). Indeed, let 
us construct a sphere of radius r' around each 
donor. Then resistances will be connected only 
between those donors one of which lies inside 
the sphere constructed around the other. Accord- 
ing to the results of Chapter 7, such donors 
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join into an infinite cluster when the condition 
B c = ( 4/3) JiA don rg=2.7±0.1 
is satisfied, whence 

r 0 = (0.865 ± 0.015) N J 0 n /3 « 1.39r don (2) 

Calculation of Resistivity 

As was mentioned in Chapter 12, the backbone 
of an infinite cluster can be pictured as a three- 
dimensional wire network (see Fig. 44) with knots 
spaced by distances of the order of the correla- 
tion radius. In the present problem each wire 
must be regarded as composed of numerous balls 
connected with resistors, and expression (1) of 
Chapter 12 must be rewritten in the form 

-ff ( r ') = r don ( | r'— °r c I ) ^ 

Expression (3) turns to infinity at percolation 
threshold, with the ..critical exponent v which, 
by virtue of universality, is the same as in other 
three-dimensional problems of percolation theory 
(v = 0.8-0. 9). When r' — r c is of the order of 
r don> we are ^ ar f rom ^e threshold, and the 
correlation radius becomes of the order of r don . 

When r' = r c , the separation between network 
sites in an infinite cluster grows to infinity. At 
precisely this point, the density of the infinite 
cluster is still zero. However, if r' > r c , the 
cluster forms channels going across the whole 
system and resulting in a finite resistivity. 

Let us continue with the procedure of adding 
resistances. If we switch on resistances between 
donors spaced by distances from r e to r c + ga§, 
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where g is a number smaller than unity (e.g. 
g = 0.2), the maximum one among the con- 
nected resistances will remain practically unal- 
tered because exp [(2r c /ag) + 2 g] & exip(2r c /a%). 
On the other hand, an infinite cluster represents 
a network with site-to-site separation of about 

y <4) 

In the theory of hopping conduction this network 
is referred to as critical. It is through this network 
that the electric current flows. 

Indeed, resistivity will not be appreciably 
changed by any further increase of r', that is, 
by adding the remaining resistances. As we find 
from formula (1), resistances for r — r c >> ag 
are very large in comparison with resistances for 
r = r c . In fact, practically no current flows through 
these resistances since it “prefers” the critical 
network whose maximum resistances are 

^max = #0 ex P ( “ 2r J a D 

The last step to accomplish is to calculate the 
resistivity of a critical network. By analogy to 
what we did in Chapter 12, this network will be 
represented by a wire skeleton modelling a 
simple cubic lattice (see Fig. 12) with period R c . 

Each wire connecting two neighboring sites 
of the lattice consists of a large number of resist- 
ances given by formula (1) and connected in 
series. As individual resistances differ from one 
another enormously, the resistance of one wire 
must be set equal to the maximum one of the 
resistances it comprises, that is, to The 
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cubic unit cell of a simple cubic lattice (see 
Fig. 27) consists of four wires connected in paral- 
lel, each wire belonging at the same time to 
four other cubic unit cells, so that there is one 
wire with resistance i? max per each unit cell. 

In order to find resistivity, it is necessary to 
multiply the total resistance ^? max °* a cubic 
unit cell by the area of its face ( R* ) and divide 
it by the edge length (R c ). This gives 

P = ^ m ax#c = ^o-tfcexp(2r c /aS) (5) 

We thus deal with a cubic unit cell as if it 
were filled by a homogeneous material with re- 
sistivity p. This is the meaning imparted to the 
concept of resistivity of a system that is in fact 
strongly inhomogeneous. 

Discussion of the Result 

It must be borne in mind that different elements 
of formula (5) were obtained with different de- 
grees of accuracy. The factor called the 

preexponential factor, was found to within a 
numerical coefficient. Indeed, we do not really 
know the quantity g in formula (4) for R c . Be- 
sides, the replacement of the network of an infinite 
cluster by a periodic lattice may well lead to an 
error in the numerical coefficient. It must be 
mentioned that no theory has been constructed 
yet that could reliably determine the numerical 
coefficient in the preexponential factor. 

As for the quantity in the exponent of the 
exponential expression (5), it is known to a good 
accuracy. In the true spirit of the derivation of 
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formula (5), r c in the exponent could be replaced 
by r 0 + a%g. 

Then the sum 2r c !a% - 1 - 2 g would appear in- 
stead of 2rja%. The second term of the sum ac- 
tually characterizes the uncertainty in our knowl- 
edge of the exponent. As has been mentioned, 
2 rja% » 1 at low donor concentrations, so that 
the relative error is small. In fact, it reflects the 
uncertainty in the numerical coefficient of the 
preexponential factor. 

The dependence of resistivity p on a donor 
concentration is determined solely by the expo- 
nent of the exponential. By virtue of (2) and (5), 



In p = In f 



2 - 0.86 

»ja«s 




The first term on the right-hand side of (6) 
depends on A don much weaker than the second 
term, and thus can be considered constant. A com- 
parison between formula (6) and experimental 
data, carried out for numerous semiconductors, 
demonstrated that this formula gives a very 
good fit to the In p vs N 6on dependence; this 
agreement was an important achievement of the 
theory described above. 

The method discussed in this chapter was 
successfully applied to determining the tempera- 
ture dependence of hopping electric conduction, 
and to its dependence on the external magnetic 
field and on a number of other parameters. Fur- 
thermore, this method holds for any inhomogene- 
ous system whose resistance varies in a wide range 
as a function of coordinates, 
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Chapter 14 

Final Remarks 

In this last chapter we shall briefly discuss some 
applications of percolation theory that have not 
been covered in the preceding chapters but never- 
theless appear to us as quite interesting, and we 
shall also make a summary, in an attempt to point 
at the features common for all the problems 
discussed in the book. 

Some Applications 

Flow of liquid through a maze. This problem is 
closely linked to that of gas penetration into 
carbon gas masks that served as the starting 
point for percolation theory. Imagine a porous 
body through which a liquid is forced under 
pressure. The liquid does not wet the material 
of the body so that capillary forces resist the 
penetration of the liquid into pores. 

A maze of capillaries permeates the body, but 
the capillary diameters are vastly different. They 
are wide in some places, and very narrow in 
others. 

First imagfne that the body has a single cylin- 
drical capillary. At a certain pressure produced 
by a piston 2 a liquid is on the left of the body 1 , 
and at atmospheric pressure air is on its right 
(Fig. 48). The body is fixed and cannot move. 
Because of surface tension, a nonwetting liquid 
penetrates the capillary only if its pressure exceeds 
atmospheric pressure by the quantity 2a/ R, where 
R is the capillary radius, and a is the surface 
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tension. At lower pressures, the liquid forms a 
convex meniscus but cannot pass through the 
capillary. 

Now imagine that the capillaries within the 
body 1 are of different radii. Some of them (the 
widest) let the liquid flow through at a given 
pressure, but others (narrower ones) do not. As 
pressure increases, the number of permeable 
capillaries grows. At small pressures, when only 




Fig. 48. 

the widest capillaries are accessible to the liquid, 
it cannot penetrate the body deeper than its sub- 
surface layer. However, at a certain critical 
pressure, the permeable capillaries form a system 
that permeates the whole body. Beginning with 
this pressure, the liquid can be squeezed through 
the body. 

Percolation theory calculates this pressure as 
well as some other characteristics of the process, 
which constitutes a very important practical 
application of the theory. 

Formation of polymer gels. A polymer consists 
of a very large number of elementary units (mo- 
nomers). Monomers in a solution may bind to 
one another and form a complex three-dimen- 
sional network that permeates the whole system. 
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The result is a gcZ, that is, a solid-like medium 
| resembling jelly. 

A model is known that adequately describes 
| the formation of a gel. It is essentially a problem 
| of percolation^ theory in which white sites are 
the molecules of a monomer, and black sites are 
the molecules of a solvent. The bonds between 
white sites form with a temperature-dependent 
| probability. A gel appears when an infinite cluster 
I forms out of connected white sites. 

This problem of percolation theory is said to 
be mixed because both sites and bonds are ran- 
I dom elements here. Let the probability for a site 
to be white be x 1 (it equals the concentration of 
monomer molecules), and the probability for 
a given bond to be intact be x 2 * We want to find 
the range of values of and x 2 in which an in* 
finite cluster of interconnected white sites is 
formed. 

By definition, x x and x 2 vary within the inter- 
val from zero to unity. 

If x 2 = 1, that is, if no bonds are broken, an 
infinite cluster exists for all x t in the range 
j % <* 1 J where x s is percolation threshold 
of idie site problem. If x x = 1, that is, if all 
sites are white, the condition under which an 
infinite cluster is formed is x b < x 2 < 1, where 
x . is percolation threshold of the bond problem. 

The square in Fig. 49 is the domain of the 
variables x 1 afid x 2 . The solid curve is the graph 
of the function x mln (x 2 ) that describes the bound- 
ary of the domain in which an infinite cluster 
exists. The function z m ,„ (*,) gives the minimum 
value of x x for each value of x 2 within the interval 
x b < x 2 < 1, at which the infinite cluster exists. 
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It is readily understood that a; mln (1) = x s , and 
^min (^b) = 1. The domain of the infinite cluster 
in Fig. 49 is the hatched area. 




Fig. 49. 
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in a well-defined manner (e.g. by a random-number 
generator with known properties). 

The common element of percolation theory 
problems is that all of them have identical geo- 
metry of bonded elements in the vicinity ol 
percolation threshold. In order to recognize this, 
it is necessary to disregard the small-sca e 
structure depending on the type of bonding an 
the properties of elements, and monitor only the 
connectivity of large blocks. The universal large- 
scale geometry dictates the universal properties 
of physical quantities that depend on the struc- 
ture of large clusters. This is the feature that 
unites percolation theory problems that look so 
dissimilar at first glance. 



If the probability x 2 is known for each value 
of temperature, the function ^ mln (x 2 ) allows 
us to determine the range of temperatures and 
monomer concentrations in which a gel is formed. 



What Is Percolation Theory, After All? 

So far we have carefully avoided defining per- 
colation theory. This would be far from simple. 
But having come to the last page of the book, we 
shall try to recapitulate the common element 
of all the problems outlined in the book, and 
formulate the subject of percolation theory. 

Percolation theory deals with the connectivity 
of a very large (macroscopic) number of elements 
under the condition that the bonding of each 
element to its neighbors is random but prescribed 
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Chapter 1 



. a lf i general rule for calculating the mean value 
states that each possible value of a random variable must 
he multiplied by the probability of this value, and the 
products be summed up. The probability for each of the 
taces of the cube to be on top equals 1/6. Therefore, 

a = ( 1 / 6) -1 + ( 1 / 6 ) .2 + ( 1 / 6) -3 + ( 1 / 6 ) - 4 + ( 1 / 6) -5 
+ ( 1 / 6). 6 = 21/6 T 

2. The results of individual trials change, but the 
average value x c (Jf) calculated over many trials re- 
mains the same because on the average the rightward mo- 

w,lh same probability as the downward 
motion. Correspondingly, x c = lim * c (Jf) will not 

change either. " V_>tx> 

3. Let us denote the threshold values recorded in the 
ith trial by xi and x h where x'\ corresponds to the new 

e mtion, and X( to the old definition (rightward percola- 
tion) U a decrease in x first interrupted the rightward 
percolation and then the downward percolation, then 

Xi ' Bllt ^ Ihe sequence was reverse, then x[ > x u 
the average value of percolation threshold of the new 
definition is 






X l~\~ X 2 + • • • +&Q 

Q 



and that of the old definition is 
x c ( JT ) = 3 + ^2 + — +*q 

The total number of trials, Q, is considered very large in 



both formulas. However, situations with x[ > xi will 
inevitably realize if the number of trials is large. For this 

reason, x i ( JIT ) > * c (JT)- Nevertheless, 
lim x' c (,_#*) = lim x c (J^) = x c 

Jf -V oo Jf -> OO 

The thing is that in an infinite system percolation 
threshold is not a random variable but a certain quantity, 
not varying from one trial to another. At the same time, 
differences between percolation thresholds for different 
directions represent a random phenomenon. From the 
standpoint of percolation, all directions are on the 
average equivalent. Hence, x c is independent of 
direction. 

4. The solution is quite similar to that of Problem o. 
Let us denote the threshold values obtained in the 

trial by xi and x t , with x[ corresponding to the new 
definition of threshold, and ^ to the old definition. It can 
be readily proved that xi < x t . Arguments similar to 
those of the preceding problem give 

x c {jD<x 0 (jP), but lim x" c (JT) 

JT-+OC 

= lim x c (jp^)~ x 0 
jr ~>oo 

5. Formula (8) yields 6 = 0.01. This signifies that 

“typical” deviations from the mean value are ±0.01. 
Hence, the last numeral in the estimate 0.59 reported by 
Watson and Leath was very likely in error. The proba- 
bility of error in the first numeral after the decimal point 
is much smaller. Since only one trial was run, the 
authors of the paper could only evaluate the error with 
which they determined percolation threshold on the 
utilized sequence of blocked sites (it proved to be ±0.005). 
However, they could not say anything about the possible 
change in the result if the experiment were rerun with 
a different random sequence. Later studies in which many 
trials were run at one value of Jf , and studies with large 
1 rtf' in TVlOGO etllHlVc! H 1 Clfl fl P- 
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monstrated that even the second decimal place is correct 
in the number 0.59; to a large extent, this must be re- 
garded as pure “luck”. 



Chapter 2 

1. The fraction of blocked sites is 1 — x— JIT'/JIT* 
and the fraction of nonblocked sites is x = (JT ~JT')I JT . 
If Q sites are chosen at random, Qx of them will be non- 
blocked, and Q (1 — x) will be blocked (the greater is Q , 
the better the accuracy is with which this relation holds). 
Consequently, the probability for a randomly selected 
site to be blocked is Q (1 — x)/Q = 1 — x , while the 
probability for it to be nonblocked equals Qx/Q = x. 
The sum of the probabilities equals unity because a 
site^can be either blocked or nonblocked: 1 — x + x = 

2. The probability of any sequence of three fixed num- 
bers equals 1/64/6 -1/6 = 1/216. The number of differ- 
ent sequences satisfying the formulated conditions for 
the numbers 1, 2, 3 equals 6 (123, 213, 321, 231, 132, 312), 
and for the numbers 1, 2, 2 it is 3 (122, 212, 221). The 
probability for one of the possible sequences to realize 
equals the sum of these probabilities. Therefore, in the 
first case the sought probability is 6*1/216 = 1/36, and 
in the second it is 3*1/216 = 1/72. 

3. (0.8) 3 *(0.9) 4 = 0.336. 

4. The solution is left to the reader. 

5. The distribution function of a random variable a 
is constant within the interval (—1, 1) and equals zero 




Fig. 50. 
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everywhere else (Fig. 50). The total area of the rectangle 
bounded by the curve / (y) (in the present case by a hori- 
zontal straight line), the abscissa axis, and the perpen- 
diculars erected at points — 1 and 1 must equal unity, 
consequently, / (?) = 1/2 for -1 < y < 1. By virtue 
of the general rule, the sought probability equals the 
area of the rectangle bounded by the line / (y), the abs- 
cissa axis, and the perpendiculars erected at points y ' — 
= —3/4 and y = —1/4 (the rectangle is shaded m Fig. 50). 
The probability equals 

l — 1/4 - (-3/4)] *1/2 = 1/4 

6. The variable y assumes all values from — oo to + oo. 

Consequently, we have to set A = — oo B = +°° 
in formulas (3) and (4). According to formula (3), 



oo 

a = ( yfjr(y)dy 



Substituting here formula (6), we find 



oo 

5 sexp ( 



* 



The integrand contains an odd function. Substituting 
y — —t and comparing the result with the original tor- 

mula, we note that a. — hence, a 0. 

According to formula (4), the variance is 



OO 

j y 2 f.jir (y)dy 



-6> J » ! e.,p(— * 



Substituting y 



,1, we obtain 
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62 = W 6 > i * 2e x P (~t*)dt 

— oo 

The integral over t equals ]/jx/2, so that 8 2 = 6 2 



Chapter 3 

1. By definition, P (1) = 1 . If * is near i y unity 

the sites may not belong to an infinite cluster for two 
reasons: 

(0 They may contain nonmagnetic atoms. The frac- 
tion of such sites is 1 — x. 

(ii) Magnetic atoms may be isolated from an infinite 
cluster as, for example, atom B in Fig. 9. But if x is 
close to unity, so that nonmagnetic atoms are few, this 
or is not so important because such isolation requires 

Zl n r apet p,f 0m f to gather around one atom 
Crour in the case of the plane lattice shown in Fig. 9) 
Ihe probability of such an event is low if the number of 
nonmagnetic atoms is small. The second factor can thus 
be safely neglected and we can assume that the fraction of 
atoms belonging to the infinite cluster simply equals 
tbe traction of magnetic atoms. We thus have P (x) = * 
provided 1 — x < 1. v } 

„ciaUW,l Ch i ato . m , of ? s * tn ple cubic lattice has six nearest 
neighbors located along the cube’s edges (see Fig. 12). 
The probability W 0 for all the nearest neighbors of an 
atom to be nonmagnetic equals the product of six proba- 

one^ofth ^ t~ (1 k~ X) °' The P robabilit y W of at least 
one ol the atoms being a magnetic atom equals 

W (x) = 1 — w a = 1 — (1 — x ) e 

By virtue of formula (2) of Chapter 3, 

p 2 (*) = xW (a;) = x [1 — (1 — x) d ] 

so that if x < 1, 

Pi (x) & 6x 2 



Answers and Solutions 



247 




Obviously, 

P 2 (x) — x [1 — (1 — x) z ] 

for any lattice in which each atom has z nearest neighbors, 
so that for x <C 1 

P 2 (x) = zx 2 

3. Figure 51 shows 12 atoms in the neighborhood of 
atom 0. All of them can take part in forming a three - 
atom cluster. Such a cluster can be composed of, for 
instance, atoms 7, 0, 2 if all these three atoms happen 
to be magnetic. The probability of this event equals the 
product of three probabilities: x*x*x = x 3 . The proba- 
bility for the cluster to be formed of atoms 0, 4, 12, or 
of any other trio of atoms, also equals x 3 . 

First we have to answer the question about the number 
of such groups of three atoms. We begin with counting 



10 




Fig. 51. 



how many of these groups include atoms 0 and 1. There 
are six of them: 015 , 018 , 019, 103 , 102 , 104. 

Now we consider three-atom groups including atoms 0 
and 3 , but not atom 1. There are five of them: 036, 
03 11 , 037 , 203, 304. 

Likewise, there are four three-atom groups with atoms 
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0 and 2 , but without atoms 1 and 3: 025 , 026 , 02 10, 
204, and three-atom groups with atoms 0 and 4, but with- 
out atoms 1, 2, 3. These are 047, 048, 04 12. 

There are thus 6 + 5 + 4-f-3 = 18 three-atom 
groups, each having the probability x 3 . We want to find 
the probability for at least one of them to realize. If 
x < 1, the events consisting in the formation of any one 
of the three-atom groups can be treated as incompatible. 
Indeed, the probability for the groups 102 and 015 to 
form simultaneously equals the probability of four atoms, 
0125, being magnetic, that is, the probability of forma- 
tion of a four-atom cluster. The probability of this event 
equals x A = x 3 *x x 3 . Therefore, the incompatibility 
of these events is confirmed for x <c 1. The probability 
for at least one of the three-atom groups to form is then 
equal to the sum of their probabilities, and 

P 3 = 18* 3 

4. The probability for a randomly selected atom to 
belong to a cluster of not less than two atoms can be re- 
presented as the sum of the probabilities of incompatible 
events: 

P % (x) — P3 {x) -f- P 2 (x) (1) 

where P 2 (x) is the probability for an atom to belong to 
a cluster of two atoms. Therefore, 

P 3 {x)^P 2 {x)-P 2 {x) (2) 

The function P 2 (x) is given by formula (3) of Chapter 3, 
and hence, we only need to calculate P 2 (x). 

Atom 0 (see Fig. 51) can form a two-atom cluster 
with atoms 1, 2, 3 or 4. The probability of the cluster 
being formed of atoms 0 and 1 equals the probability for 
both these atoms to be magnetic times the probability 
for atoms 2, 3, 4, 8, 9, and 5 to be nonmagnetic, that is, it 
equals x 2 (1 — x) Q . The events consisting in the two- 
atom cluster being formed of atoms 0 and 2, or 0 and 3, 
or 0 and 4 are essentially the same. All these events are 
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incompatible, and thus the probability P 2 (#) equals the 
sum of lour probabilities 

P 2 (z) = 4a; 2 (1 — x) G (3) 

Substituting formula (3) into formula (2), we obtain 

P 3 (x) = x [1 — (1 — a;) 4 ] — 4a; 2 (1 — x) G (4) 

which is the solution of the problem. 

Making use of the binomial theorem, we can easily 
show that expression (4) contains no terms with power 
below 3. If x < 1, then P 3 {x) « 18a: 3 , in accord with 
the result of the preceding exercise. 



Chapter 4 

1. 0.0085, 0.0072, 0.0051, 0.0026, 0.0006, 0.0000, 
0.0000...; 0.0067, 0.0044, 0.0019, 0.0003, 0.0000, 0.0000...; 
0.0032, 0.001, 0.0001, 0.0000, 0.0000.... 

2. The number b is in fact an rc-digit number. Con- 
sequently, b < 10 n . In order to generate the next number 
b' , we have to findfr 2 , divide it by 10 n , and take its inte- 
gral part. Therefore, b ' < b 2 /10 n . But5 2 /10 n = b-(b/i0 n ). 
Furthermore, {b 2 / 10 n ) < b because (b/10 n ) < 1. We thus 
conclude that b' < b, which was to be proved. 

3. 5, 15, 5, 15, 5, 15.... 

4. 5, 16, 9, 8, 5, 16, 9, 8, 5.... 

5. 5, 17, 13, 1, 5, 17, 13, 1, 5.... 

6. Condition (c) is nowhere satisfied. In Exercise 3, 
c = 0, in Exercise 5, condition (a) is not satisfied either, 
and so forth. 

7. Let X 0 = 0. This gives us 0, 3, 1, 4, 2, 0, 3, 1, 4, 
2. . ., and X x coincides with one of the numbers of 
this sequence whatever X 0 is. 

8. The distribution function of random numbers ob- 
tained from the generator is shown in Fig 4. These num- 
bers form the array V. The fraction of nonblocked sites 
in the array V equals the fraction of random numbers in 
the array V , which satisfy the inequality V < t. There- 
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fore, the mean fraction of nonblocked sites equals the 
probability for a random number to be less than t . By de- 
finition of distribution function, this probability equals 
the area bounded by the curve / ( y ), the abscissa axis, 
and two perpendiculars erected at points 0 and t. In the 
present case this is the area of a rectangle, equal 
to t. Hence, the mean fraction of nonblocked sites x 
equals t . 



Chapter 5 

1. When x is close to unity, almost all sites belong to 
an infinite cluster. Only those sites stay out of this clu- 
ster whose bonds connecting them with the surrounding 
system are all broken. The probability that one specific 
bond is broken equals 1 — x. If the lattice is square, it 
is necessary to interrupt four bonds originating from 
a site in order to isolate this site from the system (Fig. 52a). 
The probability of this event equals the product of the 
appropriate probabilities, that is, it equals (1 — a;) 4 . 
In order to isolate two sites from the system, six bonds 
must be interrupted (Fig. 526). The probability of this 
event is (1 — x) 6 . If (1 — #)<1, this probability is 
much smaller than that of one isolated site. Therefore, 
we can assume in the limiting case which is of interest 
now that all isolated sites are single, and the probability 



I 




Fig. 52. (a) One isolated site; (6) two isolated sites. The 
intact bond is shown by the solid line, and the broken 
bonds by the dashed lines. 



for a randomly selected site to be isolated is (1 — x ) 4 . 
The probability for a randomly selected site to be non- 

isnlated pctiipIs A — (A — that, is in a smiarp lat.t.ir.p 
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it is 

P h (x) = 1 — (1 — x) 4 

Similar arguments show that in a triangular lattice 

P b (x) = 1 — (1 — x) 6 
and in a honeycomb lattice 

P h (x) = 1 — (1 — x) 3 

These results hold if (1 — a:)<l. 

2. The problem is solved by analogy to the preceding 
one. If (1 — x) < 1, almost all sites belong to an infinite 
cluster. A site chosen at random is isolated from the in- 
finite cluster if all its nearest neighbors are nonmagnetic 
atoms (for definiteness, we use the terminology of the 
ferromagnetic problem). As in the preceding problem, 
the probability for a site to have for one of its neighbors 
a magnetic atom isolated from the infinite cluster is low. 
Therefore, it is sufficient to calculate the probability 
for all neighbors of the given site to be nonmagnetic. The 
probability that a site is occupied by a nonmagnetic 
atom is 1 — x. The number of nearest neighbors equals 
the number of bonds that originate from a given site. 
Consequently, the results are the same as in Problem 1: 

in a square lattice: P s (#) = 1 — (1 — #) 4 
in a triangular lattice: P s (x) = 1 — (1 — x) % 
in a honeycomb lattice: P s (a:) = 1 — (1 — x) 3 

Hence, P s (x) = P h (x) for (1 — x) < 1, in agreement 
with formula (2). 

3. Consider a honeycomb lattice with a fraction x of 
white bonds, and a triangular lattice with a fraction y 
of white bonds. (We remind the reader that the term 
“white bond” is synonymous to “unbroken bond”, and 
“black bond” is synonymous to “broken bond”. The phrase 
“a site is connected with another site” must be interpret- 
ed, unless expressly stated otherwise, as equivalent to 
“a site connected through unbroken, i.e. white, bonds”.) 
Let us superpose these lattices as shown in Fig. 53. Such 

sit , ps as A r r O TO n n m m r\ n fnr tho lattinoa nrliiln 
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sites such as D belong only to the honeycomb lattice. 
The idea of the arguments that follow lies in that the 
percolation problem on a honeycomb lattice is reducible 
to the corresponding problem on a triangular lattice. 




Fig. 53. Star-delta transformation. The dashed lines show 
the honeycomb lattice, and the solid lines show the trian- 
gular lattice. 

The probabilities for sites A, B , C to be connected 
with one another must be expressed in terms of the frac- 
tion x of white bonds in the honeycomb lattice. We need 
to make use of the geometry and statistical properties 
of the bonds that originate from D-type sites. After this 
we can discuss only the triangular lattice shown in Fig. 53 
by the solid lines and forget that each triangle has the 
dashed lines and D -type sites. 

This operation, known as the star-delta transforma- 
tion, is frequently used in designing electric circuits. 

In fact, we shall need the following quantities: 

(1) is the probability for site A not to be connect- 
ed with either site B or site C. This probability equals 
the sum of the probabilities of two incompatible events. 
The first event consists in bond AD being black, with 
arbitrary colors of bonds 1 3D and DC. Its probability 
equals 1 — x. The second event consists in bond AD 
being white, and both bonds BD and DC being black. 

aKil { tAT r\f +Uic ouont ommle thp rirndnct of three 
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(2) 1F+- is the probability for site A to be connected 
with B but not with C. It equals the probability of bonds 
AD and DB being white, and bond DC being black, and 
is found as the product of the probabilities of all three 
events: 



W+- = x 2 (1 — x) ( 2 ) 

(3) W 4. is the probability for site A to be connected 
with C but not with B. It is readily found that = 
= W+-. 

(4) W|+ is the probability for site A to be connected 
with both B and C . It equals the probability for all the 
three bonds AD , DB , DC to be white and is calculated as 
the product of their probabilities: 

W + + (x) = z 3 ( 3 > 

With these four probabilities known, we can leave 
aside the honeycomb lattice and solve our percolation 
problem on the triangular lattice. If this problem is solved, 
we shall have the critical value x h (H) for the honeycomb 
l&tticG 

Of course, this transformation does not in the least 
facilitate the solution. However, the same probabilities 
W can be expressed in terms of the fraction y (it repre- 
sents the fraction of white bonds on the triangular lattice). 
At percolation threshold these probabilities assume quite 
definite values, as yet unknown; however, equating the 
probabilities expressed as functions of x to the probabil- 
ities as functions of y % we can derive a relation between 
the thresholds of the honeycomb (*b (H)) and triangular 
fob (T)) lattices. 

Our next problem is thus to derive all four probabili- 
ties as functions of y. 

(1) kF__. For site A to be unconnected with both v 
and C, the color of bond BC is immaterial, but bonds AB 
and A C must be black. The probability of this event equals 
the product of probabilities: 



W— (a 



1 - y? 



( 4 ) 
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(2) W+_. Site A is connected with B but unconnected 
with C only if bond AB is white and both bonds BC and 
AC are black: 

W+„ = y( 1 - y) 2 (5) 

(If e.g. bond BC were white, site A would connect with C 
via path ABC.) 

(3) = W+ as in the preceding case. 

(4) W+ 4-. The probability for site A to be connected 
with both B and C equals the sum of the probabilities of 
four incompatible events. The first event consists in all 
the three bonds AB, BC, AC being white. Its probability 
is if. The other three events consist in only one among 
the three bonds being black. For instance, if the black 
bond is AB, then site A is connected with C via a white 
bond AC, and with B via path ACB. The probability of 
each of the three events equals y 2 (1 — y). Finally, we 
obtain 

W ++ {y)=-y 3 -\-3y 2 (i — y) (6) 

At percolation threshold all W (a;) must be equal to 
the respective W (y). This yields a system of equations 

(x)^ W^_(y); 1 — x-{-x(l — x) 2 = (i — y) 2 (7) 
W+_(x) = W + _(y); x 2 (i — x) = y(i — y) 2 (8) 

1V++ (x) = W ++ (y)\ x 3 = y 3 + 3y 2 (1 — y) (9) 

These equations must be satisfied by x = x\> (H) and 
y ~ x b (T). Furthermore, percolation thresholds satisfy 
relation (23) of Chapter 5, which yields that x b (H) = 
= 1 - *b (T). 

The substitution of 

V = x b (T), x = 1 — zb (T) 

into Eqs. (7), (8), and (9) transforms (8) into an identity, 

(1 - z b (T))^ * b (T) = * b (T) (1 - * b (T))» 

and (7) and (9) into the identical cubic equation: 

xi5(T)-3x b (T) + l = 0 

This equation has a single root within the interval 0 ^ 
m <: i • 
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Xb (T) = 2 sin (it/18) « 0.347 296 
Correspondingly, 

x h (H) = 1 — * b (T) = 1 — 2 sin (n/18) « 0.652 704 

4. The first step is to find the area per site for each 
lattice shown in Fig. 16, under the condition that the 
spacing between nearest neighbors is a. 

Square lattice. Each unit cell includes four sites, but 
each of these sites belongs to four adjacent unit cells. 
Hence, there is one site per each square unit cell; in other 
words, the area per site is that of one square, that is, 

S (S) = a 2 . # 

Triangular lattice. Each unit cell includes three sites, 
but each of these sites belongs to six adjacent unit cells. 
Hence, we find half a site per each triangular unit cell, 
and the area of two triangles per each site. The area ol 
an equilateral triangle with side a is Yda 2 lA. Therefore, 

the area per site is S (T) = Yda 2 /2. ... 

Honeycomb lattice. Each hexagon includes six sites, 
but each of these sites belongs to three adjacent unit cells. 
Hence, there is one hexagon per two sites. The area ot 
a hexagon equals six times the area of an equilateral tri- 
angle with side «, that is, equals 3/3a 2 /2. Hence, the 

area per each site is S (H) = 3 ] fZa 2 /k. . 

The period a is determined in each lattice via a pre- 
determined function a (x), with percolation threshold x b 
of the appropriate lattice used for x. This gives 

S( T) =(/3/2) [o(0.35)]* 

S(S) = [a (0.5)] 2 

5(H) =(3/3/4) [a (0.65)] 2 

Obviously, the function a ( x ) decreases with increasing 
x. (If the trees infect each other easier, the separation 
within a pair correspondingly decreases.) It is readily 
noticed, however, that this proposition is not sufficient 
in itself to derive even one of the above-written inequa- 
lities for areas. Indeed, in the honeycomb lattice the period 
a is the shortest, but the numerical coefficient in the 
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expression S (H) for area is the greatest, while in the 
langular lattice the situation is reverse. The lattice 
with minimum area thus cannot be found unless we know 
tne function a (x) in greater detail. 



Chapter 6 

TJ\ e filling factors for plane lattices are easily 
calculated from the results of Exercise 4 of Chapter 5* 

lYtticB aV w lllU n tr n te l thi f in the text for the honeycomb 
lattice. We shall, therefore, restrict the analysis to 3D 

^eckf to° the reader. tW ° ° f them ’ loaVing the remainin S 
Simple cubic lattice. Each cubic unit cell (see Fie 121 
contains eight sites, but each of the sites is shared by 
eight adjacent unit cells. Hence, there is one site per 
one unit cell, and the volume per each site equals the 

IrounH 6 fk, h -t Unit C , e11 ’ L , e ‘ a3 - The spheres drawn 

around the sites have a radius all and a volume 4na 3 /24 

thp Jffn lfV volume occupied by the spheres equals 
the ratio of the volume of one sphere to the volume per 
site. Hence, m a simple cubic lattice 1 

/ (SC) = 4ji/24 « 0.52 

(see B F?/'?SM^n d i CUbiC ' at . t ‘ Ce - E , ach cubic unit cell 
vert.v«?’nf 8 f2 co tains nine sites, eight of which are at the 
vertices of the cube and one is at the center point. The 

lv Jfc-nT ce “ te , r belongs to one cubic unit cell exclusive- 
° f th( :, vertex sites belongs to eight adjacent 
A ells - Hence, there are two sites per one unit cell, 
and the volume per each site is half the volume of the 
cubic unit cell, i.e. a 3 / 2. The nearest neighbor of each 
site is at a distance of half the body diagonal of the unit 
cell, i.e. V3a/2. The radius of a sphere drawn around each 
site is y 3a/4 and the volume of the sphere is (4 jt *3 3 / 2 /3 • 
r (n ^6) a*. The filling factor equals the ra- 
so°that 6 V °* Ume s Pb ere to the volume per site, 



/ (BCC) = nY 3/8 
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Chapter 7 

1. The first coordination group of site 0 consists of 12 
type-i sites (Fig. 54), the second consists of six type-2 




Fig. 54. Neighbors of site 0 in the fee lattice. 

sites, and the third of 24 type-3 sites (not all of them are 
shown in the figure); 12 -J- 6 -J- 24 = 42. 



Chapter 8 

1. From formulas (3) and (5) we find = 636 'KH cm 
= 636 A, N 0 = 7.8 -lO 13 cm- 3 . This is a very low 
critical concentration. The transition to metallic 
conduction takes place when there is one. impurity 
atom per 10® atoms of the semiconductor matrix! A very 
sophisticated technique of purification is required to ma- 
nufacture semiconductors with such a low impurity con- 
centration. 



Chapter 11 

1. Let us introduce a function Q (r), defining it as 
fV, r*r*ni»nhi li tv fm* mi information that has appeared at a 
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randomly chosen site to reach only a finite number of 
other sites. As before, P (x) = 1 - Q. The number of 
independent channels through which the information 
emerges from each site is q. Let us determine the proba- 
bility for one of these channels to he interrupted at some 
step. The interruption can result from one of two incom- 
patible events, a or b: a consists in the first bond of the 
channel being broken, and b consists in the first bond being 
intact, but the site to which this bond leads being capable 
of passing the information to only a finite number of other 
sites. The probability of event a is 1 — *, and that of 
even j j 1S i ( x )• T ^ e probabilities of incompatible events 

are added. Consequently, the probability for one channel 
to be interrupted equals 1 - * + xQ (x). All the channels 
being independent, the probability for all of them to be 
interrupted is [1 — x -f xQ (x)]i. 

This gives an equation for Q (x): 

Q ( x ) = [1 — x - I xQ (,r)]7 

This equation reduces to Eq. (4) by the substitution 
v =- 1 — x xQ. Finally, we obtain 

P — 2 ? 2 

Note that percolation thresholds in the site and bond 
problems on the Bethe lattice are identical (x c = 1/q) 
In fact, this could be predicted at the start. Indeed, let us 
assume that the site to which the broken bond leads is 
blocked. This allows us to assume that all bonds are in- 
tact, so that our bond problem has been transformed to 
the site problem. Hence, these problems necessarily have 
identical percolation thresholds. 



Chapter 12 

1. We need to calculate the resistance of a cube with 
unit-long edge. The number of wires connected in parallel 
is, as before, l//? 2 , but the length of each wire does not 
any more equal one unit of length. It is longer by a factor 
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equal to the ratio (X/R) — (x — x c ) C v > Correspond- 
ingly, the resistance of one wire is not p 0 but p 0 (x — 
— x c )~^~ v \ The necessary result is therefore obtained 
by replacing the quantity p 0 in formula (6) for o 3 by 
p 0 (x — x c ) ~&~ v \ The substitution gives 

o = a 3 {x — a: c ) t+v 

where o 3 = po^Z -2 . Hence, t — £ -[- v. 
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